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ABSTRACT 



FINeSSE: Fermilab Intense Neutrino Scattering Scintillator 

Experiment 

Understanding the quark and gluon substructure of the nucleon has been a prime 
goal of both nuclear and particle physics for more than thirty years and has led to 
much of the progress in strong interaction physics. Still the flavor dependence of the 
nucleon's spin is a significant fundamental question that is not understood. Exper- 
iments measuring the spin content of the nucleon have reported conflicting results 
on the amount of nucleon spin carried by strange quarks. Quasi-elastic neutrino 
scattering, observed using a novel detection technique, provides a theoretically clean 
measure of this quantity. 

The optimum neutrino beam energy needed to measure the strange spin of the 
nucleon is 1 GeV. This is also an ideal energy to search for neutrino oscillations at high 
Am 2 in an astrophysically interesting region. Models of the r-process in supernovae 
which include high-mass sterile neutrinos may explain the abundance of neutron- 
rich heavy metals in the universe. These high-mass sterile neutrinos are outside the 
sensitivity region of any previous neutrino oscillation experiments. 

The Booster neutrino beamline at Fermilab provides the world's highest inten- 
sity neutrino beam in the 0.5-1.0 GeV energy range, a range ideal for both of these 
measurements. A small detector located upstream of the MiniBooNE detector, 100 m 
from the recently commissioned Booster neutrino source, could definitively measure 
the strange quark contribution to the nucleon spin. This detector, in conjunction 
with the MiniBooNE detector, could also investigate disappearance in a currently 
unexplored, cosmologically interesting region. 
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Chapter 1 
Introduction 

The Fermilab Intense Neutrino Scattering Scintillator Experiment ("FINeSSE") is 
designed to measure the strange quark contribution to the spin of the proton, and to 
search, in conjunction with the MiniBooNE experiment, for disappearance. These 
measurements will employ a novel detection technique, and will examine a kinematic 
region inaccessible to any existing or presently-planned experiment. FINeSSE will 
be located 100 m from the Booster neutrino beamline production target, and 441 m 
upstream of the currently-running MiniBooNE experiment. The low energy Booster 
neutrino beam is crucial to achieving this experiment's goals; they can only be re- 
alized on this Fermilab beamline. The number of protons on target (POT) needed 
to reach the FINeSSE physics goals is 6xl0 20 , attainable in two years of running. 
The detector is designed to resolve both short, low energy proton tracks, and longer 
muon tracks from vN interactions. The FINeSSE detector will cost $2.25M ($2.8M 
with contingency); the FINeSSE Detector enclosure, $800K ($1.6M with contingency, 
escalation, and EDIA). 

1.1 Outline 

This proposal sets forth the experiment's goals, design, costs, and schedule in the 
following sections: 

• Chapter 121 provides the physics motivation for FINeSSE; 
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• Chapter 01 describes the flux and event rate at FINeSSE produced by the Booster 
neutrino beamline; 

• Chapter |U details the detector design, construction, and installation, as well as 
the readout and trigger systems; 

• Chapter El examines event kinematics, efficiencies, and backgrounds for the FI- 
NeSSE physics measurements; 

• ChapterElpoints out additional physics measurements the FINeSSE experiment 
can perform; 

• Chapter [7| provides a breakdown of costs for the detector, the electronics, and 
an enclosure for both, as well as a timeline to first beam in mid-2006. 

1.2 FINeSSE Physics, Detector, and Neutrino Beam 

The fundamental question of the flavor dependence of the spin of the proton is not 
understood. The still unresolved "spin crisis" [T| points to the fact that the proton's 
spin is not carried, as was expected, by the valence quarks. How much is carried by the 
light quark sea has been the subject of much controversy. In addition, measurements 
of the spin carried by the strange quarks in the nucleon have been plagued by model 
assumptions and experimental limitations. The FINeSSE experiment will measure 
the proton's strange spin, As, avoiding the pitfalls of previous measurements; our 
approach will be described in detail in Chapter |21 

FINeSSE, in conjunction with the MiniBooNE experiment, is sensitive to dis- 
appearance in an as-yet unexplored, astrophysically interesting region. Incorporating 
oscillations to 1 eV sterile neutrinos into the r-process in supernovae can explain the 
abundance of neutron-rich heavy metals in the universe [2j. Oscillations between these 
sterile neutrinos and muon neutrinos are expected over short baselines for neutrino 
energies around 1 GeV. The combination of FINeSSE and MiniBooNE, functioning 
as near and far detector, enables a disappearance search sensitive to these os- 
cillations. This sensitivity exceeds that of any existing or planned experiment, and 
permits exploration of the full allowed 3+1 region. 

The physics goals of FINeSSE can be achieved using a novel, relatively small, 
tracking detector placed 100 m from the neutrino production target on the Booster 
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Main Injector 




Figure 1.1: Proposed FINeSSE location with respect to the existing neutrino target 
building (MI- 12) and MiniBooNE detector. North is to the right in this drawing. 

neutrino beam line (Fig. Il.ljl . The detector is comprised of two subdetectors. The 
upstream Vertex Detector is a highly-segmented, liquid scintillator "bubble chamber" 
that tracks particle interactions via scintillation light read out on a grid of Wavelength 
Shifting (WLS) fiber strung throughout the volume. The Muon Rangestack down- 
stream of the Vertex Detector ranges out high energy muons produced in neutrino 
interactions. 

The Booster neutrino beam provides an intense source of muon neutrinos (with a 
small background of electron neutrinos) in the energy range of 0.5-1.5 GeV and a mean 
energy of ~ 700 MeV. This spectrum is ideal for the vp elastic scattering measurement 
as well as for the disappearance search. Using the currently estimated MiniBooNE 
neutrino flux jS], and assuming 3.0 x 10 20 protons on target per year jl], there would 
be approximately 360k neutrino scattering events in a detector of 9 ton active volume 
during the FINeSSE run of 2 years. This would provide a neutrino event sample of 
unprecedented size in this energy range. 



1.3 The FINeSSE Collaboration 



FINeSSE is currently a collaboration of 29 scientists from six universities and two 
national laboratories. Once approved, the collaboration is expected to grow. Both 
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university groups and national labs have already incorporated post-docs, graduate 
students, undergraduates, and a high school teacher into FINeSSE physics projects. 
FINeSSE has organized an executive committee comprised of senior scientists from 
each group to oversee the project. The collaboration is diverse yet balanced in a 
number of ways: 

• Collaborators are drawn from both the nuclear and particle physics communi- 
ties. This will help the collaboration to attain its physics goals, which span 
both of these communities, and will encourage cross- disciplinary interactions 
between nuclear and particle physicists, as well. 

• FINeSSE is comprised of both MiniBooNE and non-MiniBooNE scientists. This 
balance helps to ensure a good understanding of a neutrino beam already well 
studied by MiniBooNE, and at the same time bring in new ideas and perspec- 
tives on physics with this beam. 

• FINeSSE scientists carry with them a broad spectrum of experience. The new 
perspectives brought by the FINeSSE spokespeople are tempered by a well- 
seasoned executive committee comprised of experienced scientists with, again, 
both nuclear and particle physics backgrounds. 

FINeSSE's diversity in background and experience are already a great asset to 
its physics program, and will help to sustain it throughout its physics run. 

1.4 FINeSSE as part of the Fermilab program 

FINeSSE brings timely and important physics to the Fermilab program. In addi- 
tion, FINeSSE takes advantage of the investment Fermilab has already made in the 
Booster neutrino beamline, provides an excellent training ground for young Fermilab 
scientists, and already actively contributes to the growing Fermilab neutrino program. 

To achieve its physics goals, FINeSSE requires 6.0 x 10 20 total protons on target, 
received over the course of a two year run on the Booster neutrino beamline. This is 
within the Booster's capability in the era of NuMI and Run II running, as described in 
Chapter El and Reference [3]. It also takes advantage of a running beam, and another 
running experiment, adding value to already committed resources. 
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FINeSSE is a small, focused collaboration. Such groups are proven training 
grounds for graduate students and post-docs, vital to the future of Fermilab and 
of high energy physics. It is on such small-scale experiments that young scientists are 
guaranteed to get their hands on almost every aspect of design, construction, data 
taking, and data analysis. 

In summary, FINeSSE represents an important addition to Fermilab's program: 
it provides an extraordinary opportunity for physicists from a number of subfields 
and a variety of levels of experience to work together; it makes advantageous use of 
an existing beamline; it increases the physics reach of an existing experiment; and it 
uses a novel detection technique to address significant and interesting physics. 



1.5 Requests to the PAC 

Please consider the following specific requests with respect to approval and funding 
for FINeSSE: 

• Grant this experiment "stage 0" or "stage 1" approval at this time (approval 
pending response to any outstanding questions). This will allow us to submit 
an NSF proposal by a January 2004 deadline. 

• Recommend to the Fermilab directorate to support FINeSSE for the first stages 
of detector enclosure design work. 

• Recommend to the Fermilab directorate to provide FINeSSE with office and lab 
space. 
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Chapter 2 

Physics Motivation 



Two physics measurements form the foundation of the FINeSSE program: the mea- 
surement of the strange spin of the proton, As; and the search for disappearance 
in an astrophysically interesting region. Both topics are compelling, and can only 
be addressed with the Booster neutrino beam design. Along with these studies, a 
complement of other measurements and searches are open to FINeSSE. These other 
physics projects are addressed in Chapter El In this chapter we concentrate on the 
two main physics goals, which make FINeSSE unique. 

2.1 Strange Quark Contribution to Nucleon Spin 

From the time that the composite nature of the proton was discovered, physicists 
have sought to understand its constituents. The study of nucleon spin has grown into 
an industry experimentally, and opened new frontiers theoretically. Deep Inelastic 
Scattering (DIS) measurements with polarized beams and/or targets have given us a 
direct measurement of the spins carried by the quarks in the nucleon. A central mys- 
tery has unfolded: in the nucleon, if the u and d valence quarks carry approximately 
equal and opposite spins, where lies the remainder? 

One key contribution that has eluded a definitive explanation is the spin contri- 
bution from strange quarks in the nucleon sea. A large strange quark spin component, 
extracted from recent measurements [S], would be of intense theoretical interest, since 
it would require significant changes to current assumptions. Is this large value of the 
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strange spin due to chiral solitons [H] , a misinterpretation of the large gluon contribu- 
tions coming from the QCD axial anomaly 0IB], or incorrect assumptions of SU(3) 
symmetry In addition, an understanding of the nucleon spin structure is a key 
input to dark matter searches, since the couplings of supersymmetric particles and 
axions to dark matter depend upon the components of the spin. 

It has been known for some time that low energy (and low Q 2 ) neutrino mea- 
surements are the only theoretically robust technique (as robust as, e.g., the Bjorken 
sum rule) for isolating the strange quark contribution. New low energy, intense neu- 
trino beams now make it possible to take greater advantage of this method. The 
FINeSSE experiment, using these newly-available beams along with a novel detection 
technique, will resolve the presently murky experimental picture, providing results 
which are interesting both to particle physics and astrophysics. 

FINeSSE will measure As by examining neutral current neutrino-proton scatter- 
ing; the rate of this process is sensitive to any contributions from strange quarks (both 
s and s) to the nucleon spin. Specifically, As is extracted from the ratio of neutral 
current neutrino-proton (up — > up) scattering to charged current neutrino-neutron 
(un — > fi~p) scattering. The measurement will be made at low momentum transfer 
(Q 2 ~ 0.2 GeV 2 ), in order to unambiguously extract As from the axial form factor, 
Ga- FINeSSE will improve on the latest measurement of neutral current neutrino- 
proton scattering (BNL 734) by measuring this process at a lower-Q 2 , with more 
events, less background, and lower systematic uncertainty. 

In the following sections, we describe some of the previous and current experi- 
ments relevant to the question of strange quarks in the nucleon. We then describe why 
neutral current neutrino-nucleon elastic scattering is sensitive to the strange-quark 
contributions to the nucleon spin. We conclude with a summary of the sensitivity of 
FINeSSE to As (detailed more completely in Chapter |3J). 

2.1.1 Experimental Results: Strange Quarks in the Nucleon 

The role that strange quarks play in determining the properties of the nucleon is 
an area of much experimental and theoretical interest, and is not well understood. 
Deep inelastic scattering of neutrinos on nucleons indicate that strange quarks con- 
stitute a substantial fraction (20%) of the nucleon sea [TQ| . Nevertheless, the latest 
results from parity-violating (PV) electron scattering show that the contribution of 
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strange quarks to the charge and magnetic moment of the nucleon is consistent with 
zero [TTl However, results from DIS of polarized charged leptons indicate non- 
zero contributions of strange quarks to the nucleon spin. In the sections that follow, 
we describe the some of the experimental and theoretical issues which FINeSSE will 
help to address. 

Parity- Violating Electron Scattering 

There is a large and continuing effort to investigate the nucleon structure via Parity- 
Violating (PV) electron scattering. Electron scattering is sensitive to the strange (and 
anti-strange) quark contributions to the nucleon charge and magnetic moment dis- 
tributions. The recently-completed MIT/Bates SAMPLE experiment measured 
the strange-magnetic form factor, G S M , at a momentum transfer Q 2 = 0.1 GeV 2 , 
to be consistent with zero (albeit with large errors). The ongoing Jefferson Lab 
HAPPEX experiment [T2j measures a PV scattering asymmetry sensitive to a com- 
bination of the strange-electric and strange-magnetic form factors (G S E and G S M ) at 
Q 2 = 0.48 GeV 2 . This combination is also consistent with no strange-quark effects in 
the nucleon. HAPPEX continues to search, and will make a measurement at lower 
Q 2 (= 0.1 GeV 2 ) on a helium target in the near future [TH] . 

Upcoming PV electron experiments looking for strange quark effects are the PVA4 
experiment jHj at Mainz, which will measure a combination of G S E and G S M at Q 2 
from 0.1 - 0.48 GeV 2 ; and the GO experiment [T^], which will cover a large Q 2 (0.1 - 
1.0 eV 2 ) range. 

There is a large effort to look for strange quark effects via PV electron scattering. 
Unfortunately, these measurements are not sensitive to the strange axial vector form 
factor (G S A , related to the spin structure and As). This is because the parity violating 
contribution to the axial vector form factor of the nucleon couples via the very small 
vector form factor of the charged lepton (1 — 4 sm6\y)- As a result, this contribution 
to the PV asymmetry is dominated by a large PV radiative correction known as the 
anapole moment jllj . 

PV electron scattering combined with neutrino scattering would be a powerful 
approach to the study of strange quark effects in the nucleon. As will be explained 
below, neutrino scattering is very sensitive to the strange axial-vector form factor 
G S A , and much less so to the strange electric and magnetic form factors G S E and G S M . 
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The opposite is true in PV electron scattering. The measurements from FINeSSE 
could be combined with those from this thorough program of PV electron scattering 
to obtain accurate knowledge of all three strange form factors G S A , G S E , and G S M 

Polarized Lepton Deep Inelastic Scattering 

Results obtained in DIS by polarized leptons from polarized nucleons have been in- 
terpreted (from measurements of the polarized structure function, g\) as evidence 
for a non-zero and negative As. The SMC experiment, e.g., has reported As = 
—0.10 ± 0.05 ^7j. The method of extracting As from g\, however, has been subject 
to much debate, due to mo del- dependent assumptions of SU(3) symmetry EE] 
and to extrapolation of the spin structure function to x — > |7| . In addition, recent 
controversial results from HERMES fH] semi-inclusive DIS indicate a zero or small 
positive value for As (As = +0.03 ±0.03 ±0.01 where the first error is statistical and 
the second is systematic). It should be noted that the HERMES measurement is only 
sensitive to As, as opposed to the SMC measurement of As ± As [T9*] . The errors 
reported in these current "state-of-the-art" DIS measurements of As are typically in 
the range of 0.02 - 0.05. 

As will be shown in the sections below, the interpretation of neutral current 
neutrino-nucleon scattering suffers from none of the theoretical uncertainties inherent 
in the DIS measurements. 

Neutral-Current Neutrino Scattering 

The best measurement to date of neutral current neutrino scattering is Experiment 
E734 at Brookhaven National Laboratory (BNL). They measured neutrino-proton 
(up — > up) and antineutrino-proton (up — > up) elastic scattering (201 using a 170 t 
tracking detector in the BNL wide-band neutrino beam (E v = 1.3 GeV). From a 
sample of 951 up and 776 up elastic scattering events, they extracted differential cross 
sections (da/dQ 2 ) for 0.4 < Q 2 < 1.1 (GeV/c) 2 (Fig. EE}. These data were fit to a 
description of up — > up and up — > up to obtain the results shown in Figure 12.21 The 
results from this fit were often cited to support the claims from the DIS experiments 
at the time, that As was non-zero and negative. 

This experiment also simultaneously measured the neutrino and antineutrino neu- 
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Figure 2.1: BNL E734 results on up — > z/p and z/p — > up 120/ . The solid lines are fits 
to the data. 

tral current to charged current ratios: 



a (up — > up) 

R v = ZL. = 0.153 ± 0.007 ± 0.017 

a(un — > p, p) 



(2.1; 



and 



o~(up — > up) 
o~(up — > fi + n) 



0.218 ±0.012 ±0.023 



(2.2) 



calculated over the interval 0.5 < Q 2 < 1.0(GeV/c) 2 . The dominant error in these 
ratios was an 11% systematic. 
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Figure 2.2: Results from a fit of the BNL E734 vp — > vp and up — > up scattering 
data ]2(Jj indicating the preferred values of rj (= —As) and Ma- In this fit, Ma was 
constrained to the world-average value at the time Ma = 1.032 ± 0.036 GeV. 

Reanalysis of BNL E734 

The non-zero value for As as obtained by BNL E734 was later reexamined j2I] . This 
reanalysis more carefully considered the effects of strange contributions to the vector 
form factors and the Q 2 evolution of the axial form factor in the differential cross 
sections for up and up elastic scattering. A value of As = —0.21 ±0.10 was extracted. 
An additional form factor uncertainty of ±0.10, determined from fits to the data with 
different assumptions about the vector form factors and evolution of the axial form 
factor (via the axial vector mass, Ma), should be assigned. 

Another group has reexamined the BNL E734 result on the neutral current to 
charged current ratios, R v and Rp [22] • These ratios are sensitive to the axial form 
factor, and avoid the systematic uncertainty of the neutrino flux and nuclear model 
corrections. While these ratios hold the promise of a superior method to extract As, 
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the experimental errors from BNL E734 were too large to provide a definitive answer, 
and the conclusions of this analysis were consistent with the previous reanalysis of 
the data [2T]- 

Another reanalysis of these data has been recently performed ^H]- This approach 
is interesting in that it uses the latest data from HAPPEX on the electric and magnetic 
strange form factors in an effort to reduce the systematic errore from these form factors 
on this meausurement. The results show the viability of this method of combining 
the PV electron and neutrino scattering into one analysis, but are limited by the large 
systematic errors on the BNL 734 data sample. 

Note that the BNL E734 data have been reanalyzed at least three times since 
being published in 1987; this points to growing appreciation of the fact that neutral 
current neutrino scattering is an excellent probe of As. Unfortunately, though, the 
precision of the BNL E734 data limits the conclusions that can be drawn from it. 

How FINeSSE Will Improve These Measurements 

The neutrino neutral current to charged current ratio as measured by BNL E734 
(Eq. I2.1J1 is reported with a statistical error of 5% and a systematic error of 11%. 
This ratio was made by integrating over the Q 2 interval 0.5 < Q 2 < 1.0 GeV 2 . 

FINeSSE will improve upon this measurement in several ways: 

• The neutrino neutral current to charged current ratio measurement will be 
made as a function of Q 2 in the interval 0.2 < Q 2 < 1.0 GeV 2 . FINeSSE's 
measurement will be made at a lower-Q 2 value, where form factors are more 
easily interpreted. In addition, the shape of the ratio as a function of Q 2 holds 
additional information about the evolution of the form factors. 

• With the proposed FINeSSE detector and run plan, the statistical errors in the 
0.2 < Q 2 < 0.4 GeV 2 -bin alone will amount to a relative error on the ratio of 
only 2%. 

• Systematic errors have been estimated with detailed simulations of the detector 
to be 5%. Much of this reduction in systematic error is due to the greatly 
reduced background in FINeSSE. The background to the up — > up reaction is 
estimated at 26% (and is likely to be further reduced) compared to 40% in BNL 
E734 [201 • 
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In summary, FINeSSE will be able to make a measurement of the neutral-current 
to charged-current ratio with a 6% total error down to 0.2 GeV 2 . BNL E734 made 
a measurement with 12% total error down to only 0.5 GeV 2 . This will allow a sig- 
nificant improvement to the uncertainty on the extracted value of As as described in 
Section ^. 1.31 and Chapter |3J 

2.1.2 Relevance to Searches for Dark Matter 

Understanding of the spin contribution to the nucleon of the strange quarks is impor- 
tant for certain searches of dark matter [2S] . In i?-parity-conserving supersymmetric 
models, the lightest supersymmetric particle (LSP) is stable and therefore a dark- 
matter candidate; in certain scenarios, the relic LSP density is large enough to be 
of cosmological interest. Experimental searches for cosmic LSPs can be competitive 
with accelerator-based searches 

In the case where the LSP is the neutralino, cosmic LSP can be detected either 
directly, through elastic neutralino scattering in an appropriate target/detector, or 
indirectly. The indirect method involves detection of high-energy neutrinos from 
the center of the sun, where the heavy neutralinos accumulate and subsequently 
annihilate. If the neutralino mass is larger than the W mass, annihilation into gauge 
bosons dominates and this gives rise to high-energy neutrinos that can be detected 
on earth. The expected rate for this process also depends on the elastic neutralino- 
nucleon scattering cross section. 

The neutralino-nucleus elastic-scattering cross section contains a spin-dependent 
and a spin-independent part. The spin-dependent part is given by 



where Gp is the Fermi constant, m r the reduced neutralino mass, J the nucleus spin, 
and 



where the sum is over quark flavors and the coefficients on are functions of the compo- 
sition of the neutralino in terms of the supersymmetric partners of the gauge bosons. 
The factors Aqf and Ag™ are the quark contributions to the proton or neutron spin. 



a = —G 2 F m 2 r A 2 J(J+l), 



A=j(a p (S p ) + a n (S n )) ; 
here (Sp( n )) is the average proton (neutron) spin in the nucleus and 
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It is established |171 E] that Am and Ad have opposite signs. From the above, 
it should be clear that knowledge of As not only is important for the interpretation 
of any limits from such dark matter searches, but it could also influence the choice 
of detector material for direct searches [53], making nuclei with either proton- or 
neutron-spin excess optimal, depending on its value and sign. 



2.1.3 A Measurement of As via Neutral- Current Neutrino 
Scattering 

In neutral current elastic (NC) neutrino- nucleon scattering (vN — > vN), any isoscalar 
contribution (such as strange quarks) to the nucleon spin will contribute to the cross 
section. This is in contrast to the charged current quasi-elastic (CCQE) scattering 
process [yn — > fi~p) where only isovector contributions are possible. 

FINeSSE will use this feature of neutrino scattering to measure any contribution 
of strange quarks to the spin of the nucleon. The fact that NC neutrino scattering 
is sensitive to strange quark (isoscalar) spin in the proton, and CCQE neutrino scat- 
tering is not, will be exploited by measuring a ratio of these two processes; this will 
eliminate a number of experimental and theoretical errors. 



The Neutral Weak Axial Current of the Nucleon 



The axial part of the weak neutral current may be written |2*T) l2*o] , 



N\Al\N) 



G p 

G p 

71. 



N 



N 



G A (Q* 



N 



N 



(2.3) 

(2.4) 
The 



where Ga is the axial form factor and r z = ±1 for protons (+) or neutrons 
strange axial-vector form factor, G S A , is identified with the 37^755 term which is As, 
the spin carried by the strange quarks. So, the non-strange (u and d) quark axial 
current is accounted for in Ga, known at Q 2 = from neutron beta decay to be 
Ga(Q 2 = 0) = 1.2673 ± 0.0035 [27]. The (unknown) strange quark axial current is 
subsumed in G S A . A similar decomposition may be obtained for the vector part of the 
neutral weak current [5TJ [2E] in terms of the two non-strange vector form factors, F\ 
and F2, and the corresponding strange quarks parts, F( and F|. 
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As shown above, the strange axial form factor, G S A , in the limit of zero momentum 
transfer (Q 2 = 0), is identified with the strange quark contribution to the nucleon 
spin, As. It is worth mentioning that the arguments leading to the connection be- 
tween G S A at low momentum transfer and the strange spin content, As, that can 
be extracted from DIS data in the scaling limit, are essentially the same as those 
that were used to derive the Bjorken sum rule, which connects spin-dependent DIS 
structure functions to the coupling constants in neutron decay and is considered one 
of the most fundamental predictions of QCD. They are both based on the operator 
product expansion, which expresses the moments of structure functions in terms of 
matrix elements of local operators and perturbatively calculable Wilson coefficients. 



Neutrino Cross Sections and Form Factors 



The differential cross section, da/dQ 2 , for NC and CCQE scattering of neutrinos and 
antineutrinos from nucleons can be written as a function of the nucleon form factors 
F\ and F 2 (both vector) and G\ (axial vector) 



da 
dQ 2 



CW Z 



(2.5) 



with kinematic factor, W = 4E u /m p — Q 2 /m 2 . The + (-) sign is for neutrino (an- 
tineutrino) scattering. Q 2 is the squared-four-momentum transfer. The A, B, and C, 
contain the form factors: 



.4 
B 

C 



\G\{\ + r) - (F 2 - rF|)(l - r) + AtF.F, 
~G 1 {F 1 + F 2 ), 



1 m 

_L_ZfG 2 + F 2 + tF 
16Q 2lGl+ 1 + : 



(2.6) 
(2.7) 

(2.8) 



Up to this point, this formalism is valid for both NC and CCQE neutrino- (and 
antineutrino-) nucleon scattering. The difference between the NC and CCQE pro- 
cesses is accounted for with different form factors (Fx, F 2 , and G\) in each case. 

For NC neutrino scattering the axial-vector from factor (as described in Sec- 
tion |2^S1) may be written in terms of the known axial form factor plus an unknown 
strange form factor: 



G, 



G* 



(2.9) 
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The vector form factors may also be written in terms of known form factors plus 
a strange quark contribution, 



•1,2 



sin 2 # 



K 2 - K2] r z - sin 2 6 W + Ffr] - -F% 



(2.10) 



where F^ n ^ is the Dirac form factor of the proton (neutron) and F 2 ^ is the Pauli 
form factor. The CVC hypothesis allows us to write the same form factors in these 
equations as those measured in electron scattering. Therefore, the only unknown 
quantities in these equations are the strange vector form factors F( 2 . 

The differential cross section for neutrino scattering at low Q 2 (Eq. 12.5(1 is dom- 
inated by the axial form factor, G\. This can be seen be examining Equations 12.51 
and 12.61 In fact, at low-Q 2 , da/dQ 2 oc G\. This is a crucial point, and it is what 
makes NC neutrino scattering the best place to look for the effects of strange-quarks 
in the nucleon spin. It also makes the results less sensitive to the strange vector form 
factors, F( 2 . 



Q 2 Dependence of the Form Factors 

All of the form factors are, most generally, functions of Q 2 . The values of the non- 
strange form factors at Q 2 = are known from the static properties of the nucleon 
(e.g. charge, magnetic moment, and neutron decay constant). How the form factors 
change with increasing Q 2 is less well known and must be addressed. The form factors 
are commonly parameterized with a dipole form, where the dipole mass set by various 
measurements. The evolution of F^ 2 is known via numerous experiments on electron 
scattering; the vector dipole mass, My is 0.843 GeV/c 2 . The Q 2 dependence of the 
axial form factor, Ga, is measured via CCQE neutrino scattering; the world average 
data yield an axial mass Ma = 1.026 ± 0.021 [2*7] . 

Both the strange form factors, Ff 2 and G S A and their Q 2 evolution are unknown. 
It is most common to assume the same Q 2 dependence as for the non-strange form 
factors, using My for F( 2 and Ma for G S A . The uncertainty on the Q 2 evolution intro- 
duces some uncertainty in the extraction of G S A from a measurement. A measurement 
at Q 2 ~ 0.2 GeV 2 , however, would keep this contribution to the uncertainty at a low 
level, as the value of the form factor differs by only 20% from Q 2 = to Q 2 = 0.2 
GeV 2 . 
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Nuclear Physics Corrections 

The expression for the cross section (Eq. 12. 5|) is for scattering from free nucleons. 
Since FINeSSE will have a target that consists, in large part, of nucleons bound in 
carbon, consideration will need to be given of the effects of this binding [2U ED] . The 
corrections can be rather large when considering the absolute event rate, and can 
depend greatly on the model employed, because the amount and quality of available 
neutrino data to constrain such models, is lacking. In principle, the nuclear models 
can be constrained with the high-quality electron data available; this, however, is a 
work in progress. 

These effects become less of a concern when ratios of cross sections are consid- 
ered [HU 1221 HUE EH]- The initial and final states of the hadrons involved are quite 
similar in both NC and CCQE neutrino scattering; as a result, the corrections em- 
ployed for either channel should be similar as well. FINeSSE will utilize this fact with 
a measurement of As as explained in the following sections. 

The Neutral Current to Charged Current Ratio 

The NC neutrino scattering cross section depends strongly on G\ and therefore on G S A , 
the quantity of interest. An absolute measurement of the cross section, however, is an 
experimental challenge; the level of precision achievable for this measurement would 
not yield the desired precision for G S A . An absolute prediction is also a challenge 
from a phenomeno logical standpoint, as uncertainties in form factors and nuclear 
corrections can be large. 

It is possible to extract G S A to the desired precision by measuring the ratio of 
NC to CCQE event rates. A measurement of the ratio of NC to CCQE neutrino 
scattering event rates may be measured with greater precision, since many systematic 
uncertainties cancel in the ratio such as the neutrino flux and correlated reconstruction 
efficiencies. Theoretical uncertainties are also reduced in this quantity as many of 
these uncertainties are correlated between NC and CCQE scattering. FINeSSE will 
use this method to measure G S A . 

First, consider the ratio of NC neutrino-proton to NC neutrino-neutron scattering, 




o~{yp — * up) 
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This ratio is very sensitive to G S A The NC neutrino nucleon scattering cross 

section is proportional to G\ as explained above. A non-zero value for G S A will pull 
the denominator of this ratio one way, and the numerator the other, due to the r 2 
factor in Equation | 



R(p/n), however, is likely to be very difficult to measure accurately in a neutrino 
scattering experiment, because of the intrinsic difficulties and uncertainties involved 
with neutron detection. For this reason FINeSSE will focus on a measurement of the 
ratio of the NC neutrino-proton scattering (up — > up) to CCQE neutrino scattering 
(un — > fi~p). This ratio, 

R(NC/CC) = *(^"P) (2 . 12) 
a(vn — > n~p) 

can be more accurately measured and is still quite sensitive to G S A (albeit less so than 
R(p/n)). The numerator depends upon G S A as explained in the formalism introduced 
above. The denominator does not as the un — > fi~p process is sensitive only to 
isovector quark currents and not to isoscalar currents (such as that from strange 
quarks) . 

The up — * up and un — > fi~p differential cross sections (weighted by the calculated 
FINeSSE flux) as calculated with Equation 12.51 are plotted in Figure 12.31 for G S A = 
—0.1, 0.0, +0.1; this shows that the cross section for up — > up depends strongly upon 
G A . The cross section for un — > fi~p is independent of G S A , so the ratio of flux- 
weighted cross sections (and therefore the event rates) of these two processes depends 
upon G S A . 

This dependence is also shown in Figure l2~H as a function of G S A for three different 
Q 2 bins. In the Q 2 = 0.25 GeV 2 bin, the sensitivity of the NC/CC ratio on G A is 
approximately 1.2. The relative uncertainty in the NC/CC ratio, AR/R is related to 
the absolute uncertainty on G S A , cr{G A ) by 

v{Ga) = Y2~R' ^ ' 

Considering this and recalling that G S A (Q 2 = 0) = As yields the conclusion that 
a 5% relative measurement of the NC/CC ratio at Q 2 ~ 0.25 GeV 2 would enable an 
extraction of As with an error of 0.04. This error is comparable to that quoted in 
the latest extractions of As from charged lepton DIS [TH] , 
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Figure 2.3: Flux-weighted differential cross section as a function of Q 2 for up — > up 
(a) and un — > fi~p (b) scattering together with the cross section ratio of the these two 
processes (c). These quantities are shown for G S A = 0. (solid), = —0.1 (dashed), and 
= +0.1 (dotted). The un — > fi~p process does not depend upon G S A . 
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Figure 2.4: Ratio of flux-weighted cross sections of the up — > up and un — > \i~p 
processes as a function of the axial form factor G S A atQ 2 = 0.25 (solid), 0.45 (dashed), 
and 0.65 (dotted) GeV 2 . 

2.1.4 FINeSSE Sensitivity to As 

NC, CCQE, and background events in the proposed detector have been simulated in 
a detailed manner; a reconstruction procedure has been performed, to calculate an 
estimated sensitivity. All conceivable effects were considered; systematic errors were 
estimated. The experimental errors considered were those due to statistics and to 
various systematics. The statistical errors were calculated based on a 9 ton (fidu- 
cial) detector running for two years (6 x 10 20 protons-on-target). The systematic 
uncertainties considered include those due to free-to-bound up — > up scattering rate, 
reconstruction efficiencies, background estimation, and Q 2 reconstruction. In addi- 
tion, theoretical uncertainties from nuclear model dependence and from form factor 
estimation were considered. The details of this procedure are described in Chapter El 

A fit to the simulated data in the Q 2 region where the detector has reasonable 
acceptance yields an experimental uncertainty in As of ±0.039; the combined uncer- 
tainty from the axial mass, M A , and F% form factor is ±0.025. 

Based on these results, it has been determined that FINeSSE, with a design and 
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plan described in the Chapters below, can make an accurate measurement of As at 
down to Q 2 ~ 0.2 GeV 2 . This will enable FINeSSE to answer an important and 
unresolved question about the structure of the proton. 

2.2 Disappearance 

The FINeSSE detector can be used in conjunction with the MiniBooNE detector 
to substantially improve our understanding of neutrino oscillations at high Am 2 by 
looking for a neutrino energy-dependent deficit of event rates compared to a no- 
oscillation hypothesis. This search is motivated by astrophysical models for the pro- 
duction of heavy elements in super novae. If MiniBooNE observes a signal, a combined 
FINeSSE and MiniBooNE ("FINeSSE+MiniBooNE") run represents the next step in 
determining the underlying physics model of the oscillation. These two motivations 
are not directly connected: if MiniBooNE does not see a signal, the astrophysical case 
still makes this study compelling. 

In this section, we first provide a brief overview of the formalism for neutrino oscil- 
lations. Second, we introduce the LSND signal, along with theoretical interpretations 
involving sterile neutrinos. Third, we describe the astrophysical motivations for the 
disappearance search. Lastly, we describe the capability of a FINeSSE+MiniBooNE 
joint analysis of disappearance. 

2.2.1 Neutrino Oscillation Formalism 

"Neutrino oscillations" occur when a pure flavor (weak) eigenstate born through a 
weak decay changes into another flavor as the state propagates in space. This can 
occur if two conditions are met. First, the weak eigenstates can be written as mixtures 
of the mass eigenstates, for example: 

v e = cos 9 V\ + sin 9 v-i 
v ii = — sin 9 V\ + cos 9 v 2 

where 9 is the "mixing angle." The second condition is that each of the mass eigen- 
state components propagate with a different frequency, which can occur only if the 
masses are different. We define the squared mass difference as Am 2 = \m 2 — m 2 \ > 0. 
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In a two-component model, the oscillation probability for — > v e oscillations is then 
given by: 



Prob {v.. 



sin 29 sm z 



1.27 Am 2 (eV 2 ) L (km) 



(2.14) 



E (GeV) 

where L is the distance from the source, and E is the neutrino energy. 

Most neutrino oscillation analyses consider only two-generation mixing scenarios, 
but the more general case includes oscillations between all neutrino species. For the 
case of the three Standard Model species, this can be expressed as: 
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The oscillation probability is then: 



- 4 % (KiUfnUajUh) sin 2 [Am|(L/4E)] 

i>j 

+ 2^3 (U* ai U pi U aj U* pj ) sin[Am 2 .(L/2E)] , (2.15) 

i>j 

where Am 2 ^ = |m 2 — m 2 - 1 . Note that there are three different Am 2 (although only two 
are independent), and three different mixing angles. This method can be generalized 
to include more neutrino species in Beyond-the-Standard Model Theories. 

Although in general there will be mixing among all flavors of neutrinos, two- 
generation mixing is often assumed for simplicity. If the mass scales are quite different 
(e.g., m 3 » m2 » mi), then the oscillation phenomena tend to decouple and 
the two-generation mixing model is a good approximation in limited regions. In 
this case, each transition can be described by a two-generation mixing equation. 
However, it is possible that experimental results interpreted within the two-generation 
mixing formalism may indicate very different Am 2 scales with quite different apparent 
strengths for the same oscillation. This is because, as is evident from equation 12.151 
multiple terms involving different mixing strengths and Am 2 values contribute to the 
transition probability for u a — > Pp. 
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From equation 12 .141 one can see that the oscillation wavelength will depend upon 
L, E, and Am 2 . For short baseline experiments, sensitivity to oscillations is in a 
range of > O.leV 2 , which will term the "high Am 2 region. The oscillation amplitude 
will depend upon sin 2 26. 

2.2.2 Experimental Results: The LSND Signal 

One of the most exciting questions in high energy physics, at present, is whether 
the "LSND signal" is due to neutrino oscillations. If the signal is confirmed in 
the MiniBooNE experiment, then this is an indication for new physics beyond the 
Standard Model. Fermilab will want to be poised to pursue this result. A FI- 
NeSSE+MiniBooNE run will be the first step. 

The LSND experiment has observed a 4<r excess which can be interpreted as 
oscillations between muon and electron neutrinos [HH]- The beam was produced at 
LANSCE at LANL, with 800 MeV protons interacting with a water target, a close- 
packed high-Z target, and a water-cooled copper beam dump. The highest statistics 
came from neutrinos produced by decay at rest (DAR) of muons, with 20 < E < 50 
MeV. However, the lower statistics decay in flight (DIF) v^s were also analyzed. The 
liquid scintillator detector was located 30 m from the beam dump. Hence the L/E of 
the experiment was ~ 1 m/MeV. As a result, if the excess is interpreted as oscillations, 
the allowed region is located at high Am 2 . 

LSND is the only short-baseline experiment to have observed evidence for oscil- 
lations. Other short-baseline experiments have searched and seen no signal. Those 
most relevant to this proposal are Karmen jSH] and Bugey [21] • Karmen, which ran 
at the ISIS facility at Rutherford Labs, was similar in concept to LSND, using a DAR 
beam; but the detector was smaller, and the beam of lower intensity. Most impor- 
tantly, it had an L of 17 m; in this way, it can be thought of as a "near detector" 
for LSND. Playing the null signal in Karmen against the observed excess in LSND 
results in the allowed regions shown in Figure l2~B1 |38j . The Bugey reactor experiment 
rules out sin 2 26 > 0.04 at 90% CL. in a search for v e disappearance using a reactor. 
Assuming that oscillations respect time reversal, and that only the three standard 
model active species are involved, Bugey's result can be taken as an excluded region 
for LSND. 



The MiniBooNE experiment is a designed to decisively address the LSND signal. 
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Figure 2.5: LSND and Karmen joint analysis allowed region. 



This experiment will complete its Phase I (neutrino) run in mid-2005 and is expected 
to request further (Phase II) running for the period thereafter. By the time FINeSSE 
is on-line, the LSND signal will be tested in both v and V modes. 

If MiniBooNE confirms LSND, and all other oscillation results remain as they 
presently stand, then this necessarily implies new physics. The oscillation signals 
from solar neutrinos, atmospheric neutrinos, and LSND cannot be simultaneously 
fit with the three Standard Model neutrinos. A favored method for expanding the 
theory to allow LSND is to invoke sterile neutrinos (u s ). These are neutrinos which 
do not interact via the W or Z, but can couple to the Standard Model "active" 
neutrinos through oscillations. The most minimal extension is to introduce a single 
light sterile neutrino. This extra neutrino is not ruled out by cosmology Light 
sterile neutrinos can appear in supersymmetry, extra dimensions and GUTs |40j . 
These can all accommodate more than one light sterile neutrino, but we will confine 
our discussion to one for simplicity. 

Figure l2.6Fl eft) shows a cartoon of how the squared masses and mixings might 
be arranged if a single sterile neutrino is introduced into the theory. The vertical 
axis is logarithmic and arbitrary. The bars indicate the flavor content of each mass 
state. The LSND signal is explained by the largest squared mass splitting, with the 
transition — > v e . Because there is a triplet of neutrinos with nearly the same mass, 
and one large splitting, this is called a "3+1" model. Note that the transitions — > v s 




Figure 2.6: Left: an example of a 3+1 mass spectrum. Right: allowed region for 
v s in 3+1 models. The exclusion regions for CDHS and CCFR are also shown. 
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and v e — > v s must also be allowed for the same Am 2 . The allowed regions for 3+1 
for fits which include LSND, Karmen, Bugey, and two disappearance experiments 
(CDHS and CCFR84) are shown in Figure right) |JT]. These allowed regions can 
be addressed by FINeSSE+MiniBooNE as discussed below. 

2.2.3 Astrophysical Motivation for High Am 2 Disappearance 
Searches 

The existence of at least one sterile neutrino in the high ~ 1 eV mass range has 
interesting consequences for the heavy element abundance in the universe. In fact, 
oscillations were predicted on the basis of this abundance before the LSND signal 
was presented j2j- The allowed range extends beyond the region constrained by the 
LSND signal. Hence, whether or not MiniBooNE confirms LSND, searches for active- 
to-sterile oscillations at high Am 2 remain motivated by this astrophysical question. 

Active-to-sterile neutrino oscillations in the late time post-core-bounce period of 
a supernova will affect the r-process, or rapid neutron capture process. This presents 
a mechanism for producing substantially more heavy elements (A > 100), solving 
the long-standing problem of the high abundance of uranium in the universe. The 
FINeSSE+MiniBooNE search addresses allowed parameters for this solution. 

A favored mechanism for producing heavy elements is through the r-process in 
the neutrino-heated ejecta of a Type II or Type I/c supernova. In this model, during 
the period of the neutrino-driven wind, which lasts for ~ 10s, "seed" elements with 
A between 50 and 100 capture neutrons to produce the elements with A > 100. The 
problem is that in most models the neutron-to-seed ratio, R, is too low for production 
of the heaviest elements [12] • In fact, detailed simulation show that a phenomenon 
called the a-effect, in which neutrons are frozen into alpha particles that do not 
recombine to form heavier elements in the requisite time period, renders the neutron- 
to-seed ratio downright "anemic" |43j . 

Various solutions have been proposed. One option is to resort to the competing 
theory of neutron star mergers. The problem with this scenario is that mergers are too 
rare to produce the observed abundance of heavy elements jH] . Another alternative 
is to introduce physics which adjusts the neutron-to-seed ratio. This can be done 
by modifying the expansion rate, the entropy per baryon or the n/p ratio - all of 
which will affect the neutron-to-seed ratio. Adjusting any of these three requires 
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Figure 2.7: Allowed oscillation parameters for the oscillation-enhanced r-process. 

invoking new physics in the model. We explore the last alternative here: introducing 
v e ^sterile oscillations, which, when combined with matter effects, enhance the 
production of neutrons over protons. 

The idea HHj is that production of neutron-rich elements requires a neutron- 
rich environment. To the level that the processes v e + n — > p + e~ and v e + p — > n + e + 
are in balance during the neutrino-driven wind, there is no net excess of neutrons. In 
its simplest form, neutrino oscillations between an electron and sterile flavor would 
not upset the balance because v e and v e will oscillate at the same rate. However, when 
one introduces matter (or MSW jlHl) effects, neutrino and anti-neutrino oscillations 
are modified with opposite sign in an electron-rich environment. Oscillations of v e to 
sterile neutrinos are enhanced, while v e are depressed. This can produce a substantial 
neutron excess by removing the offending u^s. The a-process removes some neutrons, 
but stalls once the protons are devoured, leaving sufficient neutrons to produce the 
high- A elements. 

In this model, the neutron-to-proton ratio (usually characterized as the function 
Y e = 1/(1 + n/p))), depends upon the choice of Am 2 and sin 2 29. The condition for 
a successful r-process is Y e < 1/2. The smaller the value of Y e , the larger the high- A 
abundance. Figure l2~7l shows Y e as a function of the v e — > v s oscillation parameters. 
This shows that there are a wide range of "robust" solutions |S] • 

One can connect the allowed space for — > v s and — > v e to the allowed 
region for v e — > v s within 3+1 and 3+2 models. If MiniBooNE sees a signal, this will 
be a great victory for the oscillation-enhanced r-process model. In the case where 
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MiniBooNE does not see a signal, there remains a large parameter space open to this 
model. At present, disappearance in FINeSSE+MiniBooNE represents the only 
way to access that parameter space. 

2.2.4 FINeSSE+MiniBooNE Capability for Disappearance 

The FINeSSE detector can be combined with MiniBooNE to explore allowed regions 
for oscillations to sterile neutrinos via disappearance. In this analysis, FINeSSE 
serves as a near detector to accurately measure the flux, and MiniBooNE serves as 
the far detector where a deficit may be observed. This is a unique capability - there 
are no other short baseline disappearance experiments in the world. If Mini- 
BooNE observes a signal, FINeSSE+MiniBooNE will be a crucial next-step toward 
understanding the result. If MiniBooNE does not observe a signal, this region is still 
interesting because of its relevance to astrophysics. 

Figure 12.81 shows the FINeSSE+MiniBooNE expectation for the default design, 
with 6 x 10 20 protons on target, in neutrino mode. Also shown are the 3+1 allowed 
regions for fits to LSND, atmospheric, and solar (as described above); and the expec- 
tation for MiniBooNE prior to FINeSSE running. MiniBooNE will be able to address 
the lower 3+1 allowed regions. The largest 3+1 allowed region, however, can only be 
addressed by the FINeSSE+MiniBooNE combination. This combination of FINeSSE 
and MiniBooNE is therefore very powerful; it alone is able to address the full 3+1 
picture. 

The "standard" configuration for FINeSSE and MiniBooNE simultaneous running 
places the FINeSSE detector at 100 m from the primary beryllium target, with the 
25 m absorber installed in the beamline. The angular acceptance from the target 
to FINeSSE is 25 mrad, and to MiniBooNE is 10 mrad. In this analysis, we accept 
only neutrinos which traverse both detectors, meaning that we use only the inner 1 
m radius (10 mrad acceptance from target) of FINeSSE. Event rates are for a 9 ton 
fiducial volume and 6 x 10 20 protons in neutrino mode. 

In Chapter 5, we provide details on how the sensitivity shown in Figure |2~H1 was 
obtained. We explain why the standard configuration is best for the analysis. We 
also describe the method for determining the sensitivity, which compares the energy 
distribution of events in the near and far detector. This method accounts for both 
statistical and systematic errors. 
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Figure 2.8: The parameter space covered by FINeSSE+MiniBooNE for disappear- 
ance (labeled "FINESE"). Also shown: allowed regions in 3+1 models given the LSND 
signal (solid), existing exclusion regions from CDHS and CCFR, and the expected ex- 
clusion region for MiniBooNE. 
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2.3 FINeSSE on the Booster Neutrino Beamline 

The Booster neutrino beamline is the only existing beamline at Fermilab or around 
the world where this physics can be accomplished. The FINeSSE As measurement 
requires a clean, low energy neutrino spectrum, as is produced by the Booster neutrino 
production target. The oscillation physics goals of FINeSSE require the energies 
and baselines available to experiments on the Booster neutrino beamline. These 
requirements make it impossible to perform this measurement at Fermilab's other 
existing neutrino beamline, NuMI. 
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Figure 2.9: The flux spectrum of the Booster neutrino beam and the NuMI neutrino 
beam at their near detector locations. 

Figure |2~9"1 shows the NuMI beam flux at the MINOS near detector location, 290 m 
from the NuMI decay pipe. The Booster neutrino spectrum at 75 m from the end of 
the Booster neutrino beamline decay pipe with the 25 m absorber in position (100 m 
from the production target), is shown in Figure l3~2l As indicated, the average energy 
of neutrinos from the Booster is 700 MeV, with virtually no neutrinos beyond 3 GeV. 
This neutrino energy distribution is excellent for making the As measurement. It is 
an energy large enough to be beyond the region where low-energy nuclear corrections 
are significant, yet not so large where pion production and DIS scattering backgrounds 
are high. The NuMI flux, however, has an average energy above 7 GeV and a tail that 
extends past 30 GeV. This energy distribution creates large pion and DIS scattering 
rates that increase the background to NC neutrino elastic scattering considerably. 

In addition, this flux of neutrinos around the MINOS near detector enclosure will 
create a large flux of low-energy neutrons from neutrino interactions in the earth. This 
background is much larger in the MINOS area when compared to that in FINeSSE. 
The results of simulation of this effect are shown in Figure 12.101 Note that in the 
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Figure 2.10: Ratio of neutron background, from neutrino interactions in the sur- 
rounding earth, to NC-like detector events using the FINeSSE detector in the Booster 
neutrino beam (top) and the NuMI neutrino beam (bottom), shown as a function 
of energy energy deposited in the detector. This plot assumes 100% detection and 
reconstruction efficiency. 

lowest energy bin, the background is 14 times higher at the MINOS near detector 
location. For these reasons, the As measurement can not be made in the NuMI 
beam. 

In order to explore the oscillation regions discussed in Section l2.2.41 a two detector 
comparison is required, and L/E for the far detector must be on the order of 1 m/MeV. 
This is not possible in the NuMI beamline for two reasons. First, there is no tunnel 
for a detector far enough upstream of the MINOS near detector (equivalent to the 
FINeSSE detector) to measure the beam before the neutrinos have oscillated. Second, 
the high average neutrino energy of the NuMI beam means that the baseline for the 
far detector for an oscillation experiment would have to be on the order of 800 m, a 
much longer distance than is available in the NuMI beamline. These two consideration 
lead to the conclusion that FINeSSE must run on the Booster neutrino beamline. 
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Chapter 3 

The Neutrino Beam and Expected 
Event Rates 

3.1 The Booster Neutrino Beam 

The Booster neutrino beamline presently delivers beam to the MiniBooNE exper- 
iment; as the FINeSSE detector will be placed upstream of MiniBooNE, the one 
beamline will provide neutrinos to both experiments. The projections for protons on 
target (POT) in this section are based on a conservative interpretation of the Proton 
committee report [3]. 

3.1.1 Beam Intensity Requirements 

During the fall 2003 shutdown, several improvements in the Booster were made. These 
upgrades are expected to provide routine peak operation with 5 x 10 12 protons/batch 
and 5 Hz for the Booster neutrino beamline. The efficacy of these improvements will 
be understood prior to FINeSSE running, and there should also be sufficient time to 
implement additional improvements if the goals are not met by the end of 2004. 

By the summer of 2003, the Booster was routinely delivering more than 5xl0 12 
protons/batch for Stacking for Run II, demonstrating that Booster can achieve the 
batch intensity required for FINeSSE. The issues are reducing and controlling losses 
at this intensity, and achieving the required repetition rate. The principal improve- 
ments during the Fall of 2003 were modifications to the doglegs to reduce losses, 
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installation of two large aperture RF cavities to reduce losses at these two locations, 
the installation of collimators to control losses, and modifications to the RF and mag- 
net subsystems to allow an increase of the equipment repetition rate to 7.5 Hz. Once 
these improvements are operational, it is expected that the above ground radiation 
will be the limit on Booster operation; however, this limit is well above what is needed 
during the FINeSSE era. In addition, in 2004, Fermilab and Columbia University are 
expected to develop a robot for measuring the losses in the Booster during beam 
operation, which will help to understand these losses in detail. 

Although the Booster equipment may be able to achieve 7.5 Hz, the MiniBooNE 
horn imposes a limit of 5 Hz. If the Booster were to achieve 5xl0 12 protons/batch 
at 5 Hz for an hour, the MiniBooNE target would receive 9xl0 16 protons per hour. 
This is considered a nominal performance level, however, and it is not expected to 
persist for a week, (much less for an entire year). 

To relate a nominal performance to the number of protons delivered per year, one 
can define an annual efficiency. The analysis used here follows the same steps given in 
the Proton committee report 0j. The annual efficiency must include factors to account 
for the number of weeks actually scheduled for beam operation in a year; the reliability 
of the Proton Source (Linac, Booster, and beam transfer lines) during those scheduled 
weeks; and the operational efficiency for actually achieving 5x 10 12 protons/batch and 
5 Hz. The number of weeks scheduled per year is determined by the Director's Office 
and is taken to be in the range 42 to 44 weeks. The reliability of the Proton Source 
has been measured by MiniBooNE to be in the range 0.90 to 0.94. The operational 
efficiency is estimated to be 0.90 [17] . Combining these factors one obtains an annual 
efficiency of 0.66 to 0.72. 

By the time FINeSSE would start to run, however, NuMI will also be taking 
beam. NuMI is expected to use five Booster batches per Main Injector cycle. NuMI 
is expected to share the same Main Injector cycle as Stacking for Run II, and Stacking 
is expected to take two Booster batches per Main Injector cycle. The Main Injector 
cycle time is expected to be about two seconds. With these assumptions, NuMI plus 
Stacking will require seven batches every two seconds, which is an average rate of 
3.5 Hz. At the moment, some of the Booster equipment requires two "prepulses" 
with no beam, or 1 Hz. Thus, the bandwidth required by NuMI, Stacking, and the 
prepulses is 4.5 Hz. This leaves 3 Hz for delivering beam to the MiniBooNE target, 
assuming a total Booster bandwidth of 7.5 Hz. This is 60% of the maximum 5 Hz 



35 



which the Booster neutrino beamline should be receiving in 2004. If prepulses can 
be eliminated, then this can be used to add 1 Hz to this beamline, but this proposal 
does not count on that. 

Thus, one expects a nominal performance of the Booster neutrino beam for FI- 
NeSSE of 5xl0 12 protons/batch and 3 Hz. Given the range for annual efficiency 
calculated above, one calculates the POT/year for FINeSSE as 5xl0 12 x 3 Hz x 
3.15xl0 7 sec/yr x (0.66 to 0.72) = (3.12 to 3.40) xlO 20 POT/yr. 

This proposal assumes delivery of 3.0xl0 20 POT/yr, conservatively below the 
range quoted above. 

3.1.2 Booster Neutrino Beam Production 
The Neutrino Flux 

The neutrino beam is produced by the 8 GeV Fermilab Booster which currently 
feeds the MiniBooNE experiment. Protons from the Booster strike a 71 cm beryllium 
target inserted in a magnetic focusing horn. Protons arrive at this target in 1.6 /is 
long Booster spills. The timing structure within each spill delivers 84 2 ns wide 
bunches of beam, each separated by 18 ns. Secondary short-lived hadrons (primarily 
pions) produced in the target are focussed by the horn and enter a decay region. In 
normal MiniBooNE operation, this decay region is 50 m long, at the end of which 
region is a beam absorber to stop hadrons and low energy muons. Located 25 m from 
the proton target is an intermediate absorber which can be lowered into the beam 
for use as a systematic check on the MiniBooNE v e background from fi decays. It 
is assumed that the 25 m absorber will be in place during the period when FINeSSE 
is operational to accommodate FINeSSE and MiniBooNE physics goals. Figure 13.11 
provides a diagram of the two possible absorber positions. 

The neutrino flux resulting from this design was simulated using the same tools 
currently being employed by the MiniBooNE collaboration The beam simulation 
utilizes GEANT 4 transport code ^Sj, and the MiniBooNE JAM pion production 
model which includes all beamline elements (horn, shielding, absorbers, etc.) and 
n^, K^, K° production from proton interactions on beryllium. To better reproduce 
the energy distribution of neutrino events observed in the MiniBooNE detector, pion 
spectra were input from a Sanford- Wang-based global fit [JO] to pion production 
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data in the relevant energy range in a procedure similar to that adopted by K2K. 
Figure E21 shows the resultant muon neutrino flux expected from a 25 m decay length 
beam produced at the 100 m FINeSSE detector site. In this configuration, 9.55 x 10~ 9 
muon neutrinos per POT per cm 2 are anticipated with a mean energy of ~ 700 MeV. 
The neutrino flux is roughly 20 times larger than that expected in a comparable 
volume at MiniBooNE. Note that the flux is diminished by about a factor of 1.8 
in switching from a 50 m decay length to a shorter 25 m decay length. However, as 
will be demonstrated, the 25 m absorber location is ideal for optimizing FINeSSE's 
oscillation sensitivity. 

Figure 13.31 shows the individual contributions to the total neutrino flux expected 
at FINeSSE. Contaminations from I^'s and z/ e 's are predicted to be 6% and 0.5% 
of the total flux, respectively. Once the "wrong-sign" T7^ background events are 
weighted by their appropriate cross section, they will comprise less than 1.5% of the 
total events in the FINeSSE detector. 

Better knowledge of the incoming neutrino beam flux enables more precise cross 
section measurements at both MiniBooNE and FINeSSE. The Booster neutrino flux 
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Figure 3.2: Expected flux at a 100 m detector site assuming a 25m decay length. 

will be much more precisely known than the fluxes reported in previous low energy 
neutrino cross section measurements well in advance of FINeSSE's commissioning. 
This improved knowledge comes from two sources: data from the Brookhaven E910 
experiment jSH] and from the CERN HARP experiment jHI] • Analysis that is already 
underway of E910 proton-beryllium data taken at 6, 12, and 18 GeV beam energies 
will be instrumental in verifying the extrapolation of the Sanford-Wang parametriza- 
tion jlHl to the 8 GeV Booster beam energy. More importantly, HARP data taken 
at 8 GeV on the Booster neutrino production target slugs will provide a tighter con- 
straint on the flux. The high statistics HARP data will provide a statistical precision 
of ~ 2% on 7r + production, which is the main source of muon neutrinos at both 
the FINeSSE and MiniBooNE detectors. Therefore, with these additional inputs, the 
overall muon neutrino flux at FINeSSE should be known to roughly 5% [3]. 

3.2 Event Rates 

The number of neutrino events expected in the FINeSSE Vertex Detector is calcu- 
lated using the NUANCE Monte Carlo jHS] to generate neutrino interactions on CH 2 . 
NUANCE is open-source code originally developed for simulating atmospheric neu- 



3s 



Expected Fluxes at 100m with 25m absorber 



> 

O 




0.5 1 1.5 2 2.5 3 

E, (GeV) 

Figure 3.3: Predicted flux contributions at a 100 m detector site assuming a 25 m 
decay length. Muon neutrinos comprise 93.5% of the expected neutrino flux. 

trino interactions in the 1MB detector. NUANCE has since been further developed 
and is now used by the K2K, Super-K, SNO, MiniBooNE, and MINERvA collabo- 
rations. The neutrino interaction cross sections in NUANCE have been extensively 
checked against published neutrino data and other available Monte Carlo event gen- 
erators. In addition, the full NUANCE simulation has been recently shown to provide 
a good description of events in both the MiniBooNE detector and K2K near detector 
ensemble. 

For this specific use, NUANCE was modified to include the FINeSSE detector 
composition and geometry, as well as the incident neutrino flux at the 100 m detector 
site (Section l3.1.2j) . Using the input neutrino flux distribution, NUANCE predicts 
event rates, kinematics, and final state particle topologies that can subsequently feed 
hit-level GEANT detector simulations, or, as in this case, simply estimate the type 
and number of neutrino interactions expected at FINeSSE. 

Table O lists the NUANCE-predicted event populations at the 100 m FINeSSE 
detector site with the 25 m absorber in position. The table provides the expected 
rates per ton detector for 1 x 10 20 POT as well as the expected backgrounds from 
the VZ and v e content in the beam. In all cases, the event rates have been normalized 
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to the number of contained neutrino events observed in the MiniBooNE detector 
Roughly 1.3% (0.6%) of the total neutrino events result fromp^ (i/ e ) interactions in the 
detector. The dominant contributions to the total event rate result from quasi-elastic 
and resonant processes: 41% of the events are CC quasi-elastic [y^n — > fi~ p), 17% 
are NC elastic iy^N —> v^N; N = n,p), and 35% resonant single pion production 
(z/ M N — > iT (is^) N 7r) channels. 
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Table 3.1: Number of events expected at 100 m with a 25 m decay length for 1 x 
10 20 POT per ton detector and for the full requested FINeSSE running and detector 
(rightmost column). These predictions do not include final state effects in 12 C and 
assume 100% detection efficiency. 



A total of approximately 360,000 neutrino interactions can be expected at FI- 
NeSSE for the full request of 6 x 10 20 POT. This raw estimate assumes a 9 ton 
fiducial detector and 100% detection/reconstruction efficiency 
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Effect of Final State Interactions 

Because a large fraction of neutrino interactions at FINeSSE take place on carbon, 
the number of expected events will depend not only on the predicted neutrino cross 
sections and flux, but also on the final state interactions engendered by the local nu- 
clear environment. Particles produced via neutrino interactions in carbon nuclei will 
have a chance to reinteract before exiting the nucleus, and thus can vanish or change 
identity before being detected. Although the initial reaction might be a simple CC 
quasi-elastic interaction [y^ n — > \f~ p), the observed final state particles might include 
pions, multiple nucleons, low energy photons, or all of these combined. Examples of 
the types of nuclear rescattering that can distort the final state observed at FINeSSE 
include simple absorption, charge exchange (ir + n — > 7T°p, tt° p — > ir + n, tt° n — > 7r~p, 
ir~ p — > 7r°n), and both inelastic and elastic scattering. For example, consider the 
resonant interaction v^n — > fi~ pir . If the ir° is absorbed before exiting the carbon 
nucleus, the interaction will appear to be quasi-elastic v^n — > \i~ p. Hence, the pres- 
ence of such final state interactions demands accounting in our observed event rate 
calculations. 

Table E21 summarizes the final states expected at FINeSSE after using NUANCE 
to simulate re-interactions within carbon nuclei. The table defines "QE-like" , "NC- 
EL-like", and "NC-7r°-like" event categories, where these classes refer to final states 
that appear to be QE, NC elastic, or NC ir° events, respectively. Specifically, 

• CC QE-like: a CC event with a muon and any number of nucleons in the event 
(no 7r or K in the final state) 

• NC EL-like: a NC event with any number of nucleons 
(no fi, ir, or K in the final state) 

• NC 7r°-like: a NC event with any number of nucleons and a single ir° 
(no other 7r's, K, or \i in the final state) 

Comparison of Tables 13.11 and 13.21 reveal that the number of observed QE and NC 
elastic interactions increases by roughly 10 — 15% as a result of final state reinter- 
actions. This is largely a result of resonant processes where the final state pion is 
absorbed. By the same mechanism, the overall number of observed NC ir° events 
decreases by roughly 30% because the final state ir° is either absorbed or "charge 
exchanges" . The contributions are further differentiated by the number of neutrons 
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and protons produced. More than half of the events yield only a single nucleon in 
the final state. The non-negligible rate of NC 7r° events with no final state nucleons 
results mainly from coherent pion production processes where the nucleus remains 
intact. 
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Table 3.2: Number of events for 1 x 10 20 POT per ton detector after including the 
effects of final state interactions in 12 C. The event classes are further broken down to 
indicate the number of nucleons present in the final state (either 0, 1, or > 1 proton 
or neutron). 
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With these definitions, Table 13.31 lists the dominant contributions to each final 
state. Of the observed QE-like events, 86.7% are true QE interactions. Of the events 
appearing to be NC elastic in the detector, 85.8% are true NC elastic interactions. 
Of the events appearing to be NC tt°, 73.5% + 19.0% = 92.5% are true NC tt° 
resonant or coherent interactions, respectively. Therefore, under the assumption of 
100% reconstruction and detection efficiencies, the level of irreducible backgrounds 
from final state effects appears to be less than 15%. However, reconstruction and 
selection criteria may potentially amplify or reduce the effect of such backgrounds. 

Just as non-QE events can appear to be quasi-elastic in the detector (via pion 
absorption), the reverse can also occur, albeit at a much smaller rate. NUANCE 
predicts that less than 1% of QE (or NC elastic) events will fail to appear quasi- 
elastic. This results from the small probability that a proton will rescatter in the 
carbon nucleus and produce one or more pions, for example: 

p + p — > p + n + 7r + 
p + p — > n + p + 7r + 

P + P — > p + p + 7T° 

p + p — > p + p + n° + TT° 

P + P — > P + P + 7T + + 7T~ 

p + p — > n + n + 7r + + 7r + 

p + p — > p + n + ii + + 7T°, etc. 

The situation differs for NC ir° events. In this case, roughly 30% of true NC tt° 
interactions will not appear to be NC ir° events in the detector: ~ 20% of the true 
NC 7T° interactions have no final state 7r° (that pion is absorbed before exiting the 
nucleus), and ~ 10% instead contain a final state tt + or 7r~ due to charge exchange 
processes, 

7T° p — >• 7T + U 
71° U — > 7T~ p. 

Therefore, because a large number of the neutrino scatters occur on carbon, it is 
important that the FINeSSE Monte Carlo simulation include these secondary final 
state interactions. Such a model is provided by the NUANCE generator, and is used 
in all event simulations provided in this document. 
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nnal state 


contribution 


fraction (%) 


CC QE-like 


QE {v^n -> fi~p) 


86.7% 


CC QE-like 


CC 7T+ RES {y^N —> fi-Nir + ) 


9.2% 


LL QE-like 


CC n RES [y^ n — > /i p7r u ) 


2.3% 


r<r< piT7 HI™ 


IjL 7T LjUri [y^A — > /i ATX J 


1.0/0 


r^r 1 m™ 


CC 7) [y^n — > [i pi]) 


U.o/o 


C^C* CXW \\\rr* 


r^ri TWQ ( -, , AT v ..— V\ 

Did [Ufj, iv — > /i A J 


n 00/ 
U.z/o 


NC EL-like 


NC EL (u^N^v^N) 


85.8% 


NC EL-like 


NC 7T° RES iV — > z/ M iV tt°) 


7.6% 


NC EL-like 


NC 7T~ RES (i/^n — > v^pn^) 


2.3% 


NC EL-like 


NC vr+ RES -> ^n7r+) 


2.1% 


NC EL-like 


TV T S~~* 1 / \ T T /A A 0\ 

NC 7r u COH [Ufj, A^v^A 7r u ) 


-1 r\ fry 

1.9% 


NC EL-like 


NC 77 (z/ M n — ► v^nr)) 


0.2% 


NC EL-like 


NC DIS (i/ M JV -> i/^) 


0.1% 


NC 7r°-like 


NC tt° RES N^u^N tt°) 


73.5% 


NC 7r°-like 


NC 7T° COH {y^A -> z^Att ) 


19.0% 


NC 7r°-like 


NC DIS (z/^JV -> ^1) 


2.4% 


NC 7r°-like 


NC 7T _ RES (i^n — > v^pix") 


1.8% 


NC 7r°-like 


NC vr+ RES -> ^nvr + ) 


1.5% 


NC 7r°-like 


NC EL JV - ^ TV) 


1.3% 


NC 7r°-like 


NCr? faN^UpNri) 


0.3% 


NC 7r°-like 


NC multi-7r iV -> v» N vr° tt°) 


0.2% 



Table 3.3: Fractional contributions to each observed final state at a FINeSSE 100 m 
detector site. "RES" ("COH") refers to resonant (coherent) pion production pro- 
cesses; "DIS" refers to deep inelastic scattering. 
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Chapter 4 

The FINeSSE Detector 



The FINeSSE detector is a 13 ton (9 ton fiducial) active target, consisting of a track- 
ing scintillator detector followed by a muon rangeout stack comprised of scintillator 
planes interspersed with iron absorber. The physics goals of this experiment require 
the ability to tag both v^p — > v^p and v^n — > pTp interactions by looking for the final 
state protons and muons produced in these channels. Proton energy and angle are 
measured in the first part of the detector, called the Vertex Detector. Muon tracks 
are tagged in both the Vertex Detector and the downstream Muon Rangestack. The 
Vertex Detector is also ideal for cross section measurements (such as single ir° produc- 
tion), which require good final state particle separation and good energy resolution. 
FINeSSE is designed to meet these requirements with a novel, yet relatively simple, 
detector. 



4.1 Detector Design and Construction 

The layout of the FINeSSE detector can be seen in Figure 14.11 The upstream Ver- 
tex Detector contains a wavelength-shifting (WLS) fiber array situated in a large, 
open volume of liquid scintillator. The downstream Muon Rangestack is comprised 
of 4.1cm x 1cm scintillator strips organized into planes in alternating X and Y ori- 
entations, interspersed with iron absorber. The Vertex Detector is described below 
in Section 14.1.11 this is followed by a description of the Muon Rangestack in Sec- 
tion SJJ 
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Figure 4.1: A schematic drawing of the FINeSSE detector. The cubic volume is the 
(3.5m) 3 Vertex Detector. It is a (2.4m) 3 signal region surrounded by a veto, filled 
with scintillator oil. The larger volume downstream is the Muon Rangestack. 
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4.1.1 The Vertex Detector 

The FINeSSE Vertex Detector consists of a cube-shaped volume of liquid scintillator 
with dimensions 2.4 x 2.4 x 2.4 m 3 . Light generated by ionizing particles traversing 
the scintillator is picked up by 1.5 mm diameter, WLS fibers, submerged throughout 
the sensitive volume. The fibers are mounted on a support frame, and are connected 
on one end to multi-anode photomultipliers, mounted to the outside of that frame. 
The fiber frame, photomultipliers, and associated electronics form a unit; this unit 
is immersed in the liquid scintillator, which is contained in a cubic tank, 3.5 m on a 
side. The volume between the fiber structure and the tank wall is used to monitor 
charged particles entering and exiting from the tracking volume ("veto shield"). The 
photomultiplier signals are processed in situ and transmitted by Ethernet to the 
outside of the tank, thus minimizing the number of cables that penetrate the tank 
wall. A schematic drawing of the tracking detector is shown in Fig. 14.21 Cables 
penetrate the tank wall above the oil level to prevent leaks. 

Particle tracks can be reconstructed because the time of arrival of the light at the 
end of a fiber from a given source inside the detector is a known, continuous function 
of the distance between the source and the fiber. The detector will be calibrated 
using cosmic rays. 

The arrangement of the WLS fibers is shown schematically in Fig. 14.31 There 
are three sets of fibers, running parallel to the axes of a Cartesian coordinate system. 
Except for a rotation in space and an offset, the three fiber sets are identical, consisting 
of fibers that intercept the wall at the vertices of a quadrate grid. The distance 
between grid points is 30 mm. Thus, the closest distance between any two fibers in 
the full assembly is 15 mm. The resulting arrangement is invariant with respect to a 
rotation by 90° about any major axis. For the given dimensions, there are a total of 
80 x 80 x 3 = 19200 fibers. 

Tracking Scheme 

Consider a point source of light at some distance d from a long, WLS fiber. A fixed 
fraction of the light that intercepts the fiber is wavelength-shifted and transported to 
the photo detector. Ideally, the exiting light is proportional to the diameter of the 
fiber and inversely proportional to d. In reality, this distance dependence is faster 
than 1/d, because of light attenuation in the scintillator, and because the fraction of 




Figure 4.2: A schematic drawing of the Vertex Detector shown from the side. 
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Figure 4.3: The geometrical arrangement of WLS fibers inside the Vertex Detector. 
The arrangement consists of three orthogonal sets of parallel fibers. The geometry is 
symmetric with respect to a rotation by 9(P about any of the three major axes. 
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light captured by the fiber depends on the angle between the fiber and the incident 
ray. 

The light collection efficiency versus distance to the fiber can be determined by 
Monte Carlo and test measurements. This dependence will be the same throughout 
the detector volume, with the possible exception of regions close to the wall, where 
reflected light may contribute. 

The light from the source travels to all nearby fibers through a completely homo- 
geneous medium. Since the sharing of light by nearby fibers can be used to localize 
the source, the position and angle of tracks can then be reconstructed. Fibers along 
a given direction are only sensitive to the projection of the track onto a plane per- 
pendicular to that direction. Even if a second, orthogonal set of fibers is provided, it 
is still possible that a track may be parallel to one of two directions. This difficulty 
is avoided in the FINeSSE detector by having three orthogonal sets of fibers. 

Light Generation and Transport 

A comparative study of different scintillator fluids used in conjunction with WLS 
fibers is given in Ref. jHl]. In general, an ionizing particle excites ultraviolet fluo- 
rescence (~ 350 nm) in the scintillator. This light is normally shifted to a longer 
wavelength (~ 430 nm), to avoid problems with absorption in the scintillator. The 
shifted light propagates isotropically. The light that intercepts a WLS fiber is shifted 
once more, to typically ~ 500 nm to inject light into the acceptance cone of the fiber 
and to prevent re-absorption in the fiber. The WLS fiber consists of a polystyrene 
core (n=1.60), surrounded by cladding of a lesser index of refraction. The trapping ef- 
ficiency of the fiber is significantly enhanced (to typically 6%) when two claddings are 
used. For a 1.5 mm diameter fiber, the first cladding is 45 fim thick acrylic (n=1.49), 
and the second is a layer of 15 fim thick fluor-acrylic (n=1.42). The data given here 
are for BCF-91A-MC fibers from Saint-Gobain [S3]. Similar fibers are available from 
Kuraray (SH]. The second cladding also provides protection against possible chem- 
ical interaction between the liquid scintillator and WLS fibers. Long-term tests of 
fibers in mineral-oil-based scintillation fluid, carried out in the context of MINOS 
R&D (S3 EH , showed no discernible ill effects. Specifically, the BCF91A fibers used 
in these studies were not affected after having been suspended for six months, at 
temperatures up to 50° C, in mineral-oil-based scintillator, BC517L, or in the high 
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fluor concentration BC517H. Furthermore, single clad fibers suspended in B517L for 
more than two years also showed no aging or deterioration 

The attenuation length of the light propagating in the fiber is given as 3.5 m by 
the manufacturer, but a more complicated behavior is reported in the literature |54j . 
The attenuation length is the same whether the fiber is in liquid scintillator or in 
air jSH]- Because of attenuation in the fiber, the light collected at one end depends 
on the point of origin along the fiber. This effect can be reduced by applying an 
aluminum reflective coating at the other end jSHl, or even by just painting it white 
with Ti0 2 jSH- These coated fiber ends will be covered with a Teflon sleeve and 
will not be in contact with the mineral oil. Even in the case that they would, these 
coatings are inert and unlikely to interact with the scintillator oil. 

We are currently investigating a number of options concerning the liquid scintil- 
lator. It is true that the present tracking scheme makes use of the sharing of light 
from a given source by a number of fibers. However, it is also true that only fibers 
in the vicinity of the track contribute significantly to the determination of the track 
parameters. As a result, low- attenuation length oil may be the best choice, which 
is an unusual but easily accommodated need. We are continuing R&D on the best 
choice of scintillator oil. 

Advantages of the Proposed Design 

Our design represents a novel approach to the task of tracking ionizing particles 
in a large-volume detector. It exploits the fact that the response of a fiber versus 
the distance to the light source inside the detector volume is a universal function, 
and the three-fold symmetry makes the tracking sensitivity nearly isotropic. These 
properties are particularly important for the physics goals of the present proposal. 
To our knowledge, there is no other scheme that offers these features. 

We investigated an alternative detector design which consists of planes of bars of 
solid scintillator, oriented normal to the beam direction (z). Bars in even layers would 
run in the x direction, those in the odd layers in the y direction. Scintillating plastic 
bars are extruded polystyrene with 2cmxlcm cross section. Each bar is co-extruded 
with a TiC>2 outer layer for reflectivity and a hole down the center for a WLS fiber. 
This design is very similar to the K2K scibar detector recently commissioned [BUI and 
to that employed by the MINOS experiment [fi2~] . 
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Such a design for the Vertex Detector has a number of disadvantages for the 
FINeSSE physics goals. First, the tracking information would have to come from 
light sharing between several intercepted bars. The track resolution would therefore 
be limited by the bar dimension. Second, the amplitude information from a given 
bar reflects the overlap of the track with the bar cross section. Interpretation of 
this information is complicated and suffers from irregularities in stacking, caused by 
uneven extrusions and variations in the reflecting layer and in wrapping. Third, tracks 
that do not have a sufficiently large angle with respect to the stacking planes are lost. 
These disadvantages make reconstruction of low energy proton tracks and therefore 
low energy v — p elastic events particularly difficult. Details of our studies of this 
design can be found in Appendix 1X1 

Prototype Setup and Test Measurements 

To demonstrate the viability of the proposed tracking scheme in the Vertex Detector, 
we have constructed and tested a prototype. The following contains a description of 
the prototype including construction issues and the results of beam tests. 

The prototype Vertex Detector consists of a rectangular box made from anodized 
aluminum, 16 x 16 x 30cm 3 on the inside (Fig. 14.4)) . with a 6 x 5 array of multi-clad, 
1.5mm diameter, WLS fibers (Bicron BCF-91A-MC, [53]). The fibers are spaced 
20 mm apart, and penetrate the wall of the box. The fibers' O-ring seals hold them 
in place, as will be done in the full-scale Vertex Detector (see Section |4~3*]) . Light from 
the fibers is detected on one end by two 4x4- anode photomultipliers (Hamamatsu 
H8711). The light-tight box is filled with liquid scintillator (Bicron BCS517H [S3]). 

To test this prototype, the box was placed at the end of beam line I of the Radi- 
ation Effects Research Program (RERP) test station at IUCF [SJ. A low-intensity, 
200 MeV proton beam with a 6 x 6 mm 2 profile (defined by two scintillators in coin- 
cidence) was aimed at the center of the box. Each passing proton produced a track 
with almost uniform energy deposition. The box could be moved vertically (hence- 
forth called x-direction), or horizontally (y-direction) by remote control, enabling the 
beam position to be scanned perpendicular to or along the length of the fibers. The 
amplitude signals from all 30 fibers were digitized by conventional electronics and 
stored, event-by-event. 



53 




Figure 4.4: Photograph of the prototype detector with the mounted WLS fibers. The 
protruding fiber ends were covered with a light shield. The photomultipliers were 
mounted on the plate at the far end. The box could be pivoted to study beams that 
intercept the fiber matrix at an angle. 
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Light Yield Versus Track-to-Fiber Distance and Angle 

With the box axis normal to the beam, a sequence of runs were taken, moving the 
beam from top to bottom in steps of 5 mm. For each fiber and each run, the centroid of 
the accumulated amplitude spectrum, corrected by its pedestal value, was determined. 
The centroid as a function of beam position is shown in Fig. 14.51 Before the centroid 
values were plotted, each fiber output was scaled by an individual gain factor such 
that all peaks lined up at an arbitrary centroid value of 100. The 30 fiber pedestals 
and gains thus determined were retained for the remainder of the test. This part of 
the test demonstrates that there is a universal law that governs the collected light as 
a function of the track-to-fiber distance. 

In order to study the angular dependence of incident light, the box was rotated 
around a vertical axis by an arbitrary angle (j3 = 27°), and the vertical scan, described 
in the previous section, was repeated. The result was virtually indistinguishable 
from Fig. 14. 5[ except for an overall scaling by -^s- This shows that the amplitude 
response of the detector is indeed sensitive to the projection of the tracks onto a plane 
orthogonal to the fiber direction. 

Attenuation of Light in WLS Fibers 

With a normal beam at constant height (x = 0.5 cm), a scan along the fiber direction 
was carried out. For all fibers, we observed a light yield that slowly increased toward 
the end of the fiber that was coupled to the photomultiplier, with a sharp fall-off near 
both ends of the box. From the overall increase, we estimated a light attenuation 
length in the fiber of 150 to 200 cm. The fall-off at the ends is well explained by 
purely geometric arguments. 

Absolute Light Yield 

To absolutely calibrate the light yield, it is necessary to relate the photomultiplier 
amplitude to the number of collected photoelectrons (PEs). This can be done if 
the peaks that correspond to events with only one or two PEs are resolved. In the 
spectra obtained with the multi-anode photomultipliers this was not the case. We 
have therefore compared the tubes used in this test with a high-performance tube, 
using low-light pulses. We found that a track that intercepts the fiber (corresponding 
to 100 units in Fig. 14 .5j) yields an average of 17 ± 2 PEs. 
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Figure 4.5: Normalized centroids of the light yields of 15 fibers as a function of the ver- 
tical beam position. The insert shows the position of the fibers plotted. The light yield 
is given in arbitrary units. To convert to units of observed photoelectrons , multiply 
by 0.17 (see text). 
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Reconstruction of the Track of a Single Event 

In the following sections we describe an attempt to use a simple reconstruction algo- 
rithm to extract the angle and position of a track of a single event, and to determine 
the corresponding uncertainty. The track reconstruction algorithm used here is prim- 
itive but sufficient to prove the point. For instance, the simple algorithm does not 
take into account that the energy loss along the track through our test setup changes 
by about 40%. The track reconstruction algorithm used for the Monte Carlo sample 
for the entire FINeSSE detector, described in Chapter accounts for these effects. 

The observed dependence of the amplitude, A, on the fiber-to-track distance, 
d, (Fig. I4.5J) . is used to construct a function A(d), chosen so it could be inverted 
analytically. Then, the amplitude data are converted to measured distances. In other 
words, for any given event i, the distances d ni of that track from all fibers {n = 1 — 30) 
are known. 

Fig. 14. 61 shows the average distance of tracks from the top row, versus the average 
distance from the bottom row. Each event is represented by a point. The plotted 
events are comprised of 500 events from each of the runs of the vertical scan discussed 
above. As a result, the tracks to be analyzed more or less fill the space between the 
top and bottom rows. 

Since the sum of the two distances equals the spacing between top and bottom 
row (80 mm), the events are expected to lie on a straight line. For the same events, 
Fig. 14.71 shows the average distance of tracks to the second row versus the average 
distance of tracks to the third row. For the tracks between the rows the sum of the 
distances is constant, as in Fig. 14. 61 but for those outside the rows the difference must 
be constant. This can clearly be seen in the figures. 

We parametrize a track by its position at the downstream face of the detector 
(the intercept with the x-axis) and its angle with respect to the z-axis. For a given 
event (track), the distances to all fibers are deduced from the measured amplitudes 
via the function A(d), mentioned earlier. To reconstruct the track, one needs to find 
the straight line that is most consistent with these distances. This is accomplished by 
a regression algorithm minimizing the square of the difference between the calculated 
and the measured distance, weighted by the amount of light, summed over all fibers. 
In this way, one obtains the position x and angle a for each individual track. 

We have applied this procedure to 500 events of a run where the beam was normal 
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distance track to bottom row (mm) 



Figure 4.6: Distance to top row versus distance to bottom row, averaged over all six 
fibers in that row, for horizontal tracks that fill the space between the two rows. 

to the fibers, and aimed halfway between rows 2 and 3. The expected position is 
therefore x=30 mm, and the expected angle is a=0. In the reconstruction code it 
was possible to decide whether a given fiber should be allowed to contribute in the 
fitting procedure. The two panels in the top row of Fig. 14.81 show the reconstructed 
position and angle distributions when all fibers in the two rows on either side of the 
nominal beam position (a total of 24 fibers) are taken into account. The distributions 
represent a convolution of the intrinsic detector resolution and the phase space of the 
"beam," which is defined by the trigger scintillators. In order to illustrate how this 
result depends on the number of participating fibers, we have repeated this analysis, 
taking into account only three fibers in each of the rows straddling the beam (a total 
of 6 fibers). The result is shown in the bottom row of Fig. 14.81 

The measured track position and angle resolutions are dominated by statistical 
fluctuations in the number of PEs observed in the photomultiplier. This can be seen 
in the following way. Assuming a hypothetical track, and using the light-distance rela- 
tionship d(A) mentioned earlier, the signal amplitude (in arbitrary units) is calculated 
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distance track to tow # 3 (mm) 

Figure 4.7: Distance to row 2 versus distance to row 3, averaged over all six fibers in 
that row, for horizontal tracks that fill the space between top and bottom row. 

for all fibers. From the absolute yield calibration (100 amplitude units correspond 
to 17 PEs) we get the corresponding number of PEs. From a Poisson distribution 
with that mean, a random number is drawn to represent the "observed" number of 
electrons. This number is then converted back to amplitude units for all 30 fibers, 
and the resulting simulated event is reconstructed by the same algorithm used for 
the real events. As input, random tracks with a = and 26.5 mm < x < 33.5 mm 
(representing the beam width) were used. The simulated distributions for x and a 
are shown in Fig. 14.91 Obviously, they are remarkably similar to the real distributions 
(Fig. I4.8j) . and since statistics is the only reason for smearing in the simulation, we 
conclude that this is also the case for the measured data. 

Given the good agreement between real and simulated data, we can deduce the 
intrinsic detector resolution. We find that the position resolution achieved in our test 
is about 15 mm FWHM, and the angular resolution about 15° FWHM. The difference 
between using 24 fibers and using only six fibers is not solely due to counting statistics. 
It is also affected by the fact that the length of the track, that is supported by distance 
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Figure 4.8: Reconstructed position and angle of tracks for two rows of fibers on either 
side of the beam (top plots) and for only three fibers on either side of the beam (bottom 
plots). The rectangular box in the left panels indicates the nominal beam position and 
width. 
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Figure 4.9: Simulated events, taking into account the Poisson statistics of the observed 
number of PEs (cf. Fig. \4-&\) - 
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information from the fibers, changes, thus relaxing the constraint on the track angles. 
4.1.2 The Muon Rangestack 

Downstream of the FINeSSE Vertex Detector is the Muon Rangestack. The Ranges- 
tack is designed to range out muons, in conjunction with the Vertex Detector, with 
energies up to 1.5 GeV, and to measure muon energy to 10%. These characteristics 
are needed to enable the reconstruction of charged current events, which are used 
in both the As and the disappearance measurements. The Rangestack is located 
downstream of the Vertex Detector, because the muons from high energy CC events 
in the Vertex Detector tend to be produced at forward angles (See Section 15.2.21) . 

The iron absorber planes and tracking granularity in the Rangestack are designed 
to meet these requirements. The stack is comprised of 4 x 4 m 2 planes of scintillator 
strips and iron absorber with an overall depth of 0.85 m (0.98 m including support 
structure) in the beam direction (referred to here as the z direction), and a weight of 
100 t. Of the 0.85 m thickness, 0.24 m is scintillator and 0.61 m is iron. The design 
both meets the physics requirements and minimizes cost and space demands. The 
following sections describe the Rangestack's design, construction, and readout. 

Figure 14.101 shows isometric and orthographic views of the entire detector. The 
structural steel shown supporting the stack is a conceptual design. No PMT enclo- 
sures or clear fiber bundles surrounding the Rangestack are shown. 

Advantages of Choosing a Rangestack Design 

We have found that a Rangestack is the best choice to achieve our physics goals. 
The Rangestack is designed to range out muons, in conjunction with the Vertex 
Detector, with energies between 0.1 and 1.5 GeV. A stack of plastic scintillator bars 
interspersed with iron absorber is a relatively inexpensive and straightforward way 
to achieve these goals. This detector design has been well tested by the MINOS 
collaboration. They have shown that these plastic scintillator strips have excellent 
reflectivity, fast timing, simple design, long term stability, low maintenance and high 
reliability jH2]- Furthermore, the expertise for construction and assembly of such a 
detector already exists at Fermilab. 
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Figure 4.10: An isometric rendering of the Muon Rangestack in front of a box bound- 
ing the envelope of the Vertex Detector. 

A magnetic spectrometer was considered and rejected because it is very expensive 
compared to the Rangestack. Its obvious advantage - charge identification - is not 
necessary for the physics of the experiment. 

Rangestack Design 

The Rangestack consists of 21 pairs of a scintillating tracking plane and iron absorber. 
Each scintillating tracking plane is made up of 4.1 cm x 1 cm x 4 m extruded 
polystyrene scintillator strips. Each plane is oriented normal to the beam, with strips 
in the x direction, in even numbered planes, and in the y direction in odd numbered 
planes. The strips are arranged so that their depth in the z direction is 1 cm. Each 4 m 
plane, then, contains 96 strips. The strips are packaged in groups of 16, surrounded 
by a 1 mm thick aluminum can for protection and light tightness. There are six cans 
of strips per x or y layer. Cans are identical between x and y layers, just rotated 90° 
with respect to each other. The six cans in a layer are attached to a backing layer 
of 0.50 inch thick steel which serves to support the 4x4 m 2 plane of scintillator and 
protect it during crane handling and shipment, and functions as part of the absorber 



63 




Figure 4.11: A completed x layer with six cans of scintillator and a 0.50 inch thick 
steel backing. The assembly weighs 2.2 tons. Cans are attached to the steel backing 
by bolts at the ends of the can, and a strap ( not shown ) similar to the MINOS design 
which prevents the cans from bending away from the backing in the middle. 

layer between each x and y layer of scintillator. Figure 14.111 shows an assembled x 
plane. Assembled y planes are similar, with scintillator planes rotated by 90°. 

Scintillator strips in these planes are similar to those used in the MINOS experi- 
ment jH2]- The scintillator is co-extruded with a TiC>2 outer layer for reflectivity. A 
groove down one side of the strip contains a glued-in-place WLS fiber. Light produced 
in the scintillator strips is absorbed by the WLS fiber and re-emitted at ~ 500 nm. 
Re-emitted light within the acceptance for total internal reflection is transmitted to 
the end of the fiber and detected by the phototubes, which are identical to those used 
in the Vertex Detector. The same type of WLS fiber will be used in the Rangestack 
and the Vertex Detector (see Section 14.1.11 for a description) . 

Interspersed with the active scintillator planes are iron absorber planes. In order 
to minimize the size of the Rangestack, the thickness of the iron in the absorber 
planes increases with increasing depth into the detector. This design is chosen to take 
advantage of the fact that once lower energy particles range out (leaving only higher 
energy particles), the iron absorber thickness can be increased while maintaining the 
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same relative energy resolution. This is achieved easily by dividing the Rangestack 
into four sections with the thickness of the beam absorber increasing by 0.5 inches, 
or one iron plane, per section. 

The upstream end of the range stack (section 1) has six layers, three x and three 
y, with 0.50 inch total steel between layers. The backing plate for the scintillator 
planes provides this layer of absorber. Behind that is section 2, with six layers of x 
and y planes separated by a total of 1.00 inch of steel between modules. To make up 
the total thickness, there are 0.50 inch thick filler plates between scintillator modules. 
Section 3 has six layers of x and y scintillator planes separated by 1.50 inches of steel. 
Section 3 filler modules are 1.00 inch thick. Finally, section 4 has three absorber layers 
2.00 inches thick, using 1.50 inch thick absorbers. An x layer sits upstream of the 
first 2.00 inch thick absorber, followed by alternating layers of y, filler and x layers. 
The stack ends with a y layer. The final x and y layers in the stack are triple-thick 
to assure high efficiency for identifying stopped muons. Figure 14.211 shows a typical 
1.50 inch thick filler module for section 4. The minimum thickness of 0.5 inches used 
in the first section of the Rangestack is set so as to ensure a 10% energy resolution 
for the lowest energy muons that enter the Rangestack. 

Light Output and Readout in Active Rangestack Planes 

Tracking planes in the Muon Rangestack must produce enough PEs at the readout 
PMTs to observe a Minimum Ionizing Particle (MIP) and therefore track muons. 
The scintillating strips, WLS fiber, and readout used in the Muon Rangestack are 
very similar to those used in the MINOS far detector. Light output from a Mini- 
mum Ionizing Particle (MIP) in the FINeSSE Muon Rangestack can be inferred from 
prototype testing results from MINOS. 

As MINOS studies have shown, the efficiency for observing a muon crossing a 
strip is greater than 90% as long as at least 2.5 PEs are read out. In a prototype 
assembly of part of a scintillator plane for the MINOS far detector, a MIP passing 
through the center of a strip produces, on average, 6 PEs in total, with an average 
of 3 PEs read out on either side jH21- The design of the Muon Rangestack is such 
that there will be approximately twice this much light observed as a result of the 
differences between the MINOS strips and the FINeSSE Muon Rangestack strips. 
These differences in the design include the overall length of the strips, the WLS fiber 
size, and the fact that the readout takes place from one end only. The first two factors 
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increase the light levels in the Rangestack by more than a factor of two as compared 
to the MINOS detector. The fact that readout occurs from only one end will decrease 
the light read out, but only slightly as reflective paint on the terminal end of the fiber 
will pipe light back toward the readout end. 

WLS fibers are connected to clear fibers at the end of the scintillator strips as is 
done in the MINOS detector and as is shown here in the Muon Rangestack design. 
Clear fiber ribbon cable will then route the signals to the removable floor above the 
detector where they will be read out via the same system employed in the Vertex 
Detector as described in detail in Section 14.1.31 

4.1.3 Signal Readout: Phototubes and Electronics 

We propose a common design of signal readout to be used for the Vertex Detector, its 
veto shield, and the Muon Rangestack. Requirements for the readout system include 
independent amplitude and timing measurements for each channel. The amplitude 
information is used in the track reconstruction and for dE/dx measurement for the 
particle tracks; timing information is used to assemble the hits (rejecting spurious 
noise hits and background tracks), to correlate with the beam spills, and in the 
analysis of secondary events, such as muon decay and nuclear decay in the active 
volume following the capture of a neutron produced in the primary event. In the 
Vertex Detector the full scale signal range will extend to > 50 PEs, with an amplitude 
resolution of < 0.1 pe. The timing resolution is not crucial, but it needs to be ~ 10 ns 
to reduce spurious hits and background tracks. 

The readout must be sensitive both to the primary interaction events (which are 
in time with the beam spill) and to secondary events as described above. It is also 
useful that it be sensitive to cosmic ray muon tracks, for calibration purposes. For 
these reasons, the front end electronics must be self-triggering and have a relatively 
low deadtime after being hit. Deadtime of the proposed system is about 1 — 5 /is, 
depending on signal amplitude (see discussion below). The deadtime applies inde- 
pendently to each channel - this is a matter of recovery of the front end electronics 
only, since there is zero deadtime associated with the data readout. A external global 
trigger will also be implemented, which can be used to force an event readout and/or 
a charge pulser event, for diagnostic purposes. 

The readout system designed to meet these requirements combines multianode 
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photomultiplier tubes (MAPMTs) and custom readout electronics in a 128-channel 
"readout module." This includes two MAPMTs, front end electronics, HV bias cir- 
cuits, charge-injection test pulsers, and data acquisition and control communications. 
The module is a sealed, conductively cooled unit suitable for submerged operation 
inside the oil tank (for the Vertex Detector application). In the Muon Rangestack 
the readout modules will be housed above the removable floor above the detector. 
The Vertex Detector requires 150 modules for readout of the active volume and six 
modules for the veto shield. The Muon Rangestack readout requires an additional 17 
modules. 

Multianode PMTs 

Fibers from the Vertex Detector, veto, and Muon Rangestack are read out via Hama- 
matsu R7600-00-M64 MAPMTs. This multianode photomultiplier tube provides an 
8x8 array of optical readout in a 25.7 mm square by 20.1 mm deep metal pack- 
age. We intend to operate the MAPMTs at a gain of 8 x 10 5 , which corresponds 
to a typical cathode supply voltage of —875 V. Channel-to-channel gain variation 
in the MAPMT is 5:1 maximum. As a result, if it is assumed that a channel in the 
middle of this ratio is set to the desired gain, the required full scale range of the 
front end electronics is 8 x 10 5 x \fb x 50g = 14 pC, and the required resolution is 
(8 x lO 5 /^) x O.lg = 5.7 fC, a dynamic range of 2500. The R7600-00-M64 MAPMT 
is similar to the sixteen anode R5900-00-M16 (see in Figure E.13j) . 

Front End Electronics 

Continuously running commercial multichannel 10-bit pipeline A/D converter chips, 
such as the ST Microelectronics TSA1005-20IF, will form the heart of the FINeSSE 
front end electronics (FEE). One such converter will digitize signals from each anode 
of each MAPMT, at a rate of 20 Msamples/s. Since this sampling is, of course, 
significantly coarser than the MAPMT output pulse width, and significantly coarser 
than the desired time resolution of 10 ns, suitable analog signal processing before 
digitization, and digital signal processing after, will be employed. 

The analog signal processing must shape the MAPMT output pulse so that it 
extends over several A/D converter samples (offering a chance to measure the ampli- 
tude), and so that it does something characteristic (such as have a peak, or better yet 
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a zero-crossing) at a defined time after the MAPMT output pulse (offering a chance 
to measure the time of occurrence). In addition, the analog signal processing circuit 
must present a low and resistive impedance to the MAPMT anode, so that all the 
induced charge is swept quickly into the signal processing for that anode, and does 
not go into parasitic capacitances to ground or to other anodes. The simplest circuit 
which can achieve these things is an inverting op-amp with an RLC network in feed- 
back. In response to a ^-function current input, I(t) = Q5(t — t ), a damped cosine 
wave is output, V(t) = (Q/C) cos(u(t — to)) exp(— \(t — 1 )) for t > to- The amplitude 
is proportional to the input charge, and the phase, or equivalently the time of the 
first zero crossing, is directly related to the time of occurrence of the input charge. 

The required digital signal processing simply undoes the above. If R, L, and 
C are known, then uo and A are known, and a simple linear fit of the measured 
data points to V (t) = A cos(u;(t — t' )) exp(— \(t — t' )) + B sin(cj(t — t' )) exp(— \(t — 
ip)), where to * s the time of the first ADC sample found above threshold, determines 
the amplitude and the relative time to — t' of the input pulse. Of course, there will 
be manufacturing tolerances on the component values, and potentially a small drift 
with time or temperature. To calibrate these, a nonlinear fit involving uo and A as 
parameters is performed occasionally. 

The number of raw A/D samples to be used, and the desired value of u and A of 
the filter, needs to be carefully optimized. In particular, the response should decay 
to less than a few percent within 5 /xs. A following pulse which occurs within 5 /is 
can still be measured, provided that the amplitude of the waveform does not exceed 
the range of the A/D converter. The actual deadtime for a full scale amplitude pulse 
following another full-scale amplitude pulse may be as large as 5 /is, but for pulses of 
smaller amplitude, e.g., less than half-scale, the deadtime can be as small as 1 /is or 
even 800 ns. For pulses separated by less than that, it is still possible to discriminate 
that there is a second pulse, by applying a threshold to the x 2 °f the fit, even if it 
is not possible to estimate the amplitudes and times separately. The low deadtime 
and double-pulse discrimination capability is important to detect muon decay events, 
both for particle identification and as a method of gain calibration. 

An important difference from the conventional approach, where a discriminator 
and dedicated time-to-digital converter are used, should be noted. Conventionally, 
the discriminator is set to respond to the first, or at least the first few, PEs in the 
signal, and thus the measured hit time is insensitive to the later PEs. As a result, 
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the measured hit time is more accurately related to the time of the charged particle 
passing in the detector. In the present approach, where a timing discriminator is 
not used, by contrast, the measured hit time depends on all the PEs in the signal. 
This will not yield as good a result; nevertheless, because the scintillator oil and 
wavelength shifting fiber both have short decay times (2 ns for the oil, 3 to 10 ns for 
the WLS fiber), we can achieve the required 10 ns timing accuracy. 

The digital signal processing will be performed by a combination of field pro- 
grammable gate arrays (FPGAs) and a microprocessor (also used for communication 
and control - see below). The FEE FPGAs, handling 8 or 16 channels each, will 
monitor the digitized data from each channel, apply the trigger threshold, capture 
the required number of data points after seeing the threshold crossed, and send them 
to the microprocessor. There the linear fit calculations are performed, and the am- 
plitude and time thus determined are queued for transmission to the DAQ. 

This front-end electronics design for FINeSSE is intended to take maximum ad- 
vantage of low cost commercial electronics to provide a highly parallel measurement 
of pulse amplitude and timing with minimal deadtime. Although there are ASICs 
such as the VA-TA series from IDE, Inc. which have been applied to very high density 
readout electronics for MAPMTs, none of the presently available ASICs can match 
the low deadtime performance of the system described here when the hit patterns 
are non-sparse as in the FINeSSE Vertex Detector. Furthermore, FINeSSE does 
not strictly require such very high density electronics, since each MAPMT covers a 
240 mm x 240 mm area on the detector, making almost that much area available 
for the readout electronics. In any case, with the large commercial market driving 
down the costs of pipeline ADCs and FPGAs, the proposed approach can be more 
cost effective. 

An independent charge-injection test pulser will be included on each channel. 
This diagnostic feature will be useful in testing the front end electronics and the 
readout, and debugging the system during operation. 

The Readout Module 

The readout module is comprised of two Hamamatsu R7600-00-M64 MAPMTs de- 
scribed in Section 14.1.31 a high voltage bias generator, 128 channels of front-end 
electronics, a microprocessor, a data network interface, a system clock and trigger 
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Figure 4.12: Block diagram of the readout module. 

interface, and a power supply. (See the block diagram in Fig. 14.121 ) The micropro- 
cessor handles control and housekeeping tasks, data transmission tasks, and also the 
fitting algorithm (described above) to calculate the hit time and amplitude based on 
the raw ADC samples captured in the FEE. Normally only the time and amplitude 
will be sent to the DAQ, although for calibration of the algorithm above, the raw 
samples will be sent on selected events. 

The module will perform zero suppression based on a programmable amplitude 
threshold. This threshold is independent of the trigger threshold used in the FEE to 
detect and capture an event. Note that both thresholds are applied digitally and can 
easily be set on a channel-by-channel basis, which can be helpful in dealing with any 
hot channels. 

The module can also be programmed to suppress data not satisfying a pattern cut. 
For instance, a requirement of > 2 hits within a 20 ns window will almost completely 
suppress MAPMT dark pulses, even with the trigger threshold set very low, without 
adversely affecting real events. More complex cuts, such as requiring a track through 
the contiguous fibers of the module, can be imposed as well, although a cut requiring 
data from other modules would need to be imposed in the DAQ computer, since the 
event data will not be communicated directly between the modules. (Each module 
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reads out an 8 x 16 array of fibers on one face of the detector.) Note that the front-end 
electronics will process every hit seen on every channel; both the zero suppression and 
any other cuts that may be programmed into the readout module act on digitized 
data before transmitting it to DAQ, and the purpose is only to reduce the volume of 
unwanted data. 

The output data from the readout module will be, for each hit, a 32-bit coarse 
timestamp from a counter, 8-bit fine time, and 12-bit amplitude. For data budgeting, 
we assume here that each of the latter will be embedded in a 16-bit word, which may 
include some extra "status bits." (If data volume or data rates become an issue this 
could be revisited; also the 32-bit coarse timestamp could be compressed easily, since 
the relative time between subsequent transmitted hit data in general will not require 
32-bit representation.) At 8 bytes per hit, the worst case event, where all channels of 
a module are hit, will generate 1.02 kbytes of data to be transmitted to the DAQ. An 
elasticity buffer of 16 kbytes or greater will be implemented in the module to allow 
for bursts in the incoming data and other network activity such as packet re-tries and 
control transactions. 

Estimated power dissipation of the readout module is 20 W (156 mW/channel). 
The power input will be from a single 24 to 48 V DC supply. Each module will 
have an individual power cable from a fused power distribution panel outside the 
tank. In this way a module with a catastrophic power short can be easily isolated 
to allow the others to continue operation. High-reliability steering diodes will direct 
the input voltage to one of two independent fully isolated DC/DC converters in the 
readout module, based on polarity. The outputs of these supplies will be connected 
through further steering diodes. In this way we obtain the benefits of two almost 
fully independent power systems in each module, with no overhead for additional 
cables or connectors. Fully isolated DC/DC converters are required anyway, so that 
the modules can be powered from a common DC source and still have equal ground 
voltages, which simplifies the data network requirements (no AC coupling or wide 
common- mode transceivers are required). The internal circuitry of the module will 
work largely from a single +2.5 V supply (about 4.5 A), with an auxiliary ±3.3 V 
supply for some functions, and another auxiliary —1050 V supply for the MAPMT 
bias circuits. HV will be regulated to each MAPMT individually through a high- 
side pass transistor driven by an analog optoisolator. This will feed a conventional 
resistive voltage divider for the dynode string. A Cockcroft-Walton type of PMT bias 
circuit will be evaluated as an alternative, for reduced power, but reliability is likely 
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to be greater using the resistive divider circuit. 

The microprocessor and the data network interface are also potential single-point 
failures affecting an entire module. We will carefully evaluate a redundant system 
design around the microprocessor and its boot EEPROM in particular. We expect 
to implement a redundant network (pair of networks) for each of the five logical 
networks (described below); in particular an unpowered module shall be positively 
disconnected from the data (and clock) networks. Even when powered, it will be 
positively disconnected from the primary data network unless a watchdog timer is 
kept alive by properly functioning software. Download of new code shall be possible 
via either network, even if the old contents of the boot EEPROM are corrupted. (In 
particular the network disconnection watchdog may be defeated until a defined time 
after power- up.) 

With proper attention to the above details, reliability of the modules can be 
assured to a degree which will allow them to be mounted in the inaccessible interior 
of the oil tank. The costs involved in this redundant design are significantly less than 
would be involved in transporting either the scintillation light out of the oil tank on 
optical fibers, or the anode signals from MAPMTs mounted in the oil tank to readout 
electronics outside. Remaining single-point failures, except in the MAPMTs or the 
divider resistor string, can affect at most eight channels. 

The power-dissipating electronic components in the readout module will conduct 
heat through their leadframes or other thermal bonds into the extra thick (2 oz.) 
copper ground plane of the printed circuit boards. The printed circuit board ground 
planes are directly tied to the the outer aluminum case of the module, either by wedge- 
lock clamps or simply by screws and aluminum frames. Temperature rise between 
the module case and the printed circuit board ground planes will be about 5 - 10 K. 
Scintillation oil circulating in the Vertex Detector will cool the readout modules. 
The total heat load into the oil is 156 x 20 W = 3.1 kW, which can be removed 
by circulating the oil through an external water-cooled heat exchanger at at least 
8 gal/min, assuming a 5 K temperature rise. The 17 readout modules of the Muon 
Rangestack will be cooled by fans mounted externally to their cases. All the modules 
will have internal temperature monitors and a hardware overtemperature interlock. 

The readout module will be constructed as a black anodized aluminum box with 
an oil-tight gasketted lid. Attention will be paid to minimizing enclosed or recessed 
spaces in the outer surface of the box, the screws and other hardware, and the con- 
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nectors and cables; furthermore the box will be purged with nitrogen internally after 
assembly. These measures should minimize the contamination of the scintillating oil 
with oxygen, which can degrade the light yield. 

Network, Data Acquisition System, and Clock and Trigger Distribution 

The 173 readout modules in the system buffer data internally, and communicate 
it, as needed, to the data acquisition computer system over a serial multidrop bus 
architecture. Although other implementations, in particular 10 Mbit/s Ethernet or 
CANbus, will be considered carefully at the design stage, at the time of this writing 
we assume a synchronous serial bus (separate data & clock lines) implemented with 
RS-422, running at 2 Mbit/s, using HDLC (High Level Data Link) protocol. For 
acceptable bus loading the readout modules will be arranged in five networks of up 
to 39 modules each (four used for the Vertex Detector and a fifth for the Muon 
Rangestack). The DAQ computer will be a Linux box with two DMA-capable PCI 
quad RS-422 communication controller boards installed. These are expected to be 
commercial off-the-shelf devices. Control and status communications will run over 
the same networks as the readout data. 

With 10% overhead for protocol the worst case module, with 1.02 kbytes of data, 
will require 5.6 ms for data transfer. Furthermore, assuming that the networks are 
physically arranged to balance the data load, each will see typically up to 250 hits per 
spill, which is 2.0 kbytes of data. Since this is coming from 39 different nodes, assume 
a 100% overhead for protocol. (In actuality the data is buffered in the readout mod- 
ules, so not every module will send a data packet on every spill, but for a reasonable 
worst-case estimate we ignore that here.) Then the complete data transfer on the 
network requires 20 ms (this includes the up to 5.6 ms for the worst case module). 
Since the spills are 66 ms apart, there is adequate spare capacity to accomplish this. 

The overall data rate into the DAQ computer, 16.6 khits/s = 133 kbytes/s, will 
be quite manageable, for transfer over the PCI bus, online monitoring software, and 
local disk storage or TCP/IP network transmission. 

The 32-bit coarse timestamp provided by the readout module rolls over in 429 
seconds. Therefore the DAQ computer will augment this timestamp with further bits 
to provide an absolute time measurement without rollover. 

A 10 MHz reference clock and global command word are distributed to the read- 
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out modules over separate, but similar, serial multidrop networks. The global com- 
mand word is used to synchronize the timestamp counter of all modules, and to fire 
the global trigger and test pulser events as described above. The serial communica- 
tions network clocks will also be derived from the same 10 MHz clock for convenience. 
The reference clock will be provided from an oscillator module which may either be 
a standard crystal oscillator or an absolute GPS-referenced timebase. 

Heritage 

Indiana University Cyclotron Facility personnel have significant recent experience de- 
signing and constructing a 9216 channel MAPMT readout system for the STAR End- 
cap Electromagnetic Calorimeter shower- maximum detector jHS]- Figure 13. 131 shows 
a picture of the R5900-00-M16 MAPMT together with the Cockcroft- Walton base 
and front-end electronics. The front-end electronics for STAR uses switched-reset 
gated integrators and 12-bit pipeline A/D converters to provide a 12-bit integrated 
amplitude measurement for every channel for every bunch crossing time (period of 
~ 103 ns) in the RHIC ring. Data is stored until an external trigger decision, and 
then all 9216 channels are read out in less than 20 fis. The readout is modular in 
groups of 192 channels (12 MAPMTs); the modules being enclosed in rugged, mag- 
netically shielded, water-cooled steel boxes which bolt to the STAR magnet poletip. 
While in functional details the system proposed here differs considerably from the 
STAR readout, many aspects of the design (such as the pipeline A/D converter and 
FPGA circuits, the conductively-cooled printed circuit board designs, the choice of 
connectors and other components, and test and diagnostic procedures) will either 
carry over directly or be informed by the lessons learned in the STAR readout work. 

4.1.4 Calibration 

The FINeSSE calibration system must calibrate the response of the MAPMTs of the 
fiber network to light from charged particles in the scintillator. It must also calibrate 
the energy response of the tracking detector. The relative gains of the muon range 
stack scintillator strips must also be understood. 

To calibrate the response of scintillating fibers and MAPMTs, FINeSSE will 
utilize cosmic ray muons and their Michel decay electrons. The energy distribution of 
Michel electrons is well defined and provides a "standard candle" for the calibration. 
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The veto fibers alone can be used to trigger on through-going muons. By using 
the veto and tracker sections, a cosmic Michel decay trigger can be set up. The 
cosmic muon events will put a known amount of light into the detector volume, with 
the geometry of the tracks constrained by the veto section, allowing accurate checks 
of light production by muons in the detector. The Michel events will provide enough 
low level light to calibrate the charge response of the fibers and photomultipliers. 
They will also provide a good sample for jp- reconstruction calibration. 



The muon range stack will also be calibrated using cosmic ray muons. The 
primary goal is to calibrate the relative gains and efficiencies of the scintillator strips. 
This should be easily achievable with the high rate of cosmic muons that FINeSSE 
will see. The calibration procedures using muons and Michels will be developed using 
Monte Carlo events. 



75 




Figure 4.14: Side view of the 25 ft below grade detector enclosure to house the FI- 
NeSSE detector. 

4.2 Detector Location and Enclosure 

The FINeSSE detector will be housed in a below ground enclosure 100 m along the 
line connecting the Booster neutrino production target and the MiniBooNE detector. 
Figures f4. 141 and 14. 151 show a side and top view of the detector within the enclosure. 
A detailed study, cost estimate, and schedule for this enclosure, performed by FESS, 
can be found in Appendix |B] 

The detector is installed using a gantry crane, above the enclosure, to move 
detector components from a truck, down through a removable hatch to the enclosure 
floor. The hatch can be rolled in and out of position on a daily basis to accommodate 
changes in weather and therefore in the installation schedule. Once installation is 
complete, shield blocks are installed above this to reduce cosmic ray rate in the 
detector. 
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Figure 4.15: Top view of the detector enclosure. 

The detector is accessible via a staircase with landings at a floor above the detec- 
tor and at the bottom of the enclosure. The floor above the detector holds electronics 
racks and computing needed for readout, triggering and monitoring. Ethernet from 
both detector sections is fed through the removable floor to these racks. 

Temperature control is minimal as per the requirements of the detector. Tem- 
perature of the Vertex Detector is controlled by a standard oil temperature control 
system like that used in the MiniBooNE experiment |64j . Temperature control of the 
electronics racks is done by enclosing them in a caboose and cooling or heating the 
caboose as necessary. 

This detector enclosure is minimal, but it is practical and fulfills our needs. We 
have studied the feasibility of a below-ground enclosure in detail and in conjunction 
with FESS. A similar enclosure constructed of sheet piling with a removable roof 
housed experiments at PC4. A senior FINeSSE scientist examined this PC4 enclosure 
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for toads, raccoons, and mice, and concluded this design will be sufficient. 

4.3 Detector Fabrication and Installation 

Detector components will be fabricated at Fermilab and FINeSSE collaborating uni- 
versities. Components will be assembled at a staging area before installation in the 
detector enclosure. Lab E or the New Muon Lab at Fermilab would be a suitable 
place for this staging area. Because the enclosure is below ground, the goal is to 
minimize the actual construction at the site. Subdetectors are constructed in large 
"packages" with outer protection. What follows is a description of fabrication and 
assembly for the Vertex and Rangestack components of the detector. 

4.3.1 Vertex Detector Assembly 

The detector assembly is constructed as a unit, at the FINeSSE staging area on site. 
Once assembled, it is moved to the FINeSSE enclosure and craned into position. 
Engineering of the transport is underway. Finally, the scintillator tank is filled with 
scintillator oil. The Vertex Detector assembly is described below. 

The main component of the detector is a cubic structure of six identical grid 
panels. Each of these is comprised of a stainless steel angle frame and an 80x80 
array of interlocking thin stainless steel ribs. These ribs are notched and positioned 
at an equal spacing of 3 cm, welded together at each intersection and to the frame 
around the perimeter. Each grid panel is constructed with the ribs centered with 
respect to the frame in one direction, and offset by 1.5 cm in the other. Mounting 
fixtures position the panels so that alignment may be made to assure that positions 
of the grid intersections are accurate and orthogonal throughout the cube. The offset 
directions of the panels are rotated so that any one direction of grid and the resulting 
fiber array is centered with respect to the other two directions. 

The fibers are held in place by riders on the intersections of the grid (see Figs. 14.161 
and 14.17]) . The riders are spring-loaded, providing a tension of 1 N per fiber. The de- 
formation by the combined pull of 6400 fibers of a 2.5m by 2.5m grid plane, supported 
on the edges, is about 1mm in the center. 

The fibers will be mounted in the completed support structure (after all six 
sides of the structure have been joined to form a cube). For mounting, an insertion 
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Figure 4.16: Supporting grid for the WLS fibers. Shown are the riders (slotted cubes) 
and the locking mechanism, which is pushed away from the rider by a spring. 




Figure 4.17: Same as Fig. \4-lb\ but seen from the inside of the fiber volume. 
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tool is needed that can hand a fiber end to the opposite side of the cube, reaching 
through a partly occupied grid cell. It is conceivable to do this with a group of fibers 
simultaneously. Once the fibers are available at both sides, the riders are put in place. 
Once set up, we estimate that the mounting takes on the order of a minute per fiber. 
This amounts to about 10 weeks (assuming 8 hours/day and 5 days/week) for the 
fiber assembly. 

The weight of the support structure is about 2 tons. This does not include the 
330 MAPMTs with associated electronics, which are mounted on the outside of the 
fiber support structure. When routing the fibers from the riders to the tubes one 
must take into account that the minimum recommended bending radius for 1.5 mm 
diameter fiber is 12 cm. The section of fiber outside the support lattice and the tubes 
are encased in light-tight Teflon sleeves. 

Once fiber stringing is complete, the electronics module support is put in place, 
the modules are mounted and PMTs plugged in, and each of the 300 'cookies' are 
locked to their PMT with spring loaded mountings. Cables suitable for submersion 
in scintillator oil are connected to the electronics modules. 

The completed fiber structure will be lowered into a tank which will eventually 
contain the liquid scintillator. This outer tank provides a cubic volume, 3.5 m on its 
side. The fibers exit the top of the tank. Equipment to detect light from the veto 
shield, i.e., the volume outside the fiber cube is mounted along the walls of the tank. 

The assembly is protectively enclosed and transported to the experimental loca- 
tion by a truck able to move 5 tons with a maximum tilt of 5° and acceleration of 
0.5 g. The onsite crane removes the detector from the truck and lowers it in place 
where fixturing assures its proper position. The hatch to the detector hall is then 
closed. The remaining veto area monitoring equipment is installed and cables from all 
electronics modules are routed to feedthrough panels. After all electrical connections 
are tested, the lid is put in place and the detector is filled with scintillator oil. 

The system for circulation and temperature control of scintillator is then put 
into operation, as is the system for maintaining the nitrogen atmosphere in the tank 
above the scintillator. The liquid is circulated by feeding from the top back to outlets 
distributed around the bottom of the tank. The circulation is necessary to control 
the temperature of the liquid, and to de-oxygenize the liquid scintillator initially 
(resulting in a gain of about 25% in light yield). The tank has a gas-tight lid, in 
order to maintain an oxygen-free atmosphere above the liquid. The cables from the 
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electronics readout modules in the liquid exit the tank through gas-tight ports above 
the liquid level. 

Temperature of the submerged PMTs and electronics is controlled by a system 
which circulates and cools the oil. Additional PMTs and electronics monitor a veto 
region surrounding the six sides of the detector active volume. Electronic cables are 
routed to feedthrough panels along the upper surfaces of the walls of the tank above 
the oil level, where a clean nitrogen atmosphere is maintained during the operation 
of the system. 

4.3.2 Muon Rangestack Assembly 

For the Muon Rangestack, we must consider fabrication of the steel absorber and 
the scintillator hodoscopes. Fabrication of the pieces will occur at a Fermilab facility 
away from the detector site. Then the modules will be transported to the detector 
hall. 

The steel will be procured by competitive bid from commercial machine shops 
and steel vendors. The absorber will have to be welded together from smaller readily 
available plate sizes, then cut to final size. No attempt has been made here to show the 
plate segmentation, as that will undoubtedly be determined by each steel vendor and 
the exact configuration of weld seams is irrelevant to the operation of the detector. 
Finished plates would be shipped to FNAL by truck. The route will have to be 
planned carefully, as will the fixturing on the truck, as these will be oversize loads. 
Calculation has shown that the plates can be rotated from a horizontal to a vertical 
position by pivoting them on the lower edge without exceeding the yield strength 
of the steel. However, once scintillator cans have been installed, a lifting/rotating 
fixture will probably be needed. Section 4 filler modules weigh 5.72 tons, and are 
the heaviest single pieces to lift. This puts a lower limit on the building crane in the 
assembly area at 6 tons, probably 10 tons when fixtures are included. 

There are 2112 scintillator bars to be extruded, cut to length and staged. These 
could be extruded at Fermilab jHSl- After cutting and finishing the ends, WLS fiber 
of the appropriate length will be laid in the groove and glued to the bar. The fiber 
gluing operation can be done after the bars are glued into the can bottoms. Large 
tables at least 4.6 m long will be used for curing the glue between the bars and fibers, 
and cans and bars. With staging areas, the floor space available should be at least 
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10 x 10 m 2 . Machines for mixing glue, applying it, and laying the fibers will be needed. 

The MINOS TDR showed the aluminum cans which will hold the scintillator 
as made on-site from coils of 3003 aluminum sheet 0.040 inches thick jH2j- A rolling 
machine was used to flatten the coils, and then to put corners along the edges used for 
crimping the lower and upper sheets together. We will have to examine the economics 
of this process, as having finished pieces delivered from a vendor may be more cost 
and time effective. The 4 m long pans are no problem for delivery by truck. 132 
bottom cans and 132 top plates will be made. 

The PVC spacers and fiber routing manifolds are very straightforward machined 
parts that can be fabricated at any machine shop equipped with CNC milling ma- 
chines. PVC is stable after machining and relatively cheap. 

Clear fiber bundles with their connectors must also be assembled, glued, and 
finished. Care must be taken to appropriately label fibers to assure that connections 
between scintillator cans and PMTs are correct. 

Assembly of the can proceeds as follows. The set of scintillator strips are glued to 
an aluminum skin to hold them fixed, and then placed within a can. The WLS fibers 
are routed through grooves milled into a PVC plate at the end of the can. The plate 
forces the fibers into their recommended 12 cm bend radius and prevents them from 
being damaged during handling. The plate is glued into the can and light-sealed with 
a cover plate of PVC plastic. An optical connector is also attached to the grooved 
PVC plate. This connector joins the WLS fiber to clear fibers that go to the PMT 
enclosures. The construction technique follows Chapter 5 of the MINOS TDR |62j . 
with the aluminum skins crimped at the edges to enclose the scintillator. 

An exploded view of the fiber routing components is shown in Figure 14" . 191 Fig- 
ure shows a rendering of a can without the final light-tighting layer of aluminum 
over the strips. The PVC cover is also not shown. The schematics do not show the 
clear fiber bundles, the PMT enclosures, and the clear fiber. 

Once scintillator cans are complete, their light tightness and other performance 
features may be tested. This is one of the big advantages of the 16-fold packaging 
of the scintillator. Cans are much easier units to handle for movement around the 
assembly and for testing than full planes. 

Finished cans will be attached to 0.50 inch protective steel backing plates at the 
assembly factory. Fixtures for vertical storage will be necessary. Finished scintillator 
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Figure 4.18: A can of scintillator strips with WLS fibers and enclosing PVC and 
aluminum sheet weighs just under 100 lbs and is approximately 4-3 m long including 
the PVC plates for protection of the fibers out to the optical coupling. 




Figure 4.19: An exploded view of the end of the can that routes fibers to the optical 
coupling. 
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Figure 4.20: A close-up of the fiber routing plate at the end of the scintillator can. 
Not shown are the cover plates that light-tight the assembly. 




Figure 4.21: A typical filler module for section 4- Not shown on any of the module 
pictures are lifting lugs to place the modules in the stack. The ears that allow the 
plates to hang from the support structure are welded on tabs. 
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Figure 4.22: Close-up of the support rail showing the "ears" on the steel plates that 
allow the plates to hang from the rails. The column and cross member at this corner 
were turned off to show more details. 

planes weighing 2.2 tons and being about 15 feet square will be trucked to the detector 
installation site. 

At the detector site, each module sits between two support rails upon which 
the Rangestack hangs, as indicated by Figure 14.221 The rails must be far enough 
apart to clear the ends of the scintillator can optical connectors. The protective steel 
backing, which is wider than the scintillator cans, is attached to each module. The 
steel backing is designed to protect the optical connectors during lowering into the 
stack. These plates are designed with scalloped edges to allow access to the clear fiber 
bundle, connector, and fastener exiting the side of the detector (y view; Figure I4.23j) 
once the module is in place. 

Because 4x4 m 2 plates are larger than commercially available plate stock, this 
design assumes that a 4 m square plate will be welded together from smaller plates, 
and water jet or plasma cut to final shape. Figure 14.221 shows a closeup of the ears 
on the steel plates that support the scintillator and filler modules. Figure 14.231 shows 
a closeup of the side of the detector with the y optical connectors and the scalloped 
steel edge. 

Once all scintillator planes and absorber fillers are completed, they will be trans- 
ported to the detector site. MI-12 or MI-8 will be used for temporary storage at 
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Figure 4.23: Close-up view of the changing thickness of the layers from section to 
section, and the scallops allowing access to the optical connectors on the y planes. 
This face will probably also have structure to support PMT boxes and catwalks for 
maintenance access. 

the site. For installation, each module will be made vertical (if they are not already 
on the truck) and lowered into the detector pit to hang on the two support rails. 
The scintillator planes will require a simple beam lifting fixture to prevent excessive 
deflection during crane handling operations. 

After installation of planes, the surrounding PMT boxes and support structure 
will be assembled. Clear fiber bundles will then be connected and testing with cosmic 
rays can begin. This can proceed with the building roof re-installed. 
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Chapter 5 

Event Simulation, Reconstruction, 
and Analysis 

The detector and physics events have been simulated and reconstructed to determine 
the feasibility, strategy, and sensitivity of the experiment for its two main physics 
goals: a measurement of As and disappearance. In this chapter the general fea- 
tures of the neutrino interactions and tracks in the detector are described. Next, the 
detector simulation and reconstruction programs and strategies are explained. Back- 
grounds, both beam-related and beam-unrelated are then discussed. In the final part 
of the chapter, the methods, tools, and estimated sensitivity are reported for the As 
and disappearance measurements. 

5.1 Interactions in the Detector 

The event rates expected in the FINeSSE detector are listed in Table 13.11 These 
events will need to be identified, counted, and measured in order to do the desired 
physics. The detector as described in Chapter H] will have the characteristics needed 
to do that. The most important and most prevalent events that the detector will 
see are listed in Table 15.11 along with a description of how they will "look" in the 
FINeSSE detector. 

The first three reactions are the most important for the main physics goals of FI- 
NeSSE. Examples of typical events of these reaction types as simulated with GEANT 
are shown in Figures I5.1H5.3I In these figures, charged hadrons are shown as solid 



event name reaction # tracks description 



CCQE 


un — > 


/i p 


2 /i,p: two-body kinematics 


NCp 


z/p — > 


z/p 


1 p: two-body kinematics 


NCn 


z/n — > 


z/n 


n: extraneous visible tracks from np 








scattering 


CCtt 


z/n — * 




>m 2 Not two-body kinematics 


NCtt 


z/p, n 


— * p, nXvr 111 ' 


>?» 1 Not two-body kinematics 



Table 5.1: Summary and description of event types that the FINeSSE detector will 
see. "# tracks" means typical number of charged particle tracks of significant energy. 



lines, muons as wide dashed lines, neutrinos and neutrons as dot-dashed lines, and 
photons as dotted lines. The detector is shown from the side. 

The individual particles within a particular event will be identified via their track 
length and energy loss density, u dE/dx" , as well as their decay patterns: 

• fj, ±m . Long tracks with low dE/dx. High-energy ^ will exit the Vertex Detector 
and enter the veto and perhaps the range stack. For /i ± that stop in an active 
area, the decay (Michel) electron will be observed. 

• p: Short tracks with high dE/dx. A 100 MeV proton travels approximately 
10cm in liquid scintillator. 

• n: Extraneous tracks from up scattering. Occasionally transfer enough energy 
in one collision so as to be misidentified as a p track. Will thermalize and 
capture in the detector yielding a delayed 2.2 MeV 7. 

• n^: Longish tracks that look like fi ± . For 7r that decay in the active area of 
the detector, the subsequent ^ and e ± can be observed. 

• 7T°: 2 hit clusters from the 7r°-decay 7 showers. 

• e^: 1 "fat" track from the e ± shower. 

Event Kinematics 

The two-body kinematics of NC elastic (uN — > uN) and CCQE (un —> /i~p) interac- 
tions are shown schematically with the use of the kinematic ellipses in Figure 15.41 
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Figure 5.1: A typical un — * fi~p reaction in the FINeSSE detector. The event vertex 
is at the origin of the indicated coordinate system. In this event the fi~ leaves the 
Vertex Detector, stops in the range stack, and decays. The two neutrinos from the 
muon decay are seen exiting the apparatus. The short recoil proton track is just visible 
at the origin. 




Figure 5.2: A typical up — > up reaction expected in the FINeSSE detector, generated 
at the origin. The short proton track can be seen, as can the final state neutrino which 
exits the apparatus. In this view, only the Vertex Detector is shown. 
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Figure 5.3: A typical un — > un reaction expected in the FINeSSE detector. Several 
interactions of the final state neutron can be seen. The final state neutrino exits the 
apparatus. In this view, only the Vertex Detector is shown. 

Due to the energy of the incident neutrino and the relatively low mass of the muon, 
the kinematics of the NCp and CCQE events are almost identical above E v « 300 
MeV. This simplifies the analysis and interpretation of the NC and CCQE event 
reconstruction. As can be seen in Figure 15.41 the final-state lepton may have any 
angle; the outgoing hadron has a maximum lab angle of 90°. The correlations can be 
seen in this figure as well. An event with a lower-energy high-angle lepton is paired 
with a low-angle high-energy hadron (and vice versa). 

Plotting energy vs. angle (Fig. 15.5)1 reveals strong correlations between the two 
variables for such two-body reactions. The events of most interest for the As analysis 
are low Q 2 events, where the proton has a low energy and a high angle. The lepton 
in these events will be in the forward direction at high energy. The events of most 
interest for the oscillation analysis are CC interactions with outgoing muons whose 
kinetic energy lies in the 0.1 to 1.5 GeV range. Many of these lower energy muons 
range out in the Vertex Detector and veto; the most energetic ones, at small angles, 
enter and range out in the Muon Rangestack. 

The effects of Fermi momentum of the bound nucleons can be seen in Fig. 15.51 
which compares the final-state protons produced in CCQE scattering from nucleons 
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Figure 5.4: Kinematic ellipses for the NC elastic (up — > vN) reaction. These are 
also valid for the CCQE (vn — > [i~p) reaction. The circles (ellipses) of increasing 
size indicate the CM (lab) momentum for the event at E v =500, 1000, 1500 MeV. 
The longitudinal component (parallel to the beam) of the particle momentum is plotted 
on the x-axis, and the transverse component on the y-axis. In a particular event, the 
particle momentum vector is constrained to lie on the appropriate ellipse. 
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bound in carbon (Fig. 15.5b ) to those from free nucleon scattering (Fig. 15.5b ). Fermi 
momentum widens the angular distribution of the outgoing proton and suppresses 
the number of nucleons at low momentum ("Pauli blocking"). The effect of this 
additional Fermi energy (« 25MeV) on the energy of the outgoing proton is small 
and will have minimal impact on measuring the Q 2 of the reaction (via Q 2 = 2m p T p ). 
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Figure 5.5: Kinetic energy vs. cosO for the (a) fi, (b) proton in CCQE scattering 
from bound nucleons, and (c) proton in CCQE scattering from free nucleons. Protons 
in NCp reactions will show the same distribution as that in (b). 



The correlations shown in Fig. 15.51 will be used to reduce backgrounds from NC 
and CC single pion reactions by requiring that the reconstructed tracks obey the 
illustrated kinematic constraints. 
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Figure 5.6: Diagram of the detector geometry as simulated by GEANT with a super- 
imposed CCQE scattering event. 

5.2 Simulation of the Detector 

The full FINeSSE detector, including the Vertex and the Muon Rangestack subde- 
tectors, has been simulated using the the GEANT [|J7| simulation package. A diagram 
of the apparatus as modeled by the program is shown in Fig. 15.61 

The Vertex Detector is simulated as an (2.4m) 3 of liquid scintillator with a 80 x 
80 grid of embedded 1.5 mm fibers with a 3.0 cm spacing in three orientations (XY, 
XZ, YZ; Fig. 14 .3J) for a total of 19,200 fibers. This inner volume containing the fibers 
is embedded in a larger volume of (3.5m) 3 of (p = 0.85 g/cm 3 ) liquid scintillator. The 
fiber support structure and tank walls are also included in the simulation. The area 



94 



outside of the fiber area in the scintillator tank (the "veto" ) is also assumed to be 
active. 

The Muon Rangestack, downstream of the Vertex Detector, is implemented with 
geometry as described in Section 14.1.21 which consists of alternating planes of scintil- 
lator and iron. 

In the inner region, scintillation light and Cerenkov radiation from passing charged 
particles is simulated. Photons thus produced are tracked until they impinge upon a 
WLS fiber or the edge of the detector volume, or are absorbed. In the outer 50 cm 
(veto) region of the liquid scintillator and in the Muon Rangestack, individual pho- 
tons are not tracked - only energy loss is recorded. However, this is not important in 
these regions as photostatistics are not crucial. 

In the active region of liquid scintillator, it is estimated that 5000 scintillation 
photons are produced in the liquid scintillator for every 1 MeV of energy deposited by 
a charged particle [53] • The absorption and capture efficiency of photons intersecting 
a WLS fiber is estimated to be 5% These fibers have typical attenuation lengths 
of 2.5 m, and the quantum efficiency of the PMT is approximately 20%. As a result, 
approximately 10% of the photons emitted in the capture cone of the fiber will make 
it to the PMT and produce a photoelectron. Combining these two factors yields 
an overall efficiency of 0.5% that an optical scintillation photon that strikes a WLS 
fiber will be detected at the PMT. To aid the speed of the simulation, the photon 
detection efficiency and production was combined so that 25 (5000 x 0.005) photons 
are produced per MeV. An effective efficiency of 100% for the photons that hit a 
1.5 mm WLS fiber was then assumed. In this way, the effects of photostatistics were 
properly simulated and the efficiency of the simulation was kept high. The attenuation 
length (5m) [S3] of the scintillator is fairly large compared to the size of the detector 
and is not a significant effect. The saturation due to large localized energy deposits 
("Birks' " Law) is also modeled. The production of Cerenkov photons is simulated 
but is negligible as the number of Cerenkov photons is only about 1% of that for 
scintillation. 

Using these factors, the simulation predicts that a proton track track passing near 
a fiber will create on average 10 photoelectrons in the PMT. This is consistent with 
the prototype test results reported in Chapter 0] (factoring in the difference in fiber 
length). Employing this method of tracking individual optical photons in the Vertex 
Detector assures that the photo statistics (with fluctuations) are properly simulated. 



95 



The simulation program can track single-particle events to study the detector 
response for each particle type. Alternatively, it can accept event descriptions as 
generated by the NUANCE MC program as described in Chapter El The later class 
are events that were used to predict sensitivities as described below. The output of 
the detector simulation includes: a list of all the "hit" fibers in the Vertex Detector, 
a list of all the bars in the Muon Rangestack that recorded energy loss, and the total 
amount of energy and time of deposit in each area of the detector. These data are 
subsequently passed through the event reconstruction program. 

Event Reconstruction 

The simulated event data from the Vertex Detector are analyzed with a reconstruction 
program employing a Hough transform technique [HH]. The Hough transform is a 
global track finding method that uses the hit fiber information from the XZ and YZ 
orientations. (The information from the XY orientation has not yet been used. This 
information can only improve the reconstruction.) The coordinates of each fiber (that 
recorded an amount of light over an adjustable threshold value) are used to calculate 

R = X(Y) sin a + Z cos a , (5.1) 

where a is a track angle and R is the perpendicular distance from the track to the 
origin. The track angle a is varied in a loop from —90° to +90°, and the R and a 
values for each hit used to make an entry (weighted by the amount of light in the hit) 
in a histogram. 

The task of track finding then reduces to finding "peaks" in this histogram. 
Finding single tracks is quite easy with this method. For events with multiple tracks, 
alternate methods had to be developed and tuned to subtract the light from the first 
track before the algorithm was employed to find subsequent tracks. At present, the 
reconstruction program is limited to finding a maximum of two tracks in each of the 
two 2D-orientations (XZ,YZ). 

The 2D-tracks were then combined to form 3D-tracks. The total energy and 
length of each track was also calculated, from which the dE/dx of the track could be 
determined. 

The simulated detector energy, angle, and position resolution for 50-500 MeV KE 
protons and muons are shown in Figure lo~71 These particle energies are representative 
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of the tracks that will be contained in the Vertex Detector for physics events. These 
proton kinetic energies correspond to Q 2 values of 0.1-1.0 GeV 2 in NCp and CCQE 
reactions. A Gaussian fit to the energy and angle resolution yields AE = 13(16) MeV 
and AO = 100(80) mrad for protons (muons). The quantity Av plotted in Figure 15771 
is the distance from the calculated track origin from the true origin. The simulations 
of single particles predict a mean Av = 9(10) MeV for protons (muons). The distri- 
bution is slightly wider for muons as they are longer tracks. These results indicate 
detector performance that will meet the physics goals of the experiment. The ef- 
fect of detector resolution on the physics distributions will be shown in the following 
sections. 



Comparison with IU Cyclotron Facility Prototype Tests 

The results of these resolution studies are consistent with the prototype tests per- 
formed at the IUCF (described in Chapter HJ), factoring in the differences in the two 
studies. For example, the prototype detector was much smaller (and so did not fully 
contain tracks); had only one fiber orientation; and used a slightly different fiber 
spacing. Nevertheless, the tracks reconstructed from the prototype test data yielded 
an angular resolution of approximately 6° (15° FWHM). The resolution quoted in 
Section I5~!?l for protons is also 6°. Similarly, the amount of light detected in the IUCF 
tests was 17 ±2 photoelectrons for a 200 MeV passing adjacent to a fiber. These sim- 
ulations, with light collection and detection efficiencies estimated from the individual 
detector properties, are in agreement with these results. 



5.2.1 Events in the Vertex Detector 

Several examples of the tracks obtained with the Hough transform reconstruction 
method from simulated data are shown in Figures [5.81 and [5.91 In these figures, the 
particle directions are indicated by the light-colored arrows. The reconstructed tracks 
and endpoints are indicated by dark lines and dots. As indicated, this method results 
in accurate reconstructed tracks for muons and protons down to kinetic energies of 
100 MeV. 
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Figure 5.7: TTie energy (a,b), angle (c,d), and position (e,f) resolution of the Vertex 
Detector as simulated and reconstructed for a sample of 1000 single particle events. 
The plots on the left (a,c,e) are for 50-500 MeV KE protons, those on the right (b,d,f) 
for 50-500 MeV KE muons. Only tracks that were fully contained in the Vertex 
Detector were selected. This effectively limited the upper muon KE to ~ 300 MeV. 
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Figure 5.8: A XZ (top) and YZ (bottom) projection view of a CCQE event in the 
simulated Vertex Detector with the reconstructed muon (long line) and proton (short 
line) tracks superimposed. In this event, T M = 820 MeV and T p = 150 MeV. 
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Figure 5.9: A XZ (top) and YZ (bottom) projection view of a NCp event in the 
simulated Vertex Detector with the reconstructed proton track superimposed. In this 
event, T p = 100 MeV. 
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Figure 5.10: The energy distribution of muons from CCQE interactions within the 
Vertex Detector fiducial volume, that enter and range out in the Muon Ranges- 
tack(white) overlaid on the distribution of all muons(red). 



5.2.2 Events in the Muon Rangestack 



Events in the Muon Rangestack considered here, have reconstructed vertices in the 
Vertex Detector fiducial volume, exit the Vertex Detector, and enter the Rangestack. 
These higher energy muons tend to be emitted at small angles as indicated by Fig- 
ure 15.5b . which shows the outgoing angle of the muon as a function of the muon 
energy. The energy distribution of these muons events that enter and stop in the 
Muon Rangestack is shown in figure 15.101 superimposed over the distribution for all 
muons. The distribution of hits as a function of transverse position in the range stack 
created by these higher energy muons is shown in Figure 15.111 Note that there are 
few hits near edges of the Rangestack. A muon is considered to have ranged out in 
this subdetector if its last hit is not at an edge of the detector, either in the last 
plane, or one of the edge strips of an earlier plane for higher angle muons. Most 
muons below 1.5 GeV range out within the Rangestack as indicated by hits deposited 
in each rangeout layer (Figure 15.12)) . 
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Figure 5.11: Number of hits in scintillator strips in the Muon Rangestack from CCQE 
muons that originate in the Vertex Detector fiducial volume. The "Bar Number" 
measures the hit coordinate transverse to the beam direction. 
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Figure 5.12: Number of hits in each scintillator layer in the Muon Rangestack from 
CCQE muons that originate in the Vertex Detector fiducial volume. The "Layer 
Number" measures the coordinate along the beam direction. 
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5.3 Backgrounds 

We recognize two major classifications of backgrounds for FINeSSE: beam related 
and beam unrelated. 

FINeSSE can expect to see background events from neutrino interactions in the 
detector through channels other than CC and NC events, as described in Sec- 
tion 13.21 Some of these will be misidentified as signal channel interactions. There 
will also be neutrino interactions in the dirt around the detector, sending charged 
particles and neutrons into the detector. These are the most important beam related 
backgrounds. 

FINeSSE will also see a large flux of cosmic ray muons and neutrons. These 
events can be used for detector calibration, as described in Section 14.1.41 but they can 
also be misidentified as beam events and therefore represent a potential background. 
However, the small probability of an interaction in time with the short beam spill of 
the booster neutrino source (1.6 /xs) reduces these backgrounds greatly. In addition, 
they can be accurately measured and subtracted. 



5.3.1 Beam Related Backgrounds 

The NC and CC measurements will have very different backgrounds, due to differences 
in their event topologies. The CC events have an outgoing muon and a recoil proton. 
These two particles emitted from the interaction vertex will usually provide a very 
clean signal. However, at low Q 2 , where As is measured, the recoil proton will have 
very little energy. In these cases, charged particles from other interactions with two 
final states could mimic this short, proton-like track. The NC events have only the 
recoil proton. Interactions producing only low energy short tracks can mimic this 
recoil proton. These will be the hardest events to identify. 

Tables l3~Tl and l5~3l list the expected event rates from the neutrino beam. FINeSSE 
expects to see 24,435 CCQE events and 9864 NC elastic events per 1 x 10 20 POT, 
assuming 100% detection and reconstruction efficiencies. Using the present selection 
cuts, the expected event totals are 2990 CCQE events and 1604 NC elastic signal 
events per 1 x 10 20 POT. Because the As analysis will use a ratio to eliminate many 
systematic uncertainties, background events become a significant challenge to the 
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measurement. Table I5~3l summarizes the contributions to the signal from beam related 
backgrounds that originate in the detector. 

Neutrons from beam neutrino interactions in the material around the detector 
are an additional potential background. These events arise from neutrino NC and 
CC interactions in the material around the detector that create secondary particles 
that impinge on and interact in the detector. The majority of the events cannot be 
eliminated by taking advantage of the Booster neutrino beamline timing structure, as 
they will be in time with the beam. The majority of the particles that mimic a signal 
are neutrons that are produced in the surrounding dirt. They have a fairly large 
range and can travel into the central area of the detector without leaving a signal in 
the edge. Muons, can not do this and, therefore, are not a background. 

Initial Monte Carlo studies, employing the cuts described in Section l5.3.41 indicate 
that the number of neutrons that interact in the dirt around FINeSSE and pass the 
particle ID cuts are insignificant compared to the signal. 

An additional background to consider are neutrons directly produced by protons 
in the booster neutrino source beam dump. These neutrons will constitute a negligible 
background at the FINeSSE detector enclosure 50m downstream from the dump. 
There is no significant neutron flux from the beam dump after 25m downstream due 
to absorption in the earth. This was determined in shielding assessment studies for 
the Booster neutrino source [p^] . 

5.3.2 Beam Unrelated Backgrounds 

The main beam unrelated backgrounds stem from cosmic rays. After selection cuts, 
these are negligible in comparison with the NC elastic signal. The cosmic ray muon 
rate in the FINeSSE detector fiducial volume, beneath 1 m of concrete and 0.5 m of 
veto region scintillator, will be about 900 Hz. The muons are a potential background 
both for CC and for NC events. The relative rate of coincidences with the beam timing 
window is 0.0014. This can be reduced to 1.4 x 10~ 6 by factoring in a conservative 
veto efficiency of 0.999 for muons. At 0.5 m thick, the veto for the Vertex Detector 
will be more efficient than the 0.3 m thick MiniBooNE veto, which has an efficiency 
better than 1 in 10 4 for through-going muons. 

Unlike the beam related backgrounds, the number of background events con- 
sidered here depends on the Booster performance. If we assume that the Booster 
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delivers 4 x 10 12 protons per pulse, then we expect to see 36 non-vetoed cosmic ray 
muon events per 1 x 10 20 POT in coincidence with the beam. Event analysis will 
reduce this background down further; since there will be no recoil proton and the 
muons will all be headed downward. Therefore, cosmic ray muons will present a 
negligible background to CCQE or NC events. 

Cosmic ray neutrons have also been considered as a background. Their average 
rate above 50 MeV during periods of normal solar activity at sea level and ~ 40° 
geomagnetic latitude is approximately 9 x 10~ 3 sec _1 cm~ 2 [701 They have a 

momentum spectrum that falls very steeply with energy. Fewer than 5% of these can 
traverse the concrete shielding and the veto region. We therefore expect a cosmic 
neutron rate of < 1Hz in the FINeSSE Vertex Detector, for neutrons above 100 
MeV. This rate will be reduced to a negligible level because of the small coincidence 
probability with the short beam spill. 

5.3.3 Charged and Neutral Current Event Identification 

Neutral-current elastic scattering events (up — > up) will be identified in the FINeSSE 
detector by looking for single proton tracks consistent with elastic scattering kinemat- 
ics. A track is identified as a proton by a large dE/dx. Charged-current quasi-elastic 
scattering events (un — > yTp) will be identified by looking for events with two tracks 
each consistent with the expected dE/dx. In addition, other cuts are employed to 
reject backgrounds. The strategy will be to maintain a compromise between large ef- 
ficiency for low-Q 2 events while keeping backgrounds as low as possible. The squared 
four-momentum transfer, Q 2 , will be determined event by event, by measuring the 
energy of the proton in both NC and CC events. Q 2 is determined from the energy 
via Q 2 = 2m p T p . 

5.3.4 Simulation Results 

The GEANT simulation of the FINeSSE detector and reconstruction program as de- 
scribed above is used to estimate efficiencies and backgrounds for the NC/CC ratio 
measurement. It is also relied upon to determine the experimental error on a mea- 
surement of G S A . 
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Cuts 

The experimental cuts employed for up — > up (NCp) and un — > pTp (CCQE) event 
identification are summarized in Table 15.21 A reasonable compromise of efficiency 
and background rejection could be obtained out to within 15 cm of the edge of the 
active Vertex Detector volume. This results in a fiducial volume of (2.2 m) 3 = 10.6m 3 
or 68.4% of the total (2.4 m) 3 "active" volume of the Vertex Detector. With p = 
0.85 g/cm 3 liquid scintillator, a fiducial mass of 9.1 t is obtained. 



cut # 


NCp cuts 


CCQE cuts 





edge distance < 15cm 


edge distance < 15cm 


1 


# 3d tracks = 1 


# 3d tracks. eq. 2 


2 


dE/dx(p) > 2.5 


dE/dx(p) > 2.5, dE/dx(p) < 2.5 


3 


0(p) > 0.5 


6(p) + 6(n) > 1.5 


i 


no "late" light in vertex det. 


no "late" light in vertex det. 


5 


no veto or muon stack energy 


low "remaining" energy 



Table 5.2: Cuts used to identify "NCp" and "CCQE" events in the simulated sample. 

Cuts 1 and 2 identify and separate NCp and CCQE based on the number and 
type of charged particle tracks in the event. Cut 2 accepts events with tracks of larger 
angle, a feature of the two-body kinematics of these events, and rejects background 
events which tend to have tracks at smaller angles. Cut 4 lowers background by 
eliminating pions that stop and decay in the Vertex Detector into muons (which 
decay "late"). This also has a slight effect on the efficiency for accepting CCQE 
events, but it is not large since most of the muons in CCQE scattering leave the 
Vertex Detector. Cut 5 lowers backgrounds for both channels. For NCp events, it 
further lowers the pion (charged and neutral) backgrounds by cutting those events 
where the pions leave little signal in the Vertex Detector. For CCQE, events cut 5 also 
reduced pion backgrounds as a cut on low "remaining" energy (energy not assigned 
to tracks after tracking is complete) reduces the number of events with extra energy 
(due to pions). 

Efficiencies and Purities 

A sample of 215k NUANCE-generated events were tracked through the detector with 
the GEANT simulation with vertices evenly distributed within the nominal volume 
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of the Vertex Detector (2.4 m) 3 . These events are then reconstructed using the 
algorithms described above. The results are summarized in Table I5~31 



reaction channel 





IN up 


IN V^yll 








raw events 


91910 




1 Q0R9 

ryuoz 


1 nm 09 

1UU1UZ 


£\A 1 07 


passed events 


3929 


1162 


167 


48 


4 


efficiency (%) 


18.5 


5.7 


0.9 


0.0 


0.0 


fid. eff. (%) 


27.1 


8.3 


1.3 


0.1 


0.0 


purity (%) 


74.0 


21.9 


3.1 


0.9 


0.1 


CCQE cuts 


NCp 


NCn 


NCtt 


CCQE 


CCtt 


raw events 


21219 


20487 


19062 


100102 


54107 


passed events 


165 


76 


581 


7323 


1322 


efficiency (%) 


0.8 


0.4 


3.0 


7.3 


2.4 


fid. eff. (%) 


1.1 


0.5 


4.5 


10.6 


3.6 


purity (%) 


1.7 


0.8 


6.1 


77.4 


14.0 



Table 5.3: Summary of events that passed the NCp and CCQE cuts along with ef- 
ficiencies and purities: "efficiency" is the reconstruction efficiency throughout the 
(2.4 rn) 3 volume, "fid. eff." is the reconstruction efficiency within the (2.2 m) 3 fidu- 
cial volume. The simulation data set contained 215k events. To scale to the total 
expected FINeSSE event count, multiply these numbers by 2.45. 

The efficiencies and Q 2 resolutions for the NCp and CCQE "signal" samples are 
shown in Figures 15.131 and 15.141 Note the large efficiency for reconstructing NCp, 
and reasonable efficiency for CCQE, in the Q 2 = 0.2 — 0.4 GeV 2 region. At low Q 2 , 
proton tracks from NCp are short and difficult to reconstruct, causing a fall off in the 
efficiency. At high Q 2 , final state particles tend not to be contained for both the NCp 
and CCQE samples, causing a fall of in the efficiency in these regions also. In addition 
to this, very forward events in the NCp sample are cut to remove backgrounds from 
pions. This further decreases the NCp efficiency at high Q 2 . 

5.3.5 Method to Extract As from the Data 

To extract As (or G S A ) from the data, the NC to CC ratio must be formed and 
compared to predictions obtained from Eq. 12.51 This ratio must be corrected for 
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Figure 5.13: (a) Efficiency for reconstructing true NCp events with the NCp cuts 
within the fiducial volume as a function of generated Q 2 and (b) the distribution of 
the difference between generated and reconstructed Q 2 for this sample. The a of this 
distribution is 0.053 GeV 2 
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Figure 5.14: (a) Efficiency for reconstructing true CCQE events with the CCQE cuts 
within the fiducial volume as a function of generated Q 2 and (b) the distribution of 
the difference between generated and reconstructed Q 2 for this sample. The o of this 
distribution is 0.079 GeV 2 
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backgrounds and detector efficiency and may be written as: 

R(NC/CC) = £gc (^Nc-Eie N c,^Nc,t) ) (5 2) 

£nc (A"cc — l^j e cc,jNcc,j) 

where e^c ( e cc) is the NC (CC) reconstruction efficiency, and N NC (Nqc) is the 
number of NC (CC) identified events in the detector. The quantities eNc,i> e cc,j, ^Ync.i, 
and Ncc,i are the background efficiencies and calculated number of background events 
for the NC and CC channels, respectively. 

Each of the terms in Equation l5.2l will contribute to the total error on R(NC /CC). 
The relationship between the error on G S A (cr(G A )) and R is quantified in Eq. 12.181 
so we can estimate the error on G S A due to each of the terms in Eq. 15.21 These errors 
are discussed and quantified in the sections below, assuming that the measurement is 
performed in the Q 2 = 0.2 — 0.4 GeV 2 bin, and that an estimate of cr(G A ) is obtained. 
After this treatment, a fitting procedure that uses the entire range of Q 2 is used, to 
illustrate an improved method for extracting G S A . 



Statistical Errors 

The event sample of 215k NUANCE-generated events is 0.408 of the total number 
of neutrino-scattering events expected in the (2.4 m) 3 active volume of the Vertex 
Detector with 6 x 10 20 POT. So the number of events in Table IS~3*1 mav be multiplied 
by 1/0.408 = 2.45 to obtain the total number expected. Taking these scaled numbers, 
and considering the efficiencies for the NCp and CCQE cuts, yields 7341 (9936) 
accepted events in the NCp (CCQE) sample in the Q 2 = 0.2—0.4 GeV 2 bin. Statistical 
errors from these two sources will contribute 1.2% and 1.0% to the relative error on 
R(NC/CC). 



Neutral Current Neutrino Neutron Scattering 

One of the backgrounds in the numerator of Equation 15.21 is that due to the NC 
process vn — > vn (labeled "NCn" in Table I5~3*|) . This reaction will need to be treated 
separately, since it is the largest background to the NCp channel, and also since the 
rate depends upon G S A . 
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The dependence of cross section for un — > un on G S A is different than that for 
up — > up, due to the isospin factor r z (cf. Equation I2.9j) . This will "dilute" the ex- 
perimental sensitivity. For this reason it is important to keep the efficiency of the 
un — > un channel as low as possible. This is a difficult task as neutrons in a liquid 
scintillator detector often undergo charge-exchange reactions which create energetic 
protons. These reactions can be indistinguishable from the up — > up reaction. Neu- 
trons, though, often create two or more protons of lower energy and these topologies 
will not be misidentified as up — > up. 

As can be seen from Table 15.31 the reconstruction efficiency for the NCn back- 
ground is 8.3% and comprises 21.9% of the NCp sample. It is non- negligible. It would 
be desirable to improve the detector and reconstruction to lower this background. It 
is likely that this is possible by observing the delayed light due to neutron capture. 
The detector as currently designed would have the ability to observe this delayed 
light, and this would improve the neutron background rejection considerably. This 
technique, however, has yet to be investigated fully, and any potential gains are not 
assumed in this analysis. 

The effect of NCn events in the NCp sample dilutes the sensitivity of the measured 
NC to CC ratio to G A by a factor (1 — 2eNCn/ e NCp)- The impact of this on the 
analysis can be determined using the ratio of efficiencies, eNCn/eNCp — 0.19, in the 
Q 2 = 0.2-0.4 GeV 2 bin. Applying this efficiency factor, the sensitivity of R(NC/CC) 
is reduced from 1.2 to 0.75. This is a fairly substantial reduction but not large enough 
to prevent a measurement. The uncertainty in G S A becomes: 

a{G s A ) = 1.3AR/R. (5.3) 

A 5% relative measurement of the NC/CC ratio at Q 2 = 0.25 GeV 2 would enable an 
extraction of G S A with an error of 0.07. 

Scattering from Free Protons 

The FINeSSE neutrino scattering target is mineral-oil-based liquid scintillator (mainly 
CH2). The hydrogen will provide two additional free protons for up — > up scattering. 
This entails a correction to the measured NC to CC ratio. The measured ratio (after 



Ill 



background subtraction) will be: 

N cc N cc { ' J 

where N$ c and N^c are the number of events from up — > up on bound (in carbon) 
and free protons. This can be written: 

NnC - N " c x(1 + T), (5.5) 



N CC N cc 
where: 

/yfree 
JV NC 

Naively, one would expect T to be equal to the free proton to bound ratio in the 
target. This is ~ 0.34 [72] for mineral-oil-based liquid scintillator. In actuality, JF 
will be slightly larger than these numbers, as the the cross section per nucleon is 
slightly suppressed (by « 10%) for bound nucleons. In addition, the value used for 
T will depend upon the Q 2 region accepted for the analysis. So for the purposes of a 
conservative estimate here, T = 0.34 * 1.10 = 0.38 will be used. 

The desired corrected ratio (that on pure carbon) can then be written: 

R{NC/CC) = = ; and (5.7) 

iv CC I 1 + •> ) iv CC 

R 1 + F T T K J 

The number of bound and free protons in the target will be known, so the relative 
uncertainty on T will be dominated by the uncertainty in the ratio of the bound to 
free NC cross sections. An estimate of this relative uncertainty is 5%. This will add 
a systematic uncertainty in the NC to CC ratio of ~ 1.4%. 



Systematic Error in Efficiencies 



The flux and target density is same for both NC and CC reactions and does not enter 
as an additional source of uncertainty in a measurement of the ratio, R(NC/CC). 
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However the reconstruction efficiencies will not be the same for both reactions, and 
the systematic error in the ^ term of Eq. 15.21 must be considered. 

At a given Q 2 , the proton in the NC and CC reactions has the same kinematics. 
Because of this, some of the systematic error is highly correlated and the dominant 
uncertainty is in the efficiency to detect the muon in the CC reaction. This is a 
well- understood process and it should be possible to determine it with high accuracy. 
For these reasons, a relative error of 3% is estimated for ^ and will contribute at 
this level to the relative error on R(NC/CC). 

Systematic Error in Backgrounds 

The background subtraction in both the numerator and denominator of Equation 15. 21 
also contributes systematic error to R(NC/CC). This error has a contribution from 
both the absolute rate of the background and the efficiency to be reconstructed as 
a NCp or CCQE event. The relative error on R(NC/CC) due to these uncertainties 
is diluted to the reasonably small background subtraction. This dilution factor is 
equal to the ratio of background to (true) signal for NCp and CCQE events. From 
the numbers in Table l5~3l this dilution factor is 0.35 for NC events and 0.29 for CC 
events. 

It is estimated that the relative error on the product of efficiency times back- 
ground can be kept at or below 10% for both NC and CC reactions. This will be 
possible by means of detailed simulations and studies of the detector, and having 
developed a good model for the background processes. These backgrounds have been 
measured in other neutrino experiments, and likewise, will be measured in FINeSSE. 
In addition, relevant electron scattering data exist. 

The relative error of 10% on the background contributions contributes 3.5% 
(2.9%) via NC (CC) backgrounds to the relative error on R(NC/CC). 

Systematic Error Due to Q 2 Reconstruction 

A large error in the Q 2 reconstruction of events will distort the overall distribution and 
add to the error in R(NC/CC). For example, if low-Q 2 (< 0.2 GeV 2 ) CCQE events 
(of which there are many; see Figure I2.3|) are reconstructed in the Q 2 = 0.2 — 0.4 
GeV 2 bin, a larger error will result in R(NC/CC) that is not quantified in the errors 
considered above. 
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From examination of Figures \5. 131 and 15. 141 it can be seen that the Q 2 resolution 
of this detector is quite good for both NCp and NCn events. Also, there is no evidence 
of a large skew in these distributions. The Q 2 resolution is smaller than the bin width 
and will contribute a negligible error to R(NC/CC). 

This hypothesis is further checked in the fitting procedure described below as the 
reconstructed Q 2 is used to bin the data. 

Nuclear Model Uncertainties 

The majority of protons and neutrons participating in up — > up and un — > fi~p 
interactions are not free, but are bound within a carbon nucleus. Any nuclear effects 
that are not symmetric between protons and neutrons can lead to a misinterpretation 
of the data, when analyzed using a free nucleon model. For this reason, a model with 
a realistic treatment of nuclear effects will be employed for the final analysis. 

The sensitivity of the NC and CC cross sections to the specific nuclear model 
employed has been investigated by several different groups EU IH2] • The effects 
on the absolute cross sections can be quite dramatic. However, the distortion in the 
ratio R(NC/CC) has been shown to be small. It has been estimated to introduce an 
error on 0.005 in G S A 

This will be an area of continuing work. Results up to now, however, indicate 
that R(NC/CC) is insensitive to the details of the nuclear model employed for the 
calculations. 

Form Factor Uncertainties 

The (non-strange) form factors in Eq. 12.51 have been measured in a variety of exper- 
iments at many different Q 2 values. This is not true for the strange form factors. 
Information on the strange-axial form factor G S A is almost completely non-existent. 
Accumulating more information on this form factor is one of the main goals of this ex- 
periment. The strange vector form factors influence the NC to CC ratio only weakly 
and they will be precisely measured by the upcoming GO experiment ^H] to take 
place a Jefferson Lab in the near future. 

In addition, the contribution from these unknown form factors is minimized by 
performing these measurements at low-Q 2 . The sensitivity of the final result to these 
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other form factors is small as will be shown in the results from the fitting procedure 
below. 

Summary of Errors on R(NC/CC) 

The experimental systematic errors discussed in the proceeding sections are summa- 
rized in Table I5~4l These contributions, added in quadrature, yields a relative error 
on the ratio, R(NC/CC), of 5.8%. This, combined with Equation 15.31 yields a total 
error on G S A of 0.075. 



quantity, 




relative error, 


contribution to 


Q 


prefactor 


AQ/Q 


AR/R 




1 


1.2% 


1.2% 


N C c 


1 


1.0% 


1.0% 


JV NC 


0.27 


5% 


1.4% 




1 


3% 


3.0% 


Si e NC,i-^NC,i 


0.35 


10% 


3.5% 


Si e cc,iNcc,i 


0.29 


10% 


2.9% 


Total experimental 


error 


5.8% 



Table 5.4: Estimated experimental errors on a measurement of R(NC/CC) at Q 2 = 
0.2 — 0.4 GeV 2 with FINeSSE. The prefactor relates the relative error on the quantity 
to the relative error on R(NC/CC). The details are explained in the text. 

This is a conservative estimate. This error will be reduced considerably by re- 
ducing the vn — > vn background (via detection of photons from neutron capture), 
which both contributes to the error in R(NC/CC) and lowers the sensitivity of G S A to 
R(NC/CC). 

5.3.6 A Fit Procedure to Extract As from R(NC/CC) 

In the error analysis explained above, only one Q 2 bin was considered. This yielded 
a conservative estimate for the error on G S A . As can be seen in Figures \h . 1 31 and 15 . 141 
high reconstruction efficiency is obtained from Q 2 = 0.2 to 0.8 GeV 2 , and a mea- 
surement of R(NC/CC) may be made in this range. This larger Q 2 range will accept 
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more events and will allow for the systematic study of the behavior of R(NC/CC) 
and, consequently, the form factors as a function of Q 2 . 

To utilize the entire data set, a fit procedure was developed to extract G S A , to 
determine the effects of systematic errors, and to investigate the dependence of the 
extracted G S A on the other form factors. In this procedure, a corrected NC to CC ratio 
is formed using Equation 15.21 for each Q 2 bin of width 0.1 GeV 2 using the simulated 
data. The backgrounds are subtracted and efficiency corrections are applied. The 
predicted ratio is then calculated using the free nucleon cross section (Eq. I2.5J1 . cor- 
recting for the NC scattering from free protons and the estimated contribution from 
vn — > vn. This prediction is calculated for a range of G S A values. The dependence on 
the other form factors (F(, F 2 2 ) and the axial mass M A is also studied. 

The sensitivity of the data to G S A is determined by calculating: 



where the sum runs over Q 2 bins. The la error on G S A is determined by the G S A value 
where the \ 2 value changes by 1. 

This procedure yields an error on G S A of 0.039. This is considerably smaller than 
the error quoted above due to the fact that this fit procedure uses more of the data 
through a range in Q 2 . 

The sensitivity of R(NC/CC) to other form factors is also studied with this 
method. Figure 15.151 shows how a measurement of R(NC/CC) depends on other 
form factors. An uncertainty in the axial vector mass, M A , which parameterizes the 
Q 2 evolution of the G\ form factor, effectively adds a small additional uncertainty to 
the extracted value of G S A . The current world average of Ma as measured in neutrino 
scattering is 1.026 ± 0.021 GeV [2Z1- This uncertainty will add approximately 0.015 
to the uncertainty in G S A . The sensitivity of G S A on F% is weaker, but this form factor 
is not well known directly. Linear combinations of this form factor with F( have been 
measured [12]. A conservative estimate on the uncertainty on F| is ±0.1. An uncer- 
tainty of this size in F| will contribute an uncertainty of ±0.02 to G S A . Combining 
these yields a ±0.025 uncertainty on G S A . 

Based on the studies described above, FINeSSE will be capable of a measurement 
of G S A with small errors. An uncertainty of approximately 0.04 results from consider- 
ing statistical and systematic experimental errors. Uncertainties in the form factors 




(5.9) 
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Figure 5.15: \ 2 contours corresponding to ler (solid), la (dashed), and 3a (dotted) 
resulting from the fit procedure described in the text. In (a) G S A and Ma, the axial 
vector mass, are varied, in (b ) G S A and F 2 S are varied. 
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Events for 6 x 10 20 POT 
Location 


25 m Absorber 


50 m Absorber 


100 m (9 t) 
200 m (40 t) 
550 m (650 t) 


103,000 
97,000 
205,000 


195,000 
162,000 
323,000 



Table 5.5: Estimated event sample for 5 x 10 POT with detectors located at vari- 
ous distances with fiducial volumes as given. Rates are given for a 25 m and 50 m 
absorber. 

adds an additional contribution of 0.025. This level of accuracy is unprecedented in 
v scattering will allow for stringent tests with theory. 

5.4 Details of Disappearance Sensitivity Studies 

The method for determining the disappearance sensitivity compares the energy 
distribution of events in the near (FINeSSE) and far (MiniBooNE) detector. In 
this section we justify the choice of the experimental setup and the estimates of 
the systematic errors. We studied various configurations for FINeSSE+MiniBooNE 
running in order to identify the optimal experimental set-up. Here, we reduce the 
comparisons to two representative positions of the near detector (100 m and 200 m) 
and the two possible lengths of the decay region (25 m and 50 m). The event rates 
assumed in these studies are indicated in Table 1531 

5.4.1 Tools and Assumptions for This Study 

Determination of the disappearance sensitivity requires tools from both FINeSSE 
and MiniBooNE. Assumptions related to the flux in both detector have been described 
in Chapter El As discussed there, the NUANCE Monte Carlo was used to simulate 
events in both detectors [53]. The following standard MiniBooNE detector cuts were 
used: tank hits > 100, veto hits < 6, and reconstructed fiducial radius < 500 cm.[3|. 

FINeSSE and MiniBooNE are assumed to have the same reconstruction efficiency. 
This is conservative because the finely segmented FINeSSE detector will have much 
better resolution on the reconstructed neutrino energy in quasi-elastic (QE) events. 
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2.5 3 
Energy(GeV) 

Figure 5.16: True kinetic energy distribution for all muons from CCQE interactions 
in the Vertex Detector fiducial volume shown in red. Overlaid in white is the sample 
of those with a contained muon. Also shown are contributions corresponding to muons 
stopping in the Vertex Detector (green), the Vertex Detector veto (blue), and the Muon 
Rangestack (purple). 



Energy resolution for muons ranging out in the Muon Rangestack will also be known 
to 10% (see Section r4.1.2jl . FINeSSE will also have much higher purity for the QE 
sample due to its excellent low energy reconstruction. For events to be counted in the 
FINeSSE data sample for this analysis, outgoing muons from CC interactions must 
be contained entirely in the Vertex Detector, in the Vertex Detector plus the veto, or in 
these plus the Muon Rangestack. Contributions from these three contained samples, 
as a function of energy, are shown in Figure 15.161 Wiggles in this distribution, seen 
clearly in Figure 15.171 are due to contributions from one distribution turning off as 
another turns on. Folding this acceptance into overall number of interactions expected 
yield the event rates summarized in Table 15*31 

5.4.2 The Motivation for 25 m Absorber Running: Parallax 

The main physics driver for the FINeSSE+MiniBooNE standard configuration was 
reduction of "parallax" in the near detector. In this section, we describe the cause 
of parallax, why it is a significant issue for disappearance experiments, and how we 
have mitigated the effect with our design. 
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Figure 5.17: Acceptance versus true kinetic energy for muons from the CCQE event 
sample. 
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Figure 5.18: Left, solid: Angular acceptance of the FINeSSE detector with the 50 m 
absorber in radians; Right, solid: Acceptance with the 25 m absorber in radians. 
Dashed (both plots): acceptance of the MiniBooNE detector in radians. 
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Figure 5.19: Left: Ratio of the flux in the FINeSSE to the MiniBooNE detector 
using the 50 m absorber; right: Ratio when 25m absorber is installed 



In an ideal world, there are two reasons why disappearance experiments based on 
near/far detector comparisons would be performed with point-source beams. First, 
the L of each detector is well determined. Second, the near and far detector have 
the same energy and angular acceptance for the neutrinos. Unfortunately, the only 
true point sources come from beams created by stopped pion and muon decay (for 
example at LANL or SNS). These neutrinos are so low in energy that a muon cannot 
be produced in CC interactions. Hence v ^ disappearance searches are impossible. 

To produce a higher energy neutrino beam, an extended decay path is necessary. 
In this case, the modification of L for the decay length must be taken into account. 
Also, the angular acceptance of the detector will vary depending on the point along 
the line of the decay. For a detector which is relatively far from the beamline, such as 
is the case with MiniBooNE, these corrections are negligible. But for a near detector, 
this "parallax" can be a large effect. 

Because there is a strong correlation between energy and angle of the neutrinos, 
differences in angular acceptance between the near and far detectors is a major con- 
cern. These differences translate directly into differences in the energy distribution of 
events. When taking a near/far ratio, the result are "wiggles," which are not related 
to oscillations. As long as those wiggles do not mimic an oscillation signal, this is not a 
major problem. But with realistic detector efficiencies and smearing, these variations 
can severely limit sensitivity. Therefore, a top priority of any near/far experiment 
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Figure 5.20: Comparison of the parameter space covered by FINeSSE+MiniBooNE 
with the 25 m absorber (red) and the 50 m absorber (black). The unwanted "wiggles" 
in the 50 m absorber allowed region are caused by the variations in the ratio of fluxes 
at FINeSSE and MiniBooNE. 




Figure 5.21: Comparison of the parameter space covered by FINeSSE+MiniBooNE 
for a detector at 100 m (black) from the target versus 200 m (red) from the target. 
Sensitivity to high Am 2 is lost at 200 m because neutrinos will oscillate before reaching 
the FINeSSE detector. 
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must be to mitigate parallax. 

An obvious solution is to move the near detector sufficiently far from the line 
that the effect of parallax dies away. However, for a short-baseline disappearance 
experiment, this reduces the reach in Am 2 because some neutrinos have oscillated 
prior to reaching the near detector, thereby reducing the apparent flux. 

For the FINeSSE+MiniBooNE running, we propose a configuration in which the 
25 m absorber is installed in the beamline. Figure I3~T1 shows a schematic of the beam 
design, including the position of the enclosure for the 25 m absorber. This absorber 
is presently hanging above the beamline. It can be lowered (and subsequently raised) 
on chains. The process takes approximately one week. 

Installing the 25 m absorber allows the near detector to be placed 100 m from 
the target, which would minimize the problem of "parallax." Based on the JAM flux, 
70% of the pions decay before 25 m, but the average energy of the neutrinos from 
those decays is lower. As a result, the event rate with the 25 m absorber installed is 
approximately 50% of the event rate using the 50 m decay path. However, the loss 
in rate is more than compensated by the improved systematics. 

Figure 15.181 shows the angular acceptance for the 100 m detector with a 50 m 
(left) and 25 m (right) decay length. The 25 m decay length has a much narrower 
angular acceptance. The dashed line indicates the acceptance of MiniBooNE for the 
two absorber lengths. Introducing the 25 m absorber results in good agreement. 

As discussed above, the difference in angular acceptance translates into a differ- 
ence in the flux spectrum at FINeSSE and at MiniBooNE. Figure I5.19f left) shows 
the ratio of the FINeSSE to MiniBooNE fluxes for the 50 m absorber, while Fig- 
ure d. 19F right) shows the same ratio for 25 m. The 25 m option is clearly substantially 
better. While some wiggles remain, these are not resolved as a false oscillation signal 
(see discussion of systematic studies below). Figure I5~2U1 shows the parameter space 
covered by FINeSSE with the 25 m absorber configuration (black) compared to the 
50 m configuration (red). The 25 m absorber design maintains the better reach. 

An alternative method for handling the parallax would be to move the detector 
to 200 m. This has the undesirable impact of reducing the event rate by more than a 
factor of four, which limits the other analyses. Also, it reduces the sensitivity to the 
low Am 2 allowed regions of 3+1 models, because the neutrinos would have oscillated 
before reaching the near detector. A comparison of the sensitivity for 100 m and 
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200 m is shown in Figure 13.211 At 200 m, sensitivity to the high Am 2 island is lost. 
Therefore, this is not an ideal solution. 

5.4.3 Systematics Errors and Studies 

Systematic errors from the following sources have been included in the analysis: 

5 e ff : The relative error in the total rates between the near and far detector. This 
is limited mainly by understanding of the fiducial volume of both detectors. 
The understanding of the efficiency of each detector, and the uncertainty in the 
knowledge of L for the 100 m detector, also enter here. We estimate this from 
Monte Carlo studies of the detectors. In our studies, we explore realistic (5%) 
and optimistic levels (2%) for this systematic. 

dnorm- The correlated overall rate normalization of near and far detector. This is 
limited by knowledge of the number of protons on target, and of secondary me- 
son production per proton on target. Since this is not very important for a two 
detector measurement, we have used a conservative estimate of 20% uncertainty 
for these studies. 

$shape' The systematic error associated with the near/far energy difference. This is 
the systematic associated with parallax. It was estimated by varying the beam 
secondary production within the allowed errors. The variation in rate of high 
angle neutrinos in the near detector then translates into this systematic. The 
studies of the variations were done "blind," with one person inserting variation 
and another performing the analysis to identify and minimize the effect by 
changing cuts. Six studies were performed and the range of variation was up to 
20%. 

However, we believe that this error can be substantially reduced by data from 
the HARP experiment |51j and the LMC detector in the beamline. HARP is 
an experiment studying production of secondary mesons by 8 GeV protons on 
replicas of the MiniBooNE beryllium target. HARP results will easily observe 
or eliminate the wildest variations used in the above study. The Little Muon 
Counter (LMC) spectrometer also will provide an in situ test of excess wide- 
angle production. The LMC is located just beyond the 25 m absorber; it consists 
of a long pipe at 7° degrees, followed by a muon spectrometer. If there is an 
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anomalously high rate of high angle secondary production, it will be observed 
in the LMC during the MiniBooNE 50 m running. Assuming HARP and LMC 
limit the variations, a systematic of 10% is obtainable. 

$e scale- The relative energy scale error between the near and far detectors. This will 
be set using the energy distribution of Michel electrons from muon decays as 
a standard candle, along with the ir° mass and cosmic ray muon tracks. We 
believe that a 5% relative error can be straightforwardly achieved. We have 
also looked at the improvement if a 1% relative calibration were to be obtained; 
this, however, would be very difficult. 

In order to cross-compare the sensitivity for various combinations of the system- 
atic errors, we use the 90% CL limit on sin 2 29 for a null disappearance signal at Am 2 
of 1 eV 2 . Table lo~o1 shows the results. The top section of the table applies to 100 m 
running with the 50 m absorber; the middle section shows the expectation for the 
standard configuration of 25 m running. For comparison, the bottom section shows 
the sensitivity if all systematic errors are assumed to be negligible. 

Results for a range of protons on target are shown. From the table, one sees 
that statistics are important, and that in any configuration the results improve with 
further running. Henceforth we will consider only 6 x 10 20 POT, the full FINeSSE 
run. 

From the table, one can see that either the 50 m or 25 m running conditions allow 
for a disappearance search; the 25 m configuration, however, is optimal. Installing 
the 25 m absorber gives a 10% better limit then the 50 m counterpart. Assuming the 
25 m absorber, the most important systematic becomes 5 e jf and not 8 shape, showing 
that the absorber is "doing its job" in reducing the parallax. 

For the bottom line, we assume the following: 6 x 10 20 POT, 5 e ff = 2%, 5 norm = 
20%, 5 s hape = 10%, and S sca i e = 5%. These are the assumptions for the sensitivity 
shown in Figure l2~Hl 

5.4.4 Method for Determining the Sensitivity 

To obtain the sensitivity, we use a standard method applicable to near/far detector 
experiments [73] . A x 2 is formed from Monte Carlo estimates of the number of events 
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Table 5.6: Comparisons of various experimental setups and assumed systematic un- 
certainties. 
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in the near and far detectors without and with oscillations for a given Am 2 and 
sin 2 29. 




The terms involve the events in the far detector binned in visible energy, without and 
with oscillations, N{ ar and J\[° sc -f ar j and the events in the near detector binned in 
visible energy shifted by k Escale , with and without oscillations, N™ ear — ]\[° sc - near . (The 
terms with a tot subscript are the same but summed over all bins.) The oscillation 
event estimate correctly integrates over the length of the decay pipe since this is an 
important effect for the near detector. The first term in the expression is for the 
energy dependent shape analysis, and the second two terms provide the sensitivity to 
a total event counting technique. 

Statistical errors are included along with systematic uncertainties through the 
additional parameters, kj. The kj are systematic error fit parameters and the 8j are 
the systematic errors associated with each of the uncertainties. The final {k/5) 2 terms 
constrain the systematic parameters by their assumed uncertainties. The systematic 
uncertainties, which are described in detail in the following section, include: k norm , the 
overall normalization error common to both detectors; fc e //, the relative normalization 
between the two detectors; k s h ap e, a parameterization (AN- hape ) of a possible energy 
shape difference between the two detectors; and k Egca[e , the uncertainty in the energy 
scale between the two detectors. The energy scale error is included by introducing 
the fit parameter, kEscale, that scales the energy for the near detector binning by 
E' = (1 + k„ ,)E. 

To find the sensitivity of a given experimental setup to oscillations, the above \ 2 
at a fixed Am 2 is minimized with respect to sin 2 29 and all of the k parameters. For 
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k Parameter 


S Values 


Overall Normalization 


knorm 


0.20 


Relative Normalization 


k eff 


0.02 


Shape Uncertainty 


kshape 


0.10 


Relative Energy Scale 


kEscale 


0.05 



Table 5.7: Systematic uncertainties used in the energy dependent sensitivity fits. 

the results shown in this proposal, the fit uses 50 MeV reconstructed visible energy 
bins from 0.0 to 1.5 GeV. Table lo~71 lists the assumed systematic errors for the the 
sensitivity shown in Figure l2~Bl 

5.4.5 Coordination with MiniBooNE 

This analysis requires that the MiniBooNE and FINeSSE collaborations agree on the 
beamline configuration and on sharing data. After discussions, the collaborations have 
agreed to a Memorandum of Understanding, which is designed to prevent conflicts and 
to allow the two experiments to work together. In short, conflicts regarding the beam 
configuration will be resolved by program planning; conflicts regarding data-sharing 
will be resolved by a committee consisting of the co-spokespersons of MiniBooNE and 
FINeSSE. 

We expect this inter-group collaboration to proceed smoothly for three reasons. 
First and foremost, both the FINeSSE and MiniBooNE collaborations regard the 
disappearance search as a high priority. Second, configuring the beam with the 25 m 
absorber may have significant advantages for the Phase II physics of MiniBooNE |74j . 
Thus the MiniBooNE collaboration is expected to be supportive of this configuration. 
However, if there is strong pressure in favor of 50 m running, the disappearance 
search will still be interesting, if not optimal, as discussed above. Third, there is 
substantial overlap in membership between the two collaborations. Thus it will be 
possible for this bi-group analysis to be performed by members of both groups. 
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Chapter 6 
Additional Physics 



While the physics foundations of FINeSSE are As and disappearance studies, there 
is also a wide range of other interesting physics topics which FINeSSE can address. 
For FINeSSE's initial run, these include a significant contribution to cross section 
measurements and an exploration of neutrino magnetic moments. Potential future 
runs of FINeSSE could also permit a measurement of As using antineutrinos, and 
substantial improvement on a — > v e oscillation signal, should one be observed in 
MiniBooNE. We explore these capabilities in this chapter. 

6.1 Cross Section Measurements 

Currently, oscillation experiments rely on modeling of neutrino interactions in a 
regime that is poorly constrained by experimental data. Although accelerator-based 
neutrino beams have existed since the 1970s, our primary knowledge of neutrino in- 
teractions at low energy comes almost entirely from bubble chamber measurements 
made decades ago at ANL, BNL, CERN, and FNAL, all of which were limited both 
by low statistics and by large neutrino flux systematics. In addition to (or perhaps 
because of) these large uncertainties (typically 10 — 30%) the experimental results 
often conflict and are difficult to interpret, mainly because of nuclear corrections and 
exclusive final state ambiguities. 

Improved knowledge of low energy neutrino cross sections will become increas- 
ingly important as experiments move from discovery to precision measurements of 
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oscillation parameters. Consider the following examples. Present atmospheric con- 
straints on Am| 3 and 823 are already limited by flux and cross section systematics. 
Current uncertainties on neutral current (NC) tt° production cross sections currently 
restrict the ability to discriminate between — * v T and — > v s transitions in stud- 
ies of enriched NC samples in the atmospheric neutrino data. Furthermore, — > v e 
appearance searches are limited by the statistical and systematic errors related to 
background subtraction, again, most notably those associated with NC 7r° interac- 
tions where the final state is mis-identified as an electron. Both the kinematics and 
rate of NC 7r° production are less precisely known than most other reaction chan- 
nels, because of the need to model resonant and coherent contributions in addition 
to potential feed-down from inelastic channels. More precise cross section measure- 
ments are not only important for ensuring the success of oscillation measurements. 
Resonant cross sections are particularly relevant for p — > vK + proton decay searches, 
because poorly measured associated strange particle production reactions such as 
v^n —>■ fi~ K + A and v^p — > K + A present significant backgrounds and hence large 
resultant systematics. As a result, while present neutrino experiments could clearly 
benefit from improved knowledge of low energy neutrino cross sections, such advance- 
ment will become more crucial for the success of future neutrino experiments. 

6.1.1 Present Understanding of Quasi-Elastic and Single Pion 
Cross Sections 

Figure 16.11 shows the contributing neutrino cross sections in the region of interest for 
atmospheric and terrestrial based neutrino oscillation experiments. 

At energies near ~ 1 GeV, neutrino interactions are predominantly quasi-elastic (QE) 
and resonant single pion production; each are known to the 10 — 20% level from light 
target (H 2 , D 2 ) neutrino data. Although deep inelastic scattering (DIS) processes 
have been measured with impressive precision (< 2%) at high energies, it is challeng- 
ing both to measure and to model neutrino interactions at low energies, where there 
is substantial overlap between various poorly-understood contributing processes. 

Furthermore, because modern neutrino oscillation detectors consist of heavy nu- 
clei (C, O, Ar, Fe, Pb), the complex target adds additional complications. In this 
case, one must deal with the effects of Pauli blocking, Fermi motion of the target 
nucleons, and final state interactions (i.e. careful accounting for the fact that the 
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Neutrino Energy 

Figure 6.1: Charged current neutrino cross section contributions as a function of 
neutrino energy (in GeV): quasi- elastic (dashed), resonant single pion (dot-dash), 
and deep inelastic scattering (dotted) processes. Figure extracted from JTffi . 

outgoing hadron may reinteract before exiting the nucleus). Final state effects (nu- 
clear reinteractions, ir absorption, and charge exchange) often dominate; they can 
vary depending on the neutrino process, and certainly have not been disentangled 
experimentally. Nuclear effects significantly impact both the rate and kinematics of 
the neutrino reaction, as well as the observed final state event composition and mul- 
tiplicity (see, for example, Section l3~2|) . Although nuclear effects have been studied 
extensively using muon and electron beams, no comparable effort has been made 
using neutrinos. Neutrino cross sections have been measured on nuclear targets in 
the past 1 , but these experiments suffered from low statistics and typically published 
only free nucleon cross sections. By making dedicated, high statistics measurements 
of neutrino interaction cross sections on a scintillator-based target, FINeSSE could 
greatly improve the current experimental situation. 

Several efforts are already underway to more precisely measure neutrino interac- 
tions on nuclei at low energy. Measurements of NC 7r°/QE and inelastic/QE event 
ratios have been performed in the K2K water Cerenkov and scintillator-based fine 
grain near detectors, respectively [TH]. These measurements exhibit ~ 10% accuracy 
based on samples of roughly 5,000-10,000 events [77]. Although MiniBooNE can ad- 
ditionally offer improved cross section constraints with increased statistics over the 

iGargamcllc (C 3 H 8 CF 3 Br), SKAT (C 3 H 8 CF 3 Br), FNAL (Ne), CHARM and CHARM II (mar- 
ble, glass), and Serpukhov (Al) are several examples. 
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K2K near detector ensemble, such Cerenkov-based detection methods are inherently 
limited in their capabilities. The ability to disentangle the various channels (QE, res- 
onant, coherent, DIS, etc.) in a nuclear environment necessitates use of a fine-grained 
detector such as that being proposed for FINeSSE. 



6.1.2 Prospects for Measuring Cross Sections at FINeSSE 

The following subsections outline prospects for several exclusive neutrino cross section 
measurements at FINeSSE. This includes improved constraints on NC ir° production 
and strange particle production. 

Neutral Current ir° Production 

The dominant backgrounds to — > v e appearance searches result from two principal 
sources: the intrinsic v e component in the beam and NC ir° production where the 
final state is misclassified as an electron. Experiments primarily rely on Monte Carlo 
simulations to estimate their 7r° backgrounds. Such simulations must model several 
mechanisms for producing a single ir°: resonance production, coherent single pion 
production, and deep inelastic scattering in which additional hadrons are absorbed 
in the nuclear medium before being detected. The dominant means of single pion 
production at low energy arises through this first production mechanism: excitation 
of baryon resonances (A, N) that decay as: 

v^N -> / N* 

N* — ► 7T N' . 

As a result, there are seven such resonant neutrino reaction channels: three charged 
current and four neutral current: 



-> \i P^ 
u^n — > ji nit 

v a n — > is u nir 

n -> Up p rc~ 



o 



Traditionally, Monte Carlo simulations covering the low energy region have used the- 
oretical calculations by Rein and Sehgal [7H] to predict the rate and kinematics of 
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neutrino resonance production. Such models are tuned to reproduce neutrino single 
pion data, but remain poorly constrained, because of the limited availability and large 
uncertainties in existing experimental data. As an example, Figure shows the ex- 
perimental constraints on single tt° production from charged current (CC) neutrino 
data. Note that all of the data at low energy were collected from light targets. 

CC Single Pion Production 

3. B343, 285 (1990), D 2 
>21 (1979), H 2 , D 2 
1 161 (1982), H 2 , D 2 
2554 (1986), D 2 
527 (1989), CF 3 Br 

) 



E, (G°eV) 

Figure 6.2: Measurements of the CC resonant single tt° cross section, aiv^n — > 
fi~ pn°). Also shown are the Rein and Sehgal-based predictions from two publicly 
available Monte Carlo generators JTjfl. 

By contrast, Figure lB~3*l shows the only available experimental measurement of an 
absolute resonant NC 7r° cross section. These NC data result from a recent reanalysis 
of Gargamelle bubble chamber data at 2 GeV [80J . Because the NC cross sections are 
less well known, experiments typically in practice assign large 25 — 30% uncertainties 
to NC resonant production processes. 
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The data are even more sparse for the other contributing production mechanism, 
coherent single pion production. In coherent interactions, neutrinos scatter off the 
entire nucleus rather than its individual constituents. Because of the negligible energy 
transfer to the target nucleus (A), such processes distinctly provide a single forward 
scattered pion. Like in the resonant case, both NC and CC processes are possible: 

u„A -> u^Air 

Up A — > fl~A7T + 

Almost all current Monte Carlo simulations implement Rein and Sehgal's calcula- 
tion jHl] of coherent pion production cross sections and kinematics. While such 
predictions have been constrained by numerous experimental measurements at high 
energy |ST], the lowest energy data available is at 2 GeV on an aluminum spark cham- 
ber target jH2] • Figure 16.41 shows the low energy Aachen measurement compared to 
several model calculations. 



NC Coherent Pion Production Cross Section 
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Figure 6.4: Experimental measurement of the NC coherent tt° production cross section 
at low energy compared to current model predictions f&^/ ; f&^/. 

The most recent calculations [H3] yield a factor of two to six less coherent pion produc- 
tion at these energies than the earlier Rein and Sehgal based predictions [H3] . Because 
of the lack of low energy experimental data and the existence of several conflicting 
theoretical predictions, oscillation experiments typically assign a 100% uncertainty to 
coherent processes. This large uncertainty is especially important because coherent 
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production may comprise up to 20% of the overall NC 7r° rate. The ability to further 
constrain NC 7r° production at low energies would thus be of great use in achieving 
increased sensitivity to u^ — > v e oscillations, and in placing more stringent limits on 
the oscillation of standard neutrinos to sterile states. 

In two years of running, FINeSSE will collect a total sample of ~ 20, 000 NC 7r° 
resonant events and over 5,000 NC coherent ir° interactions. In addition, FINeSSE's 
superior energy resolution and event reconstruction capabilities will greatly enhance 
the ability to select ir° interactions. Figures lfi~5l and IB~7I show simulated NC resonant 
7T° events in the FINeSSE detector. A u^p — > u^pir interaction fFigure l6*3j) can be 
distinguished by the presence of three separated energy deposits corresponding to the 
final state proton and the two photons emitted from ir° — > 7 7. A charged current ir° 
p~ pir°, ( Figure l6~B|) is additionally accompanied by a final state muon 
Up n 7T° interaction (Figure 16 .7|) can only be produced by 
scattering from carbon, and will contain only two clusters of hits corresponding to 
the two photons from the 7r° decay. NC coherent ir° events (Figure l6~8|) are similar in 
signature to the former class of events; however, in this case, the energies and angles 
of the two final state photons can be used to determine if the ir° angular distribution 
is more forward peaked as one expects for coherent scattering. 
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Figure 6.5: Hit-level simulations of NC tt° interactions in the FINeSSE detector, 
Uppn . The three hit clusters correspond to the final state proton and two 
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Figure 6.6: A v^n — > pT pir° CC it interaction in the FINeSSE detector. The right 
hand figure shows the true GEANT particle trajectories in the YZ plane. In this 
case, a muon is produced in addition to the three hit clusters from the proton and two 
photons. 
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Figure 6.7: Two iVC tt interactions in the FINeSSE detector, v^n — * u^nir . 
separated hit clusters correspond to the two photons emitted from the tt° decay. 
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Figure 6.8: A NC coherent ir° interaction in the FINeSSE detector, l2 C 



12 



C 7T° . The figure on the right shows the true GEANT particle trajectories in 
the XZ plane for the same event. As can be seen, the hit clusters correspond to the 
two photons from the 7r° decay. 



While such events leave distinct signatures in the FINeSSE detector, determining 
the efficiency and accuracy with which FINeSSE can measure 7r° rates and kinematics 
will require development of a cluster algorithm and tracker designed to located gaps 
in energy deposits. While not yet optimized, preliminary studies have demonstrated 
that 7T° efficiencies and purities can be obtained that surpass those achieved with 
either the MiniBooNE or K2K Cerenkov-based detectors. Given that the largest 
sources of error arise from nuclear uncertainties, by measuring 7r° production on 
carbon, FINeSSE should therefore be able to reduce the currently assigned NC 7r° 
production uncertainty (Section 16.1.2)1 by as much as a factor of two. 

In its ability to isolate 7r° interactions, FINeSSE has an especially important 
role to play given that there is currently no NC 7r° data on nuclear targets below 
2 GeV, and only low energy CC tt° data measured on deuterium targets. FINeSSE 
is well-positioned to be the first neutrino experiment to measure NC tt° (and CC 
7T°) production on carbon at these energies (< E u >~ 700 MeV), thus providing 
an important constraint to oscillation experiments such as MiniBooNE, as well as 
to other accelerator- and atmospheric-based v e appearance experiments employing 
heavy nuclear targets. Such efforts at FINeSSE are also complementary to slightly 
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higher energy (< E v >~ 1.3 GeV) scintillator-based analyses currently being planned 
after upgrades to the K2K near detector site [SB] . 

Strange Particle Production 

Proton decay modes containing a final state kaon, p — > uK + , have large branching 
ratios in many SUSY GUT models. Because there is a non-zero probability that an 
atmospheric neutrino interaction can mimic a proton decay signature, 

v^n — > pT K + A 
v»P -> v^K + A, 

it is important to reliably estimate this background. Present uncertainties on this 
background process are as large as 100% [SH] both because neutrino strange particle 
production rates have been measured in only a few experiments, and because there 
are not many predictive theoretical models [HZ]- Figure lfi~T?l shows the only two exper- 
iments which have published cross sections on the dominant associated production 
channel, v^n — > p~K + A. Both measurements were made on a deuterium target and 
based on less than 30 events combined. 
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Figure 6.9: Experimental measurements of the v^n — > p~K + A cross section. Also 
shown are the predictions from two publicly available Monte Carlo generators JTffi . 

Based on NUANCE Monte Carlo predictions, approximately 90 such interactions 
are expected at FINeSSE in two years of running. The ability to identify and isolate 
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neutrino production of strange particles at FINeSSE would be of even larger impor- 
tance to proton decay searches, not only because of potentially large event samples, 
but also because these interactions will be taking place at energies near threshold and 
on a nuclear target, both of which are particularly relevant for present and future pro- 
ton decay studies. Identifying such events at FINeSSE will certainly be challenging. 
Although the event signature is complex, it is potentially very distinctive. Roughly 
64% of the time, a typical CC kaon event will contain two oppositely charged muons, 
one from the primary neutrino interaction and the other from K decay, K + — > 
fFigure fo.lO)) . In addition, since most of the kaons will stop and decay at rest, the 
fi + is monoenergetic with a momentum of ~ 240 MeV, and will appear roughly 12 ns 
after the initial interaction. 
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Figure 6.10: Example of a v^n — > fi K + A event in the FINeSSE detector. In this 
event, the associated strange particles decayed yielding a hadronic shower and second 
muon: h.—>pix~ (BR=63%) and K+ -> fi+u^ (BR=G4%). 



Further studies must assess FINeSSE's multi-track reconstruction capabilities and 
DIS backgrounds before any stringent limits on strange particle production can be set. 
However, such a possibility should be seriously explored as FINeSSE could potentially 
be the first experiment to measure this reaction both near threshold and on a nuclear 
target. 
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6.2 Neutrino Electron Scattering: Neutrino Mag- 
netic Moment 

In the Standard Model, massive, charged particles have intrinsic magnetic moments 
by virtue of their spin. If we introduce neutrino mass into the theory, an effective 
neutrino magnetic moment can arise. How neutrino mass fits into the Standard 
Model, be it via a standard Dirac mechanism, SUSY, or something else, affects the 
size and origin of this magnetic moment. Thus, we can use a neutrino magnetic 
moment measurement to tell us about how neutrinos fit into the larger theory. 

In the minimally extended Standard Model, massive Dirac neutrinos of mass m v 
can have a neutrino magnetic moment of: 



arising from one loop radiative corrections in diagrams with W-boson exchange as in 
Figure 16.111 

We can place neutrinos with non-zero magnetic moments in the context of a 
larger theory such as SUSY and Extra Dimensions. These extensions to the Standard 
Model such as the super symmetric left-right model predict larger neutrino magnetic 
moments of the size |HS] 



Figure 6.11: One loop radiative correction diagram with a massive neutrino. 
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Theories involving Extra Dimensions produce an effective neutrino magnetic moment 
as large as 10 _11 /i£ jHHj- 

A non-zero neutrino magnetic moment would also have important implications in 
cosmology in the development of stellar models. Astrophysical limits such as plasmon, 
or photon wave packets, decay rates from horizontal branch stars and neutrino energy 
loss rate from supernovae allow a neutrino magnetic moment as large as 10 _11 /iB 



The most straightforward way to look for a non-zero neutrino magnetic moment 
is via an excess in v — e scattering: 



A non-zero neutrino magnetic moment will give rise to an electromagnetic contribu- 
tion to this process in addition to the Standard Model Z exchange. This electromag- 
netic contribution to the cross section for neutrino-electron elastic scattering is of the 
form: 



where r is the classical electron radius, T e min is the minimum kinetic energy of the 
recoil electron, E v is the neutrino energy, and / is the neutrino magnetic moment in 
units of Bohr magnetons. 

At low y = T e /E u , the electromagnetic contribution to the neutrino-electron 
cross section increases rapidly, while the Standard Model contribution increases only 
gradually. The resulting shape dependence in the differential cross section can be 
used to look for a signal in a high statistics experiment like FINeSSE. The differential 
cross section contributions for weak and electromagnetic components of the neutrino- 
electron elastic scattering cross section are given by: 
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v b lavor 


Limit bet{HB) 


Experimental Data 


Process 


T > C 

Reference 




1.5 x lCT 10 


Super Kamiokande 


shape of differential 
cross section 


m 


V~e 


1.8 x 1CT 10 


combination of V e re- 
actor experiments 


excess in total cross 
section 




V * 


6.8 x lCT 10 


LSND 


excess in total cross 
section 


[SB! 


V T 


5.4 x 1CT 7 


DONUT 


excess in total cross 
section 





Table 6.1: Limits set on neutrino magnetic moments. 



where T e is the electron recoil energy and E v is the neutrino energy. 

Present limits set on neutrino magnetic moments are many orders of magnitude 
away from the Standard Model prediction including massive Dirac neutrinos. They 
are, however, only one to two orders of magnitude away from predictions from certain 
beyond-the-Standard-Model theories and astrophysical limits. The most stringent 
limits set by experiments for neutrino magnetic moments for each neutrino flavor are 
listed in Table 16.11 The experimental limit for the muon-neutrino magnetic moment 
coming from the LSND experiment is an upper limit of \i v < 6.8 x 10 _10 /iB set 
by measuring the total neutrino-electron elastic scattering cross section. An improved 
limit by about a factor of two, on the muon neutrino magnetic moment is expected 
from MiniBooNE [§E]. 

Number of events and, more importantly, low energy threshold for electron recoil, 
determine how sensitive an experiment is to a non-zero neutrino magnetic moment. 
FINeSSE expects a number of v — e scattering events (86 events) comparable to the 
number MiniBooNE will see (101 events). FINeSSE can improve upon MiniBooNE's 
expected measurement in identifying events at low electron recoil energy, where the 
differential cross section is most sensitive to a non-zero neutrino magnetic moment. 
Figure 16.121 shows a typical v — e interaction in FINeSSE's Vertex Detector. These 
events are very forward and have a very distinctive shower shape. The granularity of 
the Vertex Detector as compared to the MiniBooNE detector may allow FINeSSE to 
clearly identify these lowest energy events via a shower shape analysis. This would 
provide better sensitivity to non-zero muon neutrino magnetic moment than any 
previous experiment. 
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e elastic scattering interaction in the FINeSSE detector. 



The figure on the right shows the true electron shower in GEANT for the same event. 
Such events characteristically contain a single forward scattered electron. 
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6.3 Antineutrino Running and As 

If the opportunity arises to run FINeSSE with an antineutrino beam, there are some 
features of up — > up scattering that would make this channel exciting to investigate. 
The ratio in the antineutrino channel is actually slightly more sensitive to As than 
for neutrinos. This can be seen in Figure I6.13b -b. where the relative change in the 
ratio with the same change in As is larger for antineutrinos than for neutrinos. In 
addition, the change in the ratio as a function of Q 2 is different for antineutrinos. 
This would allow a simultaneous extraction of As and Ma, since (see Figure 16.13)) 
the ratio has a much different behavior in Q 2 for different values of Ma in up — > up. 

6.4 Long-term Running of FINeSSE+MiniBooNE 

If MiniBooNE sees a signal and it is confirmed by both phase II MiniBooNE and 
FINeSSE+MiniBooNE cross checks, Fermilab is likely to continue the 8 GeV program. 
The natural next step is the installation of BooNE, the multi-detector experiment 
planned as the upgrade to MiniBooNE. Given that BooNE is in the planning stages, 
it is unclear when this program would begin. However, it is safe to say that it is 
unlikely that BooNE can be fully constructed and running before 2010. During this 
construction period, it would make sense to continue running FINeSSE+MiniBooNE. 

In this section, we consider what can be achieved if FINeSSE+MiniBooNE re- 
ceives 3 x 10 21 protons on target before the BooNE detector system is inaugurated. 
This is based on the premise that upgrades to the Booster can allow the 8 GeV line 
to return to the level of 1 x 10 21 protons/year near the end of the decade. To be clear, 
we are not requesting approval for this scenario from the PAC. We are simply using 
this to illustrate the long-term value of the FINeSSE detector. 

In this case, the statistics limited measurements described in the previous section 
can be greatly improved. We believe we can reasonably achieve the systematic errors 
listed in Table lfi~4l 

At this point MiniBooNE will have acquired 4 x 10 21 protons on target (25% with 
the 50 m absorber and 75% with the 25 m absorber). The capability of such a 
measurement is shown in Figure 16.141 The colored regions show the LSND allowed 
range. The solid ellipses indicates the la measurement capability for the oscillation 
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NC 7T°: 0.03 

Radiative A Decay: 0.06 

v e from /i: 0.04 

v e from K+: 0.04 

u e from K°: 0.06 

Table 6.2: Systematic errors for the (hypothetical) scenario of 3 x 10 21 POT. 



parameters at two possible true Am 2 values. For comparison, the capability from 
MiniBooNE phase I running is indicated by the dashed ellipses. 
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Figure 6.13: The ratio of NC to CC scattering for neutrinos (left column of plots) 
and antineutrinos ( right column ) as a function of Q 2 . In each of the 4 rows one of 
the form factors is varied. In a-b) -0.2 < As < +0.2, c-d) -0.2 < F( < +0.2, e-f) 
0.95 < Ma < 1.15 GeV/c 2 . 
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Figure 6.14: The solid ellipse shows the la measurement capability for two sets 
of oscillation parameters for FINeSSE+MiniBooNE for 3 x 10 21 POT. The dashed 
ellipse shows the MiniBooNE Phase I capability. The colored segments indicate the 
LSND-allowed regions. 



Chapter 7 



Overview of Cost and Schedule 



The FINeSSE detector cost is estimated to be $2.25 M ($2.8 M with contingency). 
These estimates do not include EDIA or indirect costs. Funding for the detector will 
come from the university funding agencies. EDIA and indirect costs will be added as 
appropriate when the funding structure is in place, and will account for differences 
in each university's accounting practices. 

The FINeSSE enclosure cost is estimated to be $800K ($1.6 M including con- 
tractor O&H, EDIA, management reserve, and indirect costs) as determined through 
a Preliminary Design Report written by the Fermilab Facilities Engineering Section 
(FESS). See Appendix IB1 which contains this report. 

The dimensions of the enclosure were chosen to comfortably contain the detector 
and provide adequate working room around the detector. There is also adequate 
room for detector installation. The FINeSSE schedule incorporates detector and 
enclosure design, construction, and installation over the next two years, with data 
taking scheduled to begin in mid-2006. 

The following sections provide breakdown of detector costs, enclosure costs, con- 
struction and installation schedule, and an overview of the FINeSSE schedule from 
now until the beginning of the FINeSSE physics run. 

7.1 Detector and Enclosure Costs 

The FINeSSE detector is comprised of two subdetectors (cf. Chapter HJ): the Vertex 
Detector, followed downstream by the Muon Rangestack. Table 7.1 lists preliminary 
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costing of the components of these detectors, as well as total costs. 

The cost of the detector is driven mainly by PMT and electronics costs that 
combined total $900,000. The vertex detector tank that holds the scintillator, fiber, 
and PMTs is about $265,000. The remaining costs are for fiber, scintillator, and oil. 
The cost driver for the Rangestack detector is the iron planes, costing $327,000. Most 
of this cost comes from fabrication. An additional cost of about $100,000 is set aside 
for rigging and installation of the detector. 

Contingency is applied to the total detector costs. EDIA and indirect costs are 
not included, as it is expected that the majority of these costs will come through 
university groups, where EDIA and indirect costs are treated differently by each 
university. These costs will be included once the funding profile by universities is 
established. 

The detector enclosure design went through a series of studies within FESS. 
To meet the requirements of detector construction and operation, several enclosure 
options were studied. They are summarized in the Alternate Comparison Matrix 
given in Appendix IBI The options ranged in cost from $1.48M to $2.54M. The study 
compared the use of a gantry crane to a bridge crane, a hatch cover to a pre-engineered 
building, and a sheet pile pit to a secant pile pit. The chosen option includes a sheet 
pile pit, a hatch cover, and a gantry crane for about $1.6M. If a mobile crane were 
used, rather than a gantry crane, the construction cost could be reduced to $1.48M, 
but this would be offset by the operational costs of crane rental during the detector 
installation phase. There will be no need for crane coverage once the detector is 
installed and operating. A cross-sectional view of the enclosure is shown in Figure 4 
of the Appendix, and the plan view is given in Figure 3. 

7.2 Schedule 

The FINeSSE schedule allows for 2.5 years for enclosure and detector construction and 
installation, with the start of the physics run beginning in mid-2006. The FINeSSE 
schedule for enclosure construction calls for construction of the building to commence 
in mid to late 2004. Total time to construct this enclosure is about six months, 
as described in detail in Appendix |Bl keeping FINeSSE on schedule for detector 
installation beginning in mid-2005. Detector construction and installation is expected 
to take 1.5 years. 
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FINeSSE is on an aggressive but achievable schedule for a number of reasons. 
First, FINeSSE physics is important and timely and should be pursued as quickly as 
possible. Second, there is an existing, running neutrino beamline which should be 
utilized as much as possible. Finally, the FINeSSE physics run is timed to coincide 
with the end of MiniBooNE antineutrino running and fits well with the MiniBooNE 
schedule. The FINeSSE schedule, up to the start of our physics run in mid-2006, is 
shown in Figure [7~T1 
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Figure 7.1: TTie FINeSSE schedule, showing the milestones from winter 2004 until 
beam can be taken in 2006. Milestones are shown in red. 
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Chapter 8 
Conclusions 



We propose to construct and operate a novel neutrino detector to run in the Booster 
neutrino beamline, 100 m from the Booster neutrino production target. The physics 
motivation for this experiment is twofold. 

First, we will measure the strange spin of the proton As using a theoretically ro- 
bust method, an intense, low energy neutrino beam, and a novel detection technique. 
Second, FINeSSE will look for disappearance, in conjunction with MiniBooNE, at 
high Am 2 in an astrophysically interesting region. 

FINeSSE will achieve these physics goals using a two-part detector to measure 
both short, low energy proton tracks and longer, muon tracks. Low energy proton 
tracks will be well measured in a, liquid scintillator "bubble chamber" . Downstream 
of this Vertex Detector is a Muon Rangestack designed to range out higher energy 
muons that are not contained in the Vertex Detector. 

FINeSSE requires 6xl0 20 POT total to reach its physics goals. In the NuMI era, 
these protons can be delivered over a two year period. 
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Appendix A 

Vertex Detector Design Studies 



In determining the best detector design to address the FINeSSE physics goals, a 
number of options were considered. In addition to the "Scibath" design described 
in Chapter 01 we studied a design called "Scistack", much like the K2K Scibar 
detector [HO]- We found this detector was not optimal for identifying low energy 
proton tracks from v — p elastic scatters, crucial for the As analysis. By contrast, the 
Scibath option did not have these limitations. This appendix describes these Scistack 
studies including building and testing of a prototype detector, Monte Carlo studies, 
and conclusions. 

A.l Scistack design, Feasibility and Cost 

The Scistack detector is a 2.5 m 3 volume of scintillating plastic bars. The 2 cm x 1 cm 
extruded polystyrene bars are co-extruded with a Ti02 outer layer for reflectivity, and 
a hole down the center for a WLS fiber. Bars are organized in planes normal to the 
beam direction with bars in even-numbered planes oriented in the x direction and 
bars in odd number planes oriented in the y direction. This design has a total of 
31250 bars and associated readout channels. Readout is identical to the Scibath 
design with light piped out via WLS fibers, and recorded and digitized as described 
in Section 14.1.31 

The cost of a Scistack detector is comparable to the Scibath detector. Plastic 
and liquid scintillator as well as WLS fiber are comparable in both detectors. A cost 
driver for both detectors is channel cost. The number of channels in this design is 
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larger by a factor of 1.5. This is offset by employing channel multiplexing which 
decreases the channel cost by a factor of eight. However, channel multiplexing can be 
accommodated only with 16 anode MAPMTs (as opposed to the Scibath 64 anode 
MAPMTs) due to physical size of each anode, driving the cost per channel back up. 
Cost of WLS fiber also increases in Scistack as there are 1.5 more channels; this, 
however, is not a significant cost overall. Folding in all cost differences, costs for the 
two detectors are comparable. 

A. 2 Prototype testing 

A small prototype of the Scistack detector was designed, constructed, and tested 
during the summer of 2003. The prototype design, shown in Figure IA~T1 consists of a 
ten by three stack of 2 cm x 1 cm scintillator bars, with a WLS fiber inserted down 
the center of each. The fiber is attached to 16 anode MaPMTs and read out through 
the same electronics used for the Scibath prototype tests, described in ChapterHJ The 
Scistack prototype was exposed to the IUCF proton beam, as was Scibath (again, see 
Chapter EJ). Scistack performed as expected. Light output increased as the proton 
beam moved, normal to the bars, and closer to the readout. Tracks were easily seen, 
as expected, as the proton beam traversed the stack both with the stack oriented 
normal to the beam and at an angle with respect to the beam. This detector employs 
known technology; there were no surprises. 

A. 3 Monte Carlo Studies 

While the detection techniques proved acceptable experimentally, Scistack Monte 
Carlo work revealed a feature in the detector geometry which makes this detector 
unsuitable for FINeSSE's As measurement. 

The design goal of the FINeSSE Vertex Detector is to be able to cleanly identify 
and measure low energy protons from v — p elastic scatters. Scistack Monte Carlo 
work shows that these very events, most important to the As analysis, are frequently 
not reconstructible in Scistack due to limitations in tracking in this detector. Specif- 
ically, low Q 2 events tend to be directed at high angles with respect to the beam 
direction. The Scistack detector reconstructs these events with less success, since 
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Figure A.l: Photograph of the Scistack prototype. Scintillator bars are stacked in a 
10 x 3 array. WLS fiber, polished on both ends, routed through each bar, is read out 
on one end by Hamamatsu 16 anode MaPMTs. The entire prototype is housed in a 
light tight aluminum box with connections for Lemo and high voltage. The box is built 
on a rotatable stand to facilitate testing in the IUCF beam. 
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bars are oriented normal to the beam direction, and therefore in the same direction 
as these low Q 2 events. The Monte Carlo and these studies are described below. 



A. 3.1 Scistack Monte Carlo 

The GEANT-based Monte Carlo for the Scistack option is designed to be identical 
to the Scibath option except for the configuration of the central detector. This in- 
cludes the output ntuple which has the same format, with appropriate changes in the 
meaning of the variables, to minimize differences in the reconstruction program. 

The central detector for Scistack consists of 250 layers, each 1 cm thick, along the 
beam axis (z). Each layer has 125 scintillator bars, 2 cm in the transverse direction, 
with the bars in alternating layers aligned along the x and y axes. Each 1 cm by 2 cm 
by 2.5 m bar has a 1.5 mm OD wavelength shifting optical fiber along its axis, to 
capture and transmit photons from the scintillator to the PMTs. The outer 0.5 mm 
of each bar is a Ti02 wrapper which optically isolates the bars. The highly reflective 
inner surface of the wrapper increases the light capture in the fibers by approximately 
a factor of two. The scintillator has an attenuation length of 5.0 m. 

After allowing for fiber capture and transmission efficiencies, assumed to be the 
same as Scibath, each photon entering a fiber has its layer and bar number recorded 
in the output ntuple. The actual capture efficiency as well as the number of photons 
generated per MeV of energy loss would be higher in Scistack than in Scibath. But 
the number of photons detected is so large that this is not an important factor in the 
current studies. 

The resulting ntuple is reconstructed by the same Hough transform method used 
for Scibath, and the results are used to determine the angular and energy resolution 
of Scistack as well as its particle identification capabilities. The Monte Carlo samples 
generated for these studies include single electrons, muons, nucleons, pions, and kaons 
with kinetic energies between and 500 MeV, generated isotropically in direction. 
We also ran a large sample of neutrino interactions generated using NUANCE. Each 
sample was run with vertices near the center of the detector, and again with vertices 
spread uniformly throughout the detector. The results are discussed in the following 
section. 
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Figure A. 2: Number of tracks reconstructed versus true kinetic energy for protons 
generated within Scitrack. Events with kinetic energy below 0.15 GeV are not usually 
reconstructible (have no tracks). 



A. 3. 2 v — p Sample in Scistack 



Neutrino-proton elastic scatters at low energies produce short-track, low energy pro- 
tons. Samples of protons generated and reconstructed as described above were studied 
to determine if the Scistack detector was suitable for this measurement. Figure lA~2l 
shows the number of tracks reconstructed as a function of true kinetic energy for this 
sample. As is indicated, tracks below 150 MeV tend not to be reconstructible. These 
tracks tend to be at high angle and therefore traversed few bars. Figure IA.3I show 
reconstructed tracks as a function of true kinetic energy and angle for both Scitrack 
and Scibath. 

The neutrino interaction sample generated with NUANCE was used to under- 
stand the Scistack efficiency as compared to Scibath. Similar cuts as were applied 
to the reconstructed Scibath events to extract a NCp sample, were applied to the 
Scistack sample. The inability to reconstruct the proton events described above, 
translated into a more than 30% decrease in v — p statistics in the 0.2< Q 2 <0.4 bin, 
the crucial energy bin for the As analysis. 
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Figure A. 3: Reconstructible events versus true angle and energy in Scistack versus 
Scibath. Upper left corresponds to true kinetic energy of Scibath event; upper right 
to true kinetic energy of Scistack protons. Bottom left and right are true angle for 
Scibath and Scistack respectively. 
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A. 4 Conclusions 



While the Scistack detector technology is proven experimentally, it is not suitable for 
FINeSSE's As measurements. Scibath, by contrast, has the ability to reconstruct 
even low Q 2 events, crucial for this analysis. 
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Appendix B 



FESS Project Definition Report 
The FINeSSE Detector 



Fermilab 




November 6, 2003 



MEMO TO: B. Fleming, Particle Physics Division 



FROM: 



T. Lackowski, FESS/Engineering 
S. Dixon, FESS/Engineering 



SUBJECT: Project Definition Report 
FINeSSE Detector 

FESS/Engineering Project No. 6-7-59 



Per your request we have investigated the scope, cost and schedule required to 
provide adequate housing for the FINeSSE Detector north of the Booster Neutrino 
Target Hall. The design solution provides for a braced excavated underground 
enclosure that provides the space required to install and operate the detector. 
This report presents the developed drawings, preliminary budgetary cost 
estimates and project schedule for the conventional construction portion of the 
work. 



The objective of this project is to provide a cost effective design solution to 
provide the conventional construction required for the FINeSSE detector. Both 
initial and life cycle costs were considered in the design process. The detector's 
installation, operation and maintenance objectives formed the basis for the 
project's criteria. 

Project Criteria: 

The following describe the FINeSSE requirements used in the development of this 



• The FINeSSE detector is comprised of two parts. The upstream part is 
the vertexing detector. Located immediately following the vertexing 
detector is the muon rangestack. 

• The vertexing detector is a 3.5m x 3.5m x 3.5m (11'-6"x11'-6"x11-6") 
module located on the axis of the Booster Neutrino beamline 

• The vertexing detector is located 100 meters downstream from the 
production target located in the Booster Neutrino Target Hall. 

• Elevation of the Booster Neutrino Target is 723'-0" 

• The vertexing detector is centered at elevation 724.15'. This is the 
beam elevation 100 meters along a line connecting the Booster 
Neutrino target (elevation 723') and the center of the Booster Neutrino 
detector sphere (elevation 729.22'). 



Objective: 



study: 
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• The muon rangestack is located directly downstream of the detector. 
This equipment is 4m high x 4m wide x 3m long (13'-2" x 13'-2" x 9'- 
10"). 

• The muon rangestack is constructed of plastic scintillator and steel 
plates and weighs approximately 120 tons. 

• Removable concrete shield blocks will be installed over the detector to 
reduce cosmic and atmospheric particles interacting with the detector. 
The shielding will be comprised of existing concrete shield blocks with a 
depth of 3'-0". The detector is fully accessible during beam operations. 

• Electronic components related to the detector will be located on a 
platform above the detector. The components will be installed in a 
standard computer rack (2'-6" x 2'-6" x 7'-0" h). Racks will be installed 
in environmentally controlled cabinets. 

• No specific environmental control is required. Electric heat is provided. 

• Electrical power requirement should provide power for four (4) typical 
electronic racks and associated communication equipment. 

• An oil cooling system will be required for the vertexing detector. The 
pump, reservoir and associated piping will be located adjacent to the 
detector. 



Figure 1 (below) is a depiction of the scope of the FINeSSE detector. 




Figure 1 - FINeSSE Detector 
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Project Location: 

The project will be located 100 meters (328 feet) downstream from the Booster 
Neutrino production target. Figure 2 {bdow) indicates the project location and 
adjacent faciities. 




Fifluift 2 - 5*6 Laauon Plan 



The project $ite is outside of known wetlands a$ indicated in the Wetland 
Defi neation and Special Areas Report completed by Con soar Town send 
Ehvirodyne Engineers, Inc. in January 2002. However, additional verification and 
field investigation may be required should this project proceed to construction. 
All required NEPA documentation wil be obtained. This project wil disturb 
greater than 1 acre of surface area. 



Alternates 

FtNe-SSE b&ieciof, Pfoj&ci Defifiiticm Report 



Several possible design studies were produced during the preparation of this 
report. These alternates are contained in the attached appendix. 

Project Description: 

The civil construction portion of the project incorporates conventional construction 
methods to provide the spatial requirements for the physics apparatus. In 
general, the construction will consist of a below-grade enclosure housing the 
physics apparatus and related electronics while the above-grade portion will 
consist of a metal hatch cover and gantry crane. No above-grade building is 
anticipated. Detailed descriptions of the individual components are listed below. 

The proposed project site will be accessed from both the existing Booster 
Neutrino Target Hall parking lot and Giese Road. Construction traffic will be 
limited to the Giese Road access. The design includes the road extension and 
limited staging area. The road extension and staging area will be a gravel 
hardstand constructed utilizing salvaged rock material from the NuMI stockpile 
located at Site 12. The hardstand will be sized to allow truck access to the 
loading/unloading area beneath a 10-ton gantry crane. 

In order to provide a cost effective excavation and enclosure a braced excavation 
utilizing steel sheet piles will be installed. This sheet piling system will serve as 
the permanent wall system. The sheet pile system will be installed on three (3) 
sides of the excavated area. Conventional excavation methods will be used to 
remove the material within the braced excavation. This design assumes that 
conventional excavation equipment will utilize a ramp along the north side of the 
enclosure. The fourth side of the enclosure will be constructed on cast-in-place 
concrete after excavation is complete 

Figure 3 (shown below) is a plan view of the proposed enclosure: 




Figure 3 - Enclosure Plan 



FINeSSE Detector, Project Definition Report 



Page 4 of 7 



The design provides for the installation of 24' long precast concrete shield blocks. 
These shield blocks are assumed to be available from existing Fermilab sources 
(KTeV railhead, etc.) and will be installed above the detector components. The 
shielcfbfecks will be installed in such a manner to be removable for access to the 
detector. A removable metal hatch cover will be installed above the shield blocks 
to provide a weather barrier. The hatch cover will be installed on rails to allow it to 
be removed for access to the shield blocks and detector components. 




Figure 4 - Section Looking West 



The project will utilize a ten (10) ton, rail mounted gantry crane for material 
handling. This crane will be used to move the metal hatch, install and remove the 
precast concrete shield blocks and assemble the detector components. The 
gantry crane will be rated for exterior use and able to withstand the environmental 
conditions at Fermilab. 

This project will obtain Industrial Cooling Water (ICW) for fire fighting purposes by 
extending the existing service from the Booster Neutrino Target Hall site. Fire 
detection systems are included in the scope of work as well as appropriate 
automatic sprinkler systems are included. 

Electrical power will be extended from the Giese Road Substation at 13.8kv and 
transformed to 480V at the project site. 

Communication will be extended from the MI-12 Service building via a 5" conduit 

The attached PDR drawings provide a depiction of the design at this stage of the 
project. Addition refinement during subsequent phases will affect the overall 
project scope and cost. 
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Estimated Project Schedule 

Listed below is a proposed schedule for the remaining activities for this project: 




These durations and activities are preliminary and reflect typical durations for a 
protect of this scope. Some refinement and optimization is possible and will be 
investigated in subsequent phases. 



Estimated Protect Costs 

Listed below is a breakdown of the expected costs for this project assuming the 
entire scope of work described above is procured as one construction package: 



Construction Total $794,000 

Overhead and Profit 17% 3 135,000 

Subtotal 3929,000 

ED LA 19% 3177.000 

Subtotal 31,106,000 

Man. Reserve 22% 3243,000 

Subtotal 31,349,000 

bufrect Costs 19% 3256.000 

iTAIATED PROJECT TOTAL *1 t fl0S + 0M 
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Cost Estimate Basis 

The attached PDR drawings and detailed cost estimate describe the known scope 
of the conventional construction portion of the project. Physics apparatus and 
related S( Bxperimental equipment as assumed to be provided and installed by 
others means. 

The costs contained in this Project Definition Report are based on FY2004 dollars. 
Appropriate escalation will need to be applied once a funding source is identified. 

The included budgetary cost estimate is based on cost data taken from Means 
Cost Estimating Guides, historical data and recent construction history here at 
Fermilab. While the budgetary cost estimate can provide input for the feasibility of 
the project, further design refinement will affect the final cost of the project. 

Engineering Design Inspection and Administration (EDIA) costs are included on 
the above listed budget estimates. ED&I activities include the engineering and 
design activities in Titles 1 and II, the inspection activities associated with Title III. 
The descriptions are based on DOE Directive G430.1-1, Chapter 6. Administration 
activities include those defined by DOE Directive G430.1-1, Chapter 6 as Project 
Management (PM) and Construction Management (CM). Past historical data and 
DOE Directive G430.1-1, Chapter 25 indicate that 18%-25% of the construction 
costs is an appropriate range. 

Management Reserve costs are included in the above budget cost estimates. 
Based on DOE Directive G430.1-1, Chapter 11 DOE guidelines and the pre- 
conceptual nature of the design at this stage a Management Reserve of 20%-35% 
of the above costs is considered an appropriate range. 

The costs include Indirect Cost multipliers. Indirect Costs rates are defined by 
DOE Order 4700.1 that states indirect costs are "...costs incurred by an 
organization for common or joint objectives and which cannot be identified 
specifically with a particular activity or project. If this work should become a GPP 
project, Indirect Costs will have to be applied, but the amount will be affected by 
the rates in effect at the time this project is initiated. Currently, the Indirect Cost 
multiplier is 19% of the above costs. 

Summary 

This memo provides a budgetary cost estimate based on preliminary designs and 
conversations and is intended for planning purposes. The information produced is 
subject to refinement during the design process. 

If you would like to proceed with this project, please contact Tom Lackowski 
(x3640) or Steve Dixon (x8501). 

End: Appendix - Design Alternates 
PDR Drawings 1 through 5 
Cost Estimate, dated 11/06/03 (3 sheets) 
Existing MiniBooNE coordinates 
Gantry Crane Estimate, dated 

cc: J. Cooper, PPD (w/encl.) M/S 208 R. Stefanski, PPD (w/encl.) M/S 122 H. Jostlein PPD (w/encl.) M/S 208 
D. Nevin, FESS (w/encl.) E. Crumpley, FESS/E (w/encl.) R. Ortgiesen, FESS/IM (w/encl.) 

R. Alber, FESS/E (w/encl.) Project File 6-7-59 
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Appendix 



This appendix contains descriptions and cost estimates for alternate designs for 
the FINeSSE detector enclosure contained in the Project Definition Report.181 
These designs have not been fully explored and do not reflect a completed 
design. They are included in the appendix for reference only. 



Alternate No. 1 - Conventional Building 

The proximity of the MiniBooNE Target Hall and other Main Injector buildings 
could provide the design basis for the FINeSSE Detector building. This 
alternative provides a similar building style. 

In order to provide a cost effective excavation and enclosure a braced excavation 
utilizing steel sheet piles will be installed. This sheet piling system will serve as 
the permanent wall system. 



PZ33 SHPfcT PILING 
WITH BRACED -. 

FM1IE 



mm 

111 




23 CLEAR 



Enclosure Plan 
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The floor will be a cast-in-place concrete slab with underdrains to collect and 
transit groundwater to a sump pit for discharge. The space will be accessed via 
a steel stair. Electric unit heaters will prevent freezing, but no conditioning of the 
air is required. 

The below grade enclosure will provide space for the detector and muon 
rangestack and associated electronic equipment. 

The above grade structure will provide a loading dock, overhead bridge crane 
and space for the electronics caboose. The building will be similar in style and 
appearance to the nearby MiniBooNE Target Hall, including pre-finished metal 
wall panels, doors and roofing materials. 




Plan at Grade 
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The cross section shown below indicates the general arrangement of the facilities 
and physics components. 183 




MuidBia aiuufULui pill 



Section Looking North 



The estimated project costs for this alternative is shown below: 

Base Estimate $213,444 

Sheet Pile Pit $332,680 

Ml Style Building $333,800 

Shield Block Support $30,000 

Bridge Crane $80,000 

Subtotal $989,924 

Overhead and Profit $178,000 

Construction Subtotal $1,168,000 

EDIA $222,000 

Management Reserve $306,000 

Indirect Costs $322,000 

TOTAL $2,018,000 



3 of 11 



Appendix 

The detailed cost breakdown of this alternate is shown below: 

184 
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Alternate No. 2 - Pre-engineered Building 

The alternate investigated the possibility of providing a sheet pile enclosure 185 
similar to Alternate Number 1 while reducing the cost of the above grade building 
by utilizing a pre-engineered building system. The pre-engineered building 
system would utilize manufacturer's standard wall and roof panels. While the 
massing of the above grade structure would be similar, the detailing and 
appearance would not match the MiniBooNE Target Hall. The enclosure plan, 
plan at grade and section would be similar to Alternate Number 1 . 

The estimated project costs for this alternative is shown below: 



Base Estimate $213,444 

Sheet Pile Pit $332,680 

Pre-engineered Building $274,300 

Shield Block Support $30,000 

Bridge Crane $80,000 

Subtotal $930,000 

Overhead and Profit $167,000 

Construction Subtotal $1,097,000 

EDIA $208,000 

Management Reserve $287,000 

Indirect Costs $302,000 




TOTAL $1,894,000 



The majority of the detailed cost breakdown for this alternate can be found in the 
Alternate Number 1 description. The estimate cost breakdown for the pre- 
engineered building system is shown below: 
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Alternate No. 3 - Secant Pile Enclosure 

TheL^aternate investigated the possibility of providing a secant pile enclosure 
similar to Alternate Number 1 and a pre-engineered building system similar to 
Alternate Number 2. The secant pile enclosure consists of a circular ring of cast- 
in-place concrete piers that, after excavation, form the interior wall surface of the 
enclosure. 



fllEtLSI'AlWVA* 
iQf.LOWi 
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SHIELD HXKKS 



Enclosure Plan 



Access to the enclosure would be via steel stairs. The electronics associated 
with the detector would be housed in an above grade electronics caboose. 

This alternate provides a removable metal hatch cover to provide weather 
protection for the enclosure below. The hatch cover would be removed to access 
the precast concrete shield blocks and physics apparatus. 
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Section Looking North 

In order to facilitate final placement of the physics components a 10-ton bridge 
crane would be installed beneath the shield blocks. 

The estimated project costs for this alternative is shown below: 

Base Estimate $213,444 

Secant Pile Pit $590,000 

Hatch Cover $150,000 

Shield Block Support $30,000 

Bridge Crane $80,000 

Subtotal $1,063,000 

Overhead and Profit $191,000 

Construction Subtotal $ 1 ,254,000 

EDIA $238,000 

Management Reserve $328,000 

Indirect Costs $346,000 

TOTAL $2,166,000 
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The majority of the detailed cost breakdown for this alternate can be found in the 
Alter$&te Number 1 description. The estimate cost breakdown for the secant pile 
enclosure and the hatch cover is shown below: 
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Alternate No. 4 - Sheet Pile Enclosure/Hatch Cover 

The alternate investigated the possibility of providing a sheet pile enclosure 189 
similar to Alternate Number 1 and a metal hatch cover similar to that described in 
Alternate Number 3. 




Enclosure Plan 
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Section Looking North 
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The estimated project costs for this alternative is shown below: 

190 

Base Estimate $213,444 

Sheet Pile Pit $332,680 

Hatch Cover $150,000 

Shield Block Support $30,000 

Bridge Crane $80,000 

Subtotal $806,000 

Overhead and Profit $145,000 

Construction Subtotal $951,000 

EDIA $181,000 

Management Reserve $249,000 

Indirect Costs $262,000 

TOTAL $1,643,000 



The detailed cost breakdown of these components can be found in the previous 
alternate descriptions. 
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Alternate Comparison Matrix 

In order to provide cost input to the enclosure and building selection process, an 191 
alternate comparison matrix was developed. Each alternative is comprised of 
different components to assemble a possible design. 





Alt. No. 1 


Alt. No. 2 


Alt No 3 


Alt No. 4 


Alt No 5 


Alt No 6 


Alt. No. 7 


Base Estimate $213,444 


$213,444 


$213,444 


$213,444 


$213,444 


$213,444 


$213,444 


$213,444 


Sheet Pile Pit $332,680 


$332,680 


$332,680 




$332,680 




$332,680 


$332,680 


Secant Pile Pit $590,000 






$590,000 




$590,000 






Ml Style Building $333,800 


$333,800 














Preengineered Building $274,300 




$274,300 






$274,300 






Hatch Cover $150,000 






$150,000 


$150,000 




$150,000 


$150,000 


Shield Block Support $30,000 


$30,000 


$30,000 


$30,000 


$30,000 


$30,000 


$30,000 


$30,000 


Bridge Crane $80,000 


$80,000 


$80,000 


$80,000 


$80,000 


$80,000 






Mobile Crane 
















Gantry Crane $60,000 














$60,000 


Subtotal 


$990,000 


$930,000 


$1,063,000 


$806,000 


$1,188,000 


$726,000 


$786,000 


Overhead and Profit 18% 


$178,000 


$167,000 


$191,000 


$145,000 


$214,000 


$131,000 


$141,000 


Subtotal 


$1,168,000 


$1,097,000 


$1,254,000 


$951,000 


$1,402,000 


$857,000 


$927,000 


EDIA 19% 


$222,000 


$208,000 


$238,000 


$181,000 


$266,000 


$163,000 


$176,000 


Subtotal 


$1,390,000 


$1,305,000 


$1,492,000 


$1,132,000 


$1,668,000 


$1,020,000 


$1,103,000 


Man. Reserve 22% 


$306,000 


$287,000 


$328,000 


$249,000 


$367,000 


$224,000 


$243,000 


Subtotal 


$1,696,000 


$1,592,000 


$1,820,000 


$1,381,000 


$2,035,000 


$1,244,000 


$1,346,000 


Indirect Costs 19% 


$322,000 


$302,000 


$346,000 


$262,000 


$387,000 


$236,000 


$256,000 


TOTALS 


$2,018,000 


$1,894,000 


$2,166,000 


$1,643,000 $2,422,000 $1,480,000 


$1,602,000 



Note: 

While this matrix can provide guidance for general design approach, refinement 
during the subsequent design process will impact the eventual cost. 
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FERMILAB: FESS COST ESTIMATE 





Construction Sub-Contract Subtotal $794,000 
Sub-Contract Overhead and Profit @ 18% $135,000 
Sub-Contract Award Estimate $929,000 




DRAFT 


Project Title: 

FINeSSE Detector 


Project No. 
6-7-59 


Status: 
PDR 


Date: 
11/06/03 


Revision Date: 


ITEM NO. 


DESCRIPTION OF WORK: 


QUANTITY 


UNITS 


UNIT COST 


AMOUNT 


02 


SITE CONSTRUCTION $751,534 












0100 


General 




$27,700 












101 


Mobilize/Survey 






1 


Lot 


$10,000.00 


$10,000.00 






Soil boring 






2 


EA 


$1,000.00 


$2,000.00 






Silt Fence, installed and maintained 






800 


LF 


$14.00 


$11,200.00 






Settleing Basin 






1 


Lot 


$2,000.00 


$2,000.00 






Rock Check Dam 






1 


Lot 


$1,500.00 


$1,500.00 






IEPA Permits 






1 


Yr. 


$1,000.00 


$1,000.00 






















102 


Site Work 




$49,500 














Extend Power from Giesse Rd. 


















Poles 






4 


EA 


$2,000.00 


$8,000.00 






New Crossarms on existing poles 






5 


EA 


$500.00 


$2,500.00 






13.8 Pole mounted Air Switch 






3 


EA 


$1,000.00 


$3,000.00 






Pole mounted area light 






1 


EA 


$1,000.00 


$1,000.00 






Hang Fermi Cable / transformer (1 crew, 3 days) 


2 


Day 


$1,500.00 


$3,000.00 






Comm duct 2-5" 






400 


FT 


$80.00 


$32,000.00 
























Hardstands 




$56,160 














Excavate for Roads (500 x 28 x 1 .5) 






780 


CY 


$8.00 


$6,240.00 






Excavate for Hardstands (47 x 47 x 1.5) 




400 


CY 


$8.00 


$3,200.00 






Haul Excavated Material 






780 


CY 


$2.00 


$1,560.00 






Stone for Road 






1050 


CY 


$30.00 


$31,500.00 






Stone for Hardstand Area 






400 


CY 


$30.00 


$12,000.00 






Guardrail 






1 


Lot 


$1,000.00 


$1,000.00 






Culverts including end sections 






30 


LF 


$22.00 


$660.00 
























Earth Work 




$90,150 














Excavate Hall 






1200 


CY 


$15.00 


$18,000.00 






Excavate Ramp 






1950 


CY 


$15.00 


$29,250.00 






Backfil Ramp 






1950 


CY 


$22.00 


$42,900.00 
























Piping Systems 




$36,700 














Extend LCW piping 






400 


FT 


$80.00 


$32,000.00 






Hydrant + PIV 






1 


EA 


$2,500.00 


$2,500.00 






Ballards 






3 


EA 


$200.00 


$600.00 






Pump discharge and headwall 






1 


EA 


$1,600.00 


$1,600.00 





















Project Title: 

198 FINeSSE Detector 


Project No. 
6-7-59 


Status: 
PDR 


Date: 
11/06/03 


Revision Date: 


ITEM NO. 


DESCRIPTION OF WORK: 


QUANTITY 


UNITS 


UNIT COST 


AMOUNT 






Landscaping 




rt* A O AAA 

$12,000 














Topsoil 






zuuu 


QV 

O Y 


<C/i nn 


c;q nnn nn 






Seeding 






zuuu 


QV 
o Y 


<co nn 


<C/i nnn nn 




















03 




CONCRETE 




$138,900 














Mud Slab 






I u 


\-i Y 


5> 1 ou.uu 


<t-i enn nn 






Floor Slab 






52 


CY 


$300 00 


$15 600 00 






End Wall 






70 


CY 


$450 00 


$31 500 00 






Top Collar 






128 


CY 


$550 00 


$70 400 00 






Crane Runway 






40 


CY 


$250 00 


$10 000 00 






Stairs 






15 


Pjcprc 


$300 00 


$4 500 00 






Stair Enclosure 






15 


CY 


$300.00 


$4 500 00 






Place Fermi Shield Blocks 17 Blks (1 crew +fork + truck) 


2 


Days 


$1 200 00 


$2 400 00 
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P7 97 Pilinn Cnrtpn Stppl 79'y TV 
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55 


Tons 


$1 800 00 


$99 000 00 






VVdlCI o 






9.75 


Tons 


$2 200 00 
$4 000 00 


$21 450 00 






Tpmn Qtn itQ 
1 CI 1 l|J Oil Ulo 






2 


EA 


$8 000 00 
























Mies* MotalQ 

IVI loL- IVlCLdlo 




$110. Q3fi 














Stair 






30 


Pjcprc 


$270 00 


$8 100 00 






1 anrlinnQ 

l_Ctl IU 1 1 IUo 






24 


SF 


$44.00 


$1 056 00 






L_ ICTL-LI Ul ML>o IVI \ 1 a \ 1 1 1 1 iy 






9 


Tons 


$2 200 00 


$19 800 00 






Flprtronir^ Mp77 Cnratinn 






550 


SF 


$14.00 


$7 700 00 






60# Crane Rail for 10 ton crane 






5.4 


Tons 


$1 200 00 


$6 480 00 






40# rail for Hatch 






2.8 


Tons 


$1 000 00 


$2 800 00 






Rails 






1 


Lot 


$5 000 00 


$5 000 00 






Hatch Cover 






1 


Lot 


$60 000 00 


$60 000 00 
























Finishes 




$3,650 














3' x7' Doors 






1 


EA 


$920 00 


$920 00 






Painting Concrete 






1624 


SF 


$0.95 


$1 542 80 






Painting Hatch Cover 






950 


SF 


$1 .25 


$1 187 50 
























Conveying Systems 


















10 ton Gantry Crane 






1 


Lot 


$60 000 00 


$60 000 00 
























Mechanical 




$5,000 














DDC Controls 






5 


Pts 


$1 000 00 


$5 000 00 
























Plumbing 




$15,012 














Drain Tile 






124 


LF 


$12.00 


$1,488.00 






Duplex Sump & basin 






1 


Lot 


$15,000.00 


$15,000.00 
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FINeSSE Detector 


Project No. 
6-7-59 


Status: 
PDR 
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Revision Date: 

199 


ITEM NO. 


DESCRIPTION OF WORK: 


QUANTITY 


UNITS 


UNIT COST 


AMOUNT 






rile r 1 UltrUllUI 1 




\P i 














O r\ t-i t-t \s 1 r\ rf O 1 oi 

opi II IMcl b Z. IgvgIo 






1800 


SF 


$8.00 


$14 400 00 
























C 1 1 

LIcUll lUdl 


















^fou v Trom poie mounieu i rans 






-| 


Lot 


$fi 000 00 


$6 000 00 






tUUn rdl Icl 






-| 


EA 


$3 500 00 


$3 500 00 






99^ Amn Ponal 






-| 


EA 


$4 500 00 


$4 500 00 






*+ r iuui riAiuic 






12 


EA 


$200 00 


$2 400 00 






fZAIL OILJI lo 






4 


EA 


$200 00 


$800 00 












6 


EA 


$175 00 


$1 050 00 






Wolrlinn Oi itlot 






1 


EA 


$600 00 


$600 00 






If A" PnnHi lit W *=; if 1 9 

O/H- wUIIUUIL VV O tt I Z. 






150 


|_F 


$10 00 


$1 500 00 






1 " rnnHi lit W/ft ±£ R 

I L.UMUUIL VV/O tt D 






75 


|_F 


$18 00 


$1 350 00 






,| Rnypc 






5 


EA 


$60 00 


$300 00 






Electric Unit Heater 






3 


EA 


$2 000 00 


$6 000 00 
























Fire Detection 




$14,000 














Spot Smoke Detection 






2 


EA 


$500.00 


$1,000.00 






Fire Panel 






1 


EA 


$10,000.00 


$10,000.00 






Tie To Sitewide FIRUS 






1 


Lot 


$3,000.00 


$3,000.00 







































Station from 8GeV (Q851) to Center of MiniBooNE Detector 

O'Sheg Oshinowo June 25, 2003 



to 
o 
o 





STATION X Y Z 
ft ft ft ft 


COMMENT 


Q851 SQA 


0.00000 99791.76601 97369.55200 715.85533 


Note that Q851 is the location of M851 in the P8 beamline 
and is at the middle of Q851 


DRTGT TARGET 


678.28023 99222.44321 97485.35680 723.00000 


Target, middle 


Station from Q851 to middle of target (Beamsheet) 










CTR of Detector Sphere CTR 


1776.05055 98097.98733 98860.11297 729.22093 




Station from Beamsheet 










2454.33078| 


Station from Q851 to Center of Detector 







